RANDOM WALK AND THE THEORY OF BROWNIAN MOTION* 
MARK KAC,f Cornell University 


1. Introduction. In 1827 an English botanist, Robert Brown, noticed that 
small particles suspended in fluids perform peculiarly erratic movements. This 
phenomenon, which can also be observed in gases, is referred to as Brownian 
motion. Although it soon became clear that Brownian motion is an outward 
manifestation of the molecular motion postulated by the kinetic theory of mat- 
ter, it was not until 1905 that Albert Einstein first advanced a satisfactory 
theory. 

The theory was then considerably generalized and extended by the Polish 
physicist Marjan Smoluchowski, and further important contributions were 
made by Fokker, Planck, Burger, Fiirth, Ornstein, Uhlenbeck, Chandrasekhar, 
Kramers and others [1]. On the purely mathematical side various aspects of 
the theory were analyzed by Wiener, Kolomgoroff, Feller, Lévy, Doob, and 
Fortet [2]. Einstein considered the case of the free particle, that is, a particle 
on which no forces other than those due to the molecules of the surrounding 
medium are acting. His results can be briefly summarized as follows. 

Consider the motion of the projection of the free particlef on a straight line 
which we shall call the x-axis. What one wants is the probability 


x; t)dx 


that at time ¢ the particle will be between x; and x if it was at x» at time ¢=0. 

Einstein was then able to show that the “probability density” P(x | x; t)§ must 

satisfy the partial differential equation 
P 

( ) Ox? ’ 

where D is a certain physical constant. The conditions imposed on P are 


(a) P20 
(2) () f P(xo| x; =1 


(c) lim P(xo | = 0, for x x. 
+0 


* This is an extended version of an address delivered at the annual meeting of the Association 
at Swarthmore, Pennsylvania, December 26-27, 1946. 

t John Simon Guggenheim Memorial Fellow. 

t In what follows we shall identify this projection with the particle itself and hence consider the 
so-called one-dimensional Brownian motion. 

§ The notation P(xo| x; t) and P(n| m; s) for conditional probabilities is that used by Wang and 
Uhlenbeck [1]. It does not conform with the notation adopted in the statistical literature. Had we 
adopted the latter notation we would write P(x; t|xo) and P(m; s|n). 
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Conditions (a) and (b) are the usual ones imposed upon a probability density 
and condition (c) expresses the certainty that at t=0 the particle was at xp. 
It is well known that (1) and (2) imply that 


1 
(2-29)? /4D t 
(3) P(x0| x; t 
and that the solution (3) is unique. 

The greatness of Einstein’s contribution was not, however, solely due to 
the derivation of (1), and hence (3). From the point of view of physical ap- 
plications it was equally, or perhaps even more, important that he was able to 
show that 


(4) 


where R is the universal gas constant, T the absolute temperature, N the 
Avogadro number, and f the friction coefficient. The friction coefficient f, in the 
case the medium is a liquid or a gas at ordinary pressure, can in turn be ex- 
pressed in terms of viscosity and the size of the particle [3]. 

It was relation (4) that made possible the determination of the Avogadro 
number from Brownian motion experiments, an achievement for which Perrin 
was awarded the Nobel prize in 1926. However, the derivation of (4) belongs 
to physics proper, and presents no particular mathematical interest ; we shall 
therefore not be concerned with it in the sequel. 

As soon as the theory for the free particle was established, a natural ques- 
tion arose as to how it should be modified in order to take into account outside 
forces as, for example, gravity. Assuming that the outside force acts in the 
direction of the x-axis and is given by an expression F(x), Smoluchowski has 
shown that (1) should in this case be riggs by 

Two cases of special interest and importance are: 
F(x) = —a; field of constant force (for example, gravity). 
F(x) = —bx; elastically bound particle (for example, small pendulum). 

At this point it must be strongly emphasized that theories based on (1) and 
(5) are only approximate. They are valid only for relatively large ¢ and, in the 
case of an elastically bound particle, only in the overdamped case, that is, when 
the friction coefficient is sufficiently large. These limitations of the theory were 
already recognized by Einstein and Smoluchowski but are often disregarded 
by writers who stress that in Brownian motion the velocity of the particle is 
infinite. This paradoxical conclusion is a result of stretching the theory beyond 
the bounds of its applicability. An improved theory (known as “exact”) was 
advanced by Uhlenbeck and Ornstein and by Kramers. The Uhlenbeck-Orn- 
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stein approach was further elaborated by Chandrasekhar and Doob. 

In what follows we shall be concerned with a discrete approach to the Ein- 
stein-Smoluchowski (approximate) theory. This approach was first suggested 
by Smoluchowski himself; it consists in treating Brownian motion as a discrete 
random walk. Smoluchowski used this approach only in connection with a free 
particle but we shall also treat other classical cases. Moreover, a re-interpreta- 
tion of one of the discrete models will allow us to discuss the important question 
of recurrence and irreversibility in thermodynamics. 

The main advantages of a discrete approach are pedagogical, inasmuch as 
one is able to circumvent various conceptual difficulties inherent to the continu- 
ous approach. It is also not without a purely scientific interest and it is hoped 
that it may suggest various generalizations which will contribute to the develop- 
ment of the Calculus of Probability. 


2. The free particle. Imagine a particle which moves along the x-axis in such 
a way that in each step it can move either A to the right or A to the left, the du- 
ration of each step being r. The fact that we are dealing with a free particle is 
interpreted by assuming that the probabilities of moving to the right or to the 
left are equal, and hence each equal }. Instead of P(xo| x; t) we now consider 
P(nA| mA; ST) =P(n|m; s) which is the probability that the particle is at mA 
at time sr, if at the beginning it was at nA. Noticing that P(n| m; Ss) is alsothe 
probability that after s games of “heads or tails” the gain of a player is vy =m—n, 
we can write 


s! 


0 


if |v| Ss and |»| +siseven, 


otherwise. 


Suppose now that A and 7 approach 0 in such a way that 


(7) —=D, — Xo, St = 


It then follows from the classical Laplace-De Moivre theorem [4] that 


2 
f en "MD id 


8 ;s) = 
(8) lim P(n| m;s) abi 
and hence the fundamental result of Einstein emerges as a consequence of what 
in probability theory we call a “limit theorem.” 

It is both important and instructive to point out a striking formal connec- 
tion between the discrete (random walk) and the continuous (Einstein) ap- 
proaches. We notice that P(n| m; s) satisfies the difference equation 


(9) P(n|m;s +1) = 4P(n| m—1;s) +4P(n| m+ 1;s), 


(6) 
A? 
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which we write in the equivalent form 
P(nd | mA; (s + 1)r) — P(nA | mA; 
A? = (m + 1)A;sr) — 2P(nA | mA; st) + P(nd | (m — 
2r A? 


(10) 


In the limit (7) this difference equation goes over formally into the differential 
equation 


(11) 


which as noted before was the basis of Einstein’s theory. This formal connec- 

tion between the two approaches can be made rigorous, but we shall not go 

into this. However, we shall use it as a guiding heuristic principle in constructing 

models of Brownian motion when outside forces are taken into account. 
Finally, let us mention that it is the relation 


which is responsible for the conclusion that the velocity of a Brownian particle 
is infinite. In fact, in our model, the ratio A/r plays the role of the instantaneous 
velocity and it obviously approaches infinity as A-0. 


3. Particle in a field of constant force and in the presence of a reflecting 
barrier. We again consider random walk along the x-axis in which a particle can 
move A to the right or A to the left, the duration of each step being r. We now 
introduce the following new assumptions: 

(a) The probability of a move to the right is g=}— A, and consequently the 
probability of a move to the left is p=}+ A. Here, is a certain physical con- 
stant, and A must be chosen sufficiently small so that g>0. 

(b) When the particle reaches the point x =0 (reflecting barrier) it must, in 
the next step, move A to the right. 

Without the assumption (b) the problem would be quite simple and of no 
great physical interest. In actual experiments with heavy Brownian particles, 
like those of Perrin, the bottom of the container acts as a reflecting barrier and 
the elucidation of its influence on the Brownian motion is of considerable theo- 
retical interest. 

This problem has been solved by Smoluchowski, on the basis of his equa- 
tion (5) but we shall show that one can solve the discrete problem and obtain 
Smoluchowski’s result by passing to a limit. 

Assuming that the particle starts from nA 20 (n an integer) we seek an ex- 
plicit expression for P(n|m; s). We first notice that P(n|m; s) satisfies, for 
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mz=2, the difference equation 


(12) P(n| m;s +1) = gP(n|m—1;s) + pP(n| m+ 1;s), 
and that for m =1 and m=0 we have 

(12a) P(n| 1; 5+ 1) = P(n|0;s) + pP(n| 2; 5), 
(12b) P(n|0;s +1) = pP(n| 1; 5). 


We also have the initial condition 
(13) P(n| m; 0) = 6(m, n), 


where 5(m, m) denotes, as usual, the Kronecker delta. 

The difference equation (12) when rewritten in the form analogous to (10) 
can be shown to go over formally (in the limit A-0, r-0, A?/2r7 =D, nA—x., 
mA—x, st =t) into the differential equation 


oP oP 
(14) 


which is of the form (5) with F(x) = —a= —48Df. 

To find P(n| m; s) we use a method which is basic in the study of the so-called 
Markoff chains, of which our problem is but a particular example, and which in 
its essentials goes back to Poincaré [5]. Let (p), be the (infinite) vector 


P(n| 0; s)7 
P(n| 1; s) 
(15) = | P(n| 2; s) 
and A the infinite matrix 
07000 
10 p00 
(16) A=|0 q¢ 0 0 
00q 0 


Then, the difference equation (12), (12a) and (125), can be written in the matrix 
form as 


(17) (p)s41 A(p)s. 
Thus it follows immediately that 


(18) = A*(D)o, 


4 
fe 
on 
a 
. 
5 
j 
A 
= 
5 
- 


374 RANDOM WALK AND THE THEORY OF BROWNIAN MOTION [September, 


where (p)o is the vector 


1 ’ 


1 being the mth component, the components being numbered from zero on. In- 
terpreting (18), we see that 


(19) P(n| m; 5s) = the (m, n) element of A‘. 

To make use of (19), we notice that if R>n+s and we consider the finite matrix 
Apr, which is the upper left R by R submatrix of A, then form<R 

(20) the (m, n) element of A’ = the (m, n) element of rm 

or equivalently, 


(21) the (m, 1) element of A’ = lim of the (m, n) element of Ap. 


For each R there exist matrices Pz and Qpr such that 


(22) POe = I 
and 
do(R) 0 
(23) Ar = Pr wise Or, 
0 Ar-1(R) 
Ao(R), A(R), being the eigenvalues of the matrix Ar. To simplify 


the notation we write \; for \;(R). 
Multiplying the matrix Ar s times by itself and making use of (22), we ob- 
tain 


= 
Xo 0 
(24) Ar = Pr oh Or 
= 0 


and one can calculate the (m, m) element of A’ explicitly provided the diag- 
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onalization (23) can be performed explicitly. This indeed is the case. 


Let (x)o, (x)1,° °°, (x)ra be the “right” and (y)o, (y)1,---, (y)r-1 the 
“left” eigenvectors belonging respectively to the eigenvalues Xo, Ay, , 
In other words, for k=0, 1, -- - , R—1, we have 


Ar(x)e = 
Ar(y)e = d(y) 


where Ap is the transpose of Ar. 
Suppose furthermore that the vectors can be so normalized that 


(25) (x)e-(y)e = 1, k=0,1,---, 


where (x)x: (y)% is the inner (dot) product of the vectors. Since it is well known 
that 


(x)j-(y)e = 0 for \; Ax 


we see that, in the case when all the eigenvalues are distinct, we can take as Pr 
the matrix whose columns are the vectors (x), and for Qz the matrix whose rows 
are the vectors (y)x. 

In order to determine the eigenvalues and the right eigenvectors we consider 
the system of linear equations 


px1 = Axo 
Xo + px, = Ax 
(26) qx1 + = 
= Axr-1, 
and the extended infinite system 
px1 = AXo 
xo + = AX 


(27) 


qxr-1 + = AXR 


If we can find non-trivial solutions of (27), for which 
(28) tr = 0, 


we will have found solutions of (26). 

It turns out that (28) will yield an equation in \ which has R distinct roots 
and thus our procedure gives us all eigenvalues, and consequently all right 
eigenvectors. Multiplying the members of the equations of (27) by 1, 


| 
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z, 2*,..., and adding, we obtain formally 
1 1 0 


or, upon introducing the abbreviation, 


(29) fle) = ast, 
we have 
(30) + 20] + 29] = M0). 
From (30) we obtain 


and since this function is analytic in the neighborhood of zero the formal pro- 
cedure used above can be justified. 
Let p; and pz be the reciprocals of the roots of 


(32) Az +p =0. 
We have then 

{ 1— dz 
p(1 — prz)(1 — 2s) 


and introducing partial fractions, 


1— dz 1 A—op1 1 1 pe—X 1 


(34) i 
— piz)(1— p22) p 1 — pez 
Thus 
(35) — for k= 1, 


and, in particular, the equation xg =0 assumes the form 


A- 
(36) p2 R 


Equation (36) must now be solved for \. Assuming R to be even, and seeking 
solutions in the form 


= pq cos 0, 


‘ 
2 ’ 
q — Az + 
| 
- 
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we are led to the equation 
tan RO 1 


tan 


For R>(2p—1)" this equation is seen to have R—2 distinct roots, 
Or-2, which lie in the subintervals 


where j ranges through the integers from 1 to R—1 with the exception of j = R/2. 
Corresponding to @:, @2, Or-2 we have R—2 distinct eigenvalues, 


= 24/pq cos Ox, k=1,2,---,R—2, 


and the components of the right eigenvector belonging to \, can be written in 


the form 
m/2 


sin 


where 


q 
if 
= p 
2 q if ao 0, 
and a, is a normalizing constant which will be fixed later. For sufficiently large 


R the remaining eigenvalues Xo and Ag_1 can be shown to be given by the 
formulas 


= cosh 0, = — ro, 
where 6 is the only positive root of the equation 
tanh RO 1 
tanh 6 2p - 


The components of the corresponding right eigenvectors are given by the ex- 


pressions 
(m) q ( sinh 
x ao | — cos 


(m) 
XR-1 


sinh 
adr_1(— 1) (cosh — 2BA 


sinh 60 


It remains now to find the left eigenvectors. This can be accomplished in exactly 
the same manner and we merely quote the results. We obtain 


q 
q 
X 
k 
| 
j 
|_| = 
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m/2 © 
= (2) (cosmo — 2BA 
q 


sin 


for m=0, 1,..., R-1;Rk=1,2,---, R—2, and 


m/2 sinh mé 
yo bo (2) (cosh — 2BA ‘) 
q 


sinh 69 
(m) af sinh mo 
1) (£) (cosh mA = 2BA ). 

q sinh 4 
To satisfy the normalization conditions (25) we must have 

R-1 , 
(37) + = 1, k=1,2,---,R—-2, 

m=1 

R-1 2 
(38) aabi(a + ¥ Fn(6o) ) = 1, k=0,R—1, 

m=1 
where 

sin 
= cos mO — 28A — 
sin 
and 
sinh 
F,,(6) = cosh m@ — — . 
sinh 0 

We can, of course, put d)=a1= - - - =a@r_1=1, and determine the b’s from (37) 


and (38). Referring back to (19), (20), and (24), and recalling that columns of 


Pp are the right eigenvectors (x);, and the rows of Qe are the left eigenvectors 
(y)x, we obtain 


s (m) (n) 
(39) P(n| m;s) = axe 


k=0 
or, more explicitly, 


n/2 m/2 
P(n| m; s) = bo(24/pq cosh (4); F n(00)F n(80) [1 + (— 1) 
n/2 q m/2 — R--2 


Making use of (21), we can achieve considerable simplification by letting 
R-~, In fact, it can be shown that 


(41) = [1 + (= 1) 


2pq 
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+ T a * (2 pq J cos + 


Although in various places we have tacitly assumed that p and q are different 
from 3, the final formula (41) can easily be seen to be valid also for the case 
p=q=}. In this case (free particle in the presence of a reflecting barrier) the for- 
mula assumes the remarkably simple form 


2 
(42) P(n| m;s) = — f cos* © cos mO cos 
0 


and the right member can be expressed in terms of binomial coefficients. This 
formula can also be derived in a much simpler way using, for instance, the classi- 
cal method of images. 

In the limit 


one can show that 


lm P(n| m; s) -f P( xo | x; t)dx, 


Ze 


where 


2 
43) P = 4 f Dit y)dy, 
(43) P(xo| x; = 4Be-#* + e y 


and 
g(x, y) = cos xy — 28(sin xy)/y. 


The proof of this theorem is not elementary inasmuch as it utilizes the so called 
“continuity theorem for Fourier-Stieltjes transforms” [6]. Formula (43) can be 
shown to be equivalent with Smoluchowski’s formula given in [1]. 


4. An elastically bound particle. Again the particle can move either A to the 
right or A to the left, and the duration of each step is r. However, the probability 
of moving in either direction depends on the position of the particle. More pre- 
cisely, if the particle is at kA the probabilities of moving right or left are 


—{1-——) or — —}, 
2 R 2 R 
respectively. R is a certain integer, and possible positions of the particle are 


limited by the condition —RSkSR. The basic probabilities P(n|m; s) now 
satisfy the difference equation 


| 
4 
A? 
A—0, r— 0, = D, nh — Xo, st = 1, 
T 
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R 1 R-—-m+1 
(44) P(n| m;s + 1) = m+ 1; s) 1;s), 
which must be solved with the initial condition 
(45) P(n| m; 0) = 5(m, n). 
In the limit 
A? 
2r Rr 


St = 1, — Xo, mA — x, 


the difference equation (44) is seen to go over formally into the differential 
equation 


oP 0(xP) 
(46) D— 
ot Ox Ox? 


which is Smoluchowski’s equation (5) with F(x) = —x/yf. 

The discrete problem in a different form and in a different connection was 
first proposed and discussed by P. and T. Ehrenfest in 1907 [7]. In the next 
section we shall discuss their original formulation. A fairly detailed treatment 
was given by Schrédinger and Kohlrausch in 1926 [8] and a brief exposition can 
be found in the review article of Wang and Uhlenbeck [1]. It seems that Schré- 
dinger and Kohlrausch were the first to point out the connection between the 
Ehrenfest model and Brownian motion of an elastically bound particle. How- 
ever, an explicit solution of (44) with the initial condition (45) was apparently 
not known. I have recently found such a solution using the matrix method de- 
scribed in Section 3 [9]. Instead of the infinite matrix of that section we must 
now consider the finite matrix 


1 
2R 
2 
1 = 
2R 
(47) B= 1 3 
01-— 0 — 
2R 2R 
1 
0 0 — 0 
2k) 


and the problem is again reduced to finding the eigenvalues 
Xo, Of B and matrices P and Q such that 


4 
; 
4 
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ae 


1947] RANDOM WALK AND THE THEORY OF BROWNIAN MOTION 
(48) 
and 
0 
(49) B=P Q. 
0 
AR) 


As before, P(n| m; s) is the (m, n) element of B*, where 


) 
0 
(50) 
0 


8 
\ AR) 


ty 
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In order to perform the diagonalization (49) explicitly we start (following the 
procedure of Section 3) by trying to find the eigenvalues and the right eigenvec- 


tors of B. For this purpose we consider the system of linear equations 


: 
— 


1 1 3 
(1-5) nt 


and the auxiliary infinite system 


(52) 


2 
xo + — = Ax} 
2R 
1 
oR = AXoR, 
1 
—— = AXo 
2R 
= 2R 
| 
“A 
2R 2R 3 
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1 2R+1 
2R 2R 
2R+2 
2R 


Xor+1 = AXor 


= AXor+1 


If we can find non-trivial solutions of (52) for which 
(53) Xer+1 = 


we will have found solutions of (51). It will turn out that this procedure will 
again yield all eigenvalues and right eigenvectors. Multiplying the members of 
the equations of (52) by 1, z, 2?, - - - , and adding, we obtain formally 


k 


k=0 k=0 k=0 


or, introducing the abbreviation 


fe) = 


k=0 


2 1 


We thus get the differential equation 


A-2 
(54) f'(s) = 2R f(z), 
1 — 2? 
whose solution satisfying f(0) =x» is easily found to be 
(55) f(z) = xo(1 — + 
Since f(z) is analytic in the neighborhood of z=0 the formal procedure can be 
justified. 
We now notice that if 


f(z) is a polynomial of degree 2R, and hence x2x41;=0. The numbers (56) are 
thus seen to be eigenvalues of B and, since there are 2R+1 of them, we see that 
we have found all the eigenvalues. It also follows that the components of the 
right eigenvector belonging to the eigenvalue \;=7/R can be taken as 


(i) (i) (7) (i) 
Co =1,Ci , C2 Cor, 
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where the C’s are defined by the identity 


(57) (1—s) = Cys. 


kad 

So far we have followed very closely the procedure described in Section 3. Sur- 

prisingly enough, we encounter unexpected difficulties in trying to carry out 

the analogy still further and determine by similar means the left eigenvectors. 

To find the matrix Q we resort to a different method. Let us first recall that 

P can be taken as the matrix whose jth column (for convenience columns and 
rows are numbered from —R to R) is 

/ 


1 

Cr. 
Matrix Q must satisfy the equation 

PQ? =I, 
which is an immediate consequence of the equation PQ=I, and hence denoting 
by a_r, -**,@o, °° +, @g, the consecutive elements of the jth column of Q’, we 
must have 
(58) Cri = 5(j, r), R, R. 


From (58) it follows that 


R R R 
Rte R+ (k) (k) R+ 
=D = Cryo = Cre 


r=—R r=—R k=—R k=—R Finca 
R 2R 
(k) 
las, 
k=—R s=0 


or, by virtue of (57), 


k=—R 1=0 


Thus 


(1 —2)® 2)2k 1=0 
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Let 
so that 
and 1—z=——- 
+f 

In terms of £ (59) assumes the form 

(— 2R 

92k 


and since by (57) 


—i) 


‘ 2R 
049 Der, 


we obtain, by comparing coefficients of corresponding powers of ¢, 


(— 1) oo? 


a = 


or finally, 


(— 
(61) a, = 


Formula (61) determines explicitly the elements of Q’ (and hence of Q), and it is 
now possible to write an explicit expression for P(n|m; s). In fact, making use 


of (50), we obtain 
(62) P(n|m;s) = Cr+iCr+m 


In the limit 


1 
2r Rr 
we have 
lim P(n| m; s) -f P(xo| x; t)dx, 
22 
where 
(63) P(xo| x; = aye 
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The proof is again made to depend on the continuity theorm for Fourier-Stieltjes 
transforms. 

The frequency function (63) was first discovered by Lord Raleigh [10]. Its 
connection with Brownian motion of an elastically bound particle, in the strongly 
overdamped case, was established by Smoluchowski who arrived at it quite 
independently. 


5. The Ehrenfest model. Irreversibility and recurrence. Imagine 2R balls 
numbered consecutively from 1 to 2R, distributed in two boxes (I and II) so 
that at the beginning there are R+n, —RSnSR, balls in box I. We chose at 
random an integer between 1 and 2R (all these integers are assumed to be equi- 
probable) and move the ball, whose number has been drawn from the box in 
which it is, to the other box. This process is then repeated s times and we ask 
for the probability Q(R+n|R+m; s) that after s drawings there should be 
R+m balls in box I. 

A moment’s reflection will persuade one that this formulation (originally 
proposed by P. and T. Ehrenfest) [7] is equivalent to the random walk formu- 
lation of Section 4, if one interprets the excess over R of balls in box I as the dis- 
placement of the particle (A=1). Thus 


QR+n|R+m;)= P(n| m;s), 


where P(n| m; s) has the meaning of Section 4. 

In the present formulation we have a simple and convenient model of heat 
exchange between two isolated bodies of unequal temperatures. The tempera- 
tures are symbolized by the numbers of balls in the boxes and the heat exchange 
is not an orderly process, as in classical thermodynamics, but arandom one like 
in the kinetic theory of matter. The realistic value of the model is greatly en- 
hanced by the fact that the average excess over R of the number of balls in box 
I, namely, the quantity 


R 


> mP(n| m; s) 


m=—R 


can easily be shown to be equal to 


which in the limit R->«, 1/Rr—y, st =t, gives 


or the Newton law of cooling. 

There are several proofs of (64) [11]. The most straightforward one, which 
is not however the simplest, is based on formula (62). 

The Ehrenfest model is also particularly suited for the discussion of a famous 
paradox which at the turn of this century nearly wrecked Boltzmann’s inspired 


5 
& 
4 
x 
= 
ry 
tt 
a 
7 
i 
: 


386 RANDOM WALK AND THE THEORY OF BROWNIAN MOTION [September, 


efforts to explain thermodynamics on the basis of kinetic theory. In classical 
thermodynamics the process of heat exchange of two isolated bodies of unequal 
temperatures is irreversible. On the other hand, if the bodies are treated as a 
dynamical system the famed “Wiederkehrsatz” of Poincaré asserts that “almost 
every” state (except for a set of states which, when interpreted as points in 
phase space, form a set of Lebesgue measure 0) of the system will be, to an arbi- 
trarily prescribed degree of accuracy, again approximately achieved. Thus, 
argued Zermelo, the irreversibility postulated in thermodynamics and the “re- 
currence” properties of dynamical systems are irreconcilable. Boltzmann then 
replied that the “Poincaré cycles” (time intervals after which states “nearly 
recur” for the first time,—the word “nearly” requiring further specification) are 
so long compared to time intervals involved in ordinary experiences that pre- 
dictions based on classical thermodynamics can be fully trusted. This explana- 
tion, though correct in principle, was set forth in a manner which was not quite 
convincing and the controversy raged on. It was mainly through the efforts of 
Ehrenfest and Smoluchowski that the situation became completely clarified, 
and the irreversibility interpreted in a proper statistical manner. 

It will now be easy to discuss this explanation by appealing to the Ehrenfest 
model. Let P’(n|m; s) denote the probability that after s drawings (the dura- 
tion of each drawing is r) R+m balls will be observed for the first time in box I 
if there were R+2 balls in that box at the beginning. In particular, P’ (n| n; s) is 
the probability that the recurrence time of the state “n” (defined by the presence 
of R+ balls in box I) is sr. One can then show that 


(65) P'(n| n;s) = 1, 


s=1 


or, in other words: each state is bound to recur with probability 1. This is the statis- 
tical analogue of the “Wiederkehrsatz.” One can show furthermore that the 
mean recurrence time, namely, the quantity 


>, stP'(n| n; S) 


s=1 


is equal to 


(R + n)\(R —n)! 
(2R)! 


This is the statistical analogue of a “Poincaré cycle,” and it tells us, roughly 
speaking, how long, on the average, one will have to wait for the state “n” to 
recur. 

If R+mn and R—n differ considerably, 8, is enormous. For example, if 
R=10000, = 10000, 7 =1 second, we get 


(66) 


6 = 270000 seconds (of the order of 10° years!). 
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If on the other hand, R+-n and R—x are nearly equal, 0, is quite short. If in the 
above example we set » =0 we get (using Stirling’s formula) 


~ 100\/x seconds ~ 175 seconds. 


It was Smoluchowski who advanced the rule [12] that if one starts in a state 
with a long recurrence time the process will appear as irreversible. In our ex- 
ample if one starts with 20000 balls in one box and none in the other, one should 
observe, for a long time, an essentially irreversible flow of balls. On the other 
hand, if the mean recurrence time is short, there is no sense to speak about ir- 
reversibility. 

We now give the proofs of (65) and (66). We shall base our considerations 
on a formula which Professor Uhlenbeck used for similar purposes in some of his 


unpublished notes. The formula in question is: 
a—1 
(67) P(n| m; s) = P'(n| + P'(n| m; k)P(m| m; s — k). 
k=1 


To convince oneself of the validity of this formula we divide all possible ways 
of reaching “m” from “n” in s steps into classes according to when “m” has been 
reached for the first time. We then observe that starting from “n” one can reach 
“m” in s steps in the following s mutually exclusive ways: 

(1) “m” is reached for the first time after s steps. 

(2) “m” is reached for the first time in 1 step and then, starting from “m” it is 
again reached in s—1 steps. 

(3) “m” is reached for the first time in 2 steps and then, starting from “m”, 
it is reached again in s—2 steps, and so forth. We note furthermore that 
the probability that “m” will be reached for the first time in k steps 
and then, starting from “m,” it will be reached again in s—k steps, is 

(68) P'(n| m; k)P(m|m;s — k). 
This completes the proof of (67). 

It should be emphasized that the justification of using the product of prob- 
abilities in (68) rests upon the fact that in our process the past is independent of 
the future. In other words, once we know that the system starts, say, from “m,” 
its subsequent behavior is independent of the way in which “m” was reached in 


the first place. 
We introduce now the generating functions 


(69) h(n | m;z) = P(n| m; s)z* 
(70) m; 2) m; 


and note that (67) is equivalent to 
h(n| m; 2) = g(n| m; 2) + h(m| m; 2)g(n| m; 2), 


# 
H 
34 
Mt 
i 
$ 
ag 
a 
13 
| | 


388 RANDOM WALK AND THE THEORY OF BROWNIAN MOTION ([September, 


or 
h(n | m; 2) 


In particular, 


(72) 
1+ h(n| 1+ h(n| n;2) 
and we also note that 
dh(n| Nn; 2) 
dg(n| n; 2) dz 
7 = 
(73) dz (1+ h(n| n; 


From the definition of g(n| n; 2), we obtain 
(74) lim g(m|n;z) = > P'(n| n; S) 


dg(n|n;z) 


(75) lim = srP'(n| n; 5). 


21 dz 


It is from these formulas that we shall derive (65) and (66). We have, using (62) 


i=—R e=1 
and since 
(— 1)" (-n) 
1 = ———_ CrijCrin 
R+jUR+ 
we obtain 
— R 1 = 
(76) 1+ h(n|n;2) = 
j 
1——z 
R 


All terms in the sum on the right hand side of (76) are regular at z=1 except the 
term corresponding to j= R, which has a simple pole at that point. Thus we can 
write 


1 
1+ A(n|n;2) = pe) = 


where #(z) is regular at z=1. We see that 
lim (1 + h(n| n;2)) = © 


| 
s=1 
be 
, 
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and hence, using (72) and (74) 


P'(n| = 1. 


e=1 
It is easy to see that 


(— A (2R)! 


and, denoting this expression by w, we have (for |s| >1) 
dg(n|n;s) — +o 
dz [(1 — 2)p(z) + 


and hence 


dg(n| n; 2) 


21 dz 


This together with (75) yields (66). 

The above considerations can be extended to more general processes. How- 
ever, Markoffian processes (i.e., processes for which (68) is valid) are still the 
only ones for which one can also calculate the “fluctuation” of the recurrence 
time, namely, the quantity 


1 
w 


(77) (n| n;s) — On. 


Without going into the details, let us mention that (77) can be calculated in 
terms of 


d2g(n| n; 2) 
lim ——————_ 
dz? 


The fluctuation (77) gives us a measure of stability of the mean recurrence time 
inasmuch as it permits us to estimate how likely (or unlikely) it is to get a spec- 
ified deviation of the actual recurrence time from the mean. It may seem that 
since the generating function g(m|; z) is known explictly it should be easy to 
get an explicit expression for P’(n|n; s). This, however, is not the case. We have 
not succeeded in finding such an expression, except for P’(0|0;s), and even then 
we had to use a different method. We shall give a brief description of this 
method. Let 


P(n| m;1) = pam 
Then, 
P'(n| nN; s) = Pam,P 


va 
| | 
| 
e=1 
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where the accent on the summation sign indicates that ms¥n,j=1,2,--+,s—1. 


Now let 
0 ifi=n 
ifi Xn. 
Noticing that 


we can write 


where the summation is now extended over all m,. If B is the matrix 


((Pnm)) 
and B, the matrix 
((EnPnm€m))» 
we see that 
(78) P'(n| n; s) = (n, n) element of BB; B. 


We may note that B, is obtained from B by crossing out the mth row and the 
nth column of the latter, and replacing them by a row and column consisting 


entirely of zeros. If B, can be explicitly diagonalized, that is, written in the 
form 


M1 0 


where 
P. 101 = 7 


one can calculate P'(n| n; Ss), explicitly using (78). 

We have applied this method to the Ehrenfest model, but only in the case 
when the middle (zeroth) row and column of B are replaced by a row and col- 
umn consisting entirely of zeros have we been able to diagonalize explicitly the 
resulting matrix B,. The diagonalization proceeds very much as in Section 4, 
but it has proved necessary to distinguish between the cases when R is even or 
odd. In case R is even, we were able to derive the formula 


| R 1 
(79) P'(0| 0;s) = (2) s2=2, 


R 


2 
= Ei, 
= 
4 
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where the summation is extended over all odd integers 7 between —R and R. 
The details of the derivation are somewhat tedious and will not be reproduced 
here. Formula (79) furnishes a partial solution to a question left open by Wang 


and Uhlenbeck [1]. 
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VIBRATION MODES OF TAPERED BEAMS 
EDMUND PINNEY, University of California 


The vibration modes of harmonically oscillating thin beams are obtained 
by solving the differential equation 


(1) — = 0 


subject to certain boundary conditions. In this differential equation, x is the 
distance along the beam, y is the amplitude of the transverse vibration, E is 
Young’s modulus for the beam material, J is the cross-sectional moment of 
inertia of the beam, p is its linear density, and w is the circular frequency of vi- 
bration. 
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When EI and p are constant along the beam the solution of (1) is simple, and 
may be expressed in terms of trigonometric and hyperbolic functions. The fre- 
quency w will be a root of a simple transcendental equation. Such theories are to 
be found in most textbooks on vibration theory. 

When EI and p vary very slowly with x, solutions to (1) may be obtained 
by writing EI = Eolo[1+&(x)], p=po[1+n(x)], y=yo(x)+$(x), where and 
fo are constants, where yo(x) is a linear combination of 


(2) sh (kx), ch (kx), sin (kx), cos (kx), 
where the quantities £, , ¢/yo are much smaller than unity, and where 
(3) = pow*/Eolo. 


Substituting into (1) and neglecting terms of higher order than the first in &, 
¢, we obtain 


(4) d*t/dx* — = f(x), 
where 
(S) f(x) = — + k4yon. 


Equation (4) has a particular solution 
(6) t(x) = ve {10 [sh k(x — #) — sin k(x — 2) Jdt. 


The general solution is obtained by adding to this an arbitrary linear combina- 
tion of the functions in (2). 

If EI and p may be represented by a finite number of terms of a Fourier se- 
ries, then y may be obtained from (1) as a Fourier series whose coefficients are 
the solutions of certain algebraic equations [1; §VII]. This method is useful for 
beams of greater taper than those to which the analysis of the preceding para- 
graph may be applied. Another method of great practical importance is due to 
Myklestad [2; Ch. VI]. 

When EZ and p are certain particular functions of x, it is possible to obtain 
exact solutions of (1) in terms of known functions. One case of considerable in- 
terest is that in which 


(7) EI = EoIox**?, 


where EoIo and po are constants. The case where B=0 leads to solutions of the 
form x? where p(p—1)(p+n)(p+u—1) =k*. If one of the four conditions 


holds, the vibration modes of the beam may be expressed in terms of Bessel 
functions. 


To show this, we insert (7) into (1), getting 


(8) 6=1/2, B=-1/2, B=n/2, 
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(9) x?(d?/d x?) — = 0. 
If we write 

d 

= 

dx 

it follows that 
a? 
x? —— = — 1) 


and that 
O(x"f) = n)f. 

Using these, we may write (9) in the form 

By [3; §4.5], the equation 
(11) (0+ + a — 28)(6 + a — 26r)(0 + a — 28 — = 
has the solution 
(12) y = + BY,(cx®) + CI,(cx*) + DK,(cx')]. 


Comparing (10) and (11), we see that solutions may be obtained in the following 
cases: 


Type v B 
I 1/2 2k 
(13) II —p —1/2 p/2+1 2k 
III 1/p 1/2 


Beams of type J are of particular interest because the exponent for EI is 2 “a 
greater than that for p, a situation often encountered in practice. Examples are 
the wedge, corresponding to y=1, and the cone, corresponding to 4=2. For 
modes of type I 


(14) » y = F,(A, B,C, D; 2), 
where 
F,(A, B, C, D; x) 
= + BY,(2kx"!*) + CI,(2kx'!?) + 
By [3; §§3.2, 3.56, 3.71], it is known that 
(16) F{ (A, B, C, D; x) = kF,-1(A, B, C, —D, x), 


(15) 


a 
i 
i 
4 
| 


394 VIBRATION MODES OF TAPERED BEAMS [September, 


and 
(17) Fyyi(A, B, C, D; x) + xF,-1(A, B, —C, —D; x) = (u/k)F,(A, B, —C, D; x). 


It follows from (16) that the vibration frequencies for a free-free beam (end 
conditions: y’’ =0, y’’’ =0) of type J are the same as those of a similar clamped- 
clamped beam (end conditions: y=0, y’ =0) of type I in which yp is replaced by 
p—2. 

For beams of types I, II, III, and IV, (12) provides the basis for the solution 
of a large number of beam vibration problems. The determination of the natural 
period of vibration of such a beam usually leads to a fourth order determinant 
whose elements are Bessel functions. The roots can be calculated numerically 
without very great difficulty if the Bessel functions are tabluated. Forced vibra- 
tion problems of such beams are similar in theory and complexity. 

In the preliminary design of aircraft this theory may be useful in obtaining 
approximate vibration modes of airplane wings to be used in preliminary flutter 
calculations. Many wings have tip chords that are small in comparison with 
their root chords. Such winds may, with little error, be considered to taper to a 
point at a distance a, say, from the root of the wing. The wing should be approxi- 
mately of type I. 

To avoid singularities, take B=D=0 in (15). Also y=0, y’=0 when x=a 
for a wing that may be considered to be cantilevered at the root. From (15) and 
(16), we have 


(18) F,(A, 0, C, 0; a) — 0, F,1(A, 0, Cc, 0; a) = 0. 

By (3), (15), if 21,2, +++, are roots of 

(19) = 

then the vibration frequencies, a, we, - + + , are given by 

(20) wn = 


The ratio C/A may then be computed from (18). 

For wings that cannot be considered to taper to a point, this determination 
may still be useful in providing a first approximation for the frequency calibra- 
tion. 

By (13), the vibration modes may be expressed in terms of elementary func- 
tions when p=n+4 for beams of types I and II, and when p=  iaatl for 
beams of types III and IV, 2 being an integer. 
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TETRAHEDRONS HAVING A COMMON FACE* 
VICTOR THEBAULT, Tennie, Sarthe, France 


The main theorem of this paper is: Consider two tetrahedrons T=ABCD, 
T'=A'BCD, having a common face BCD. Let the perpendiculars at A to faces 
CDA, DAB, ABC of T meet the planes of the faces A'CD, A'DB, A'BC of T' in 
Ay, Az, As, and the perpendiculars at A’ to faces A'CD, A’DB, A'BC of T’ meet 
the planes of faces CDA, DAB, ABC of T in A’, A's, A’3. Then the perpendiculars 
from A to plane A,A2A3 and from A’ to plane A';A',A's; meet in a point in the 
common face BCD. 

If we designate by 6, y, 6 the angles made by planes A’CD, A’DB, A'’BC 
with the face BCD, and further designate by a and a’, b and b’, ¢ and c’ the di- 
hedral angles having for edges BC and DA, - - - , of T, the equations of the per- 
pendicular from the vertex A to the plane A1A2A3 are: 


y t 


sin c’ cot (c’ — sin cot (b’ — y) sin a cot (a — 4) 


This perpendicular meets the plane of face BCD in a point A” whose normal co- 
érdinates with respect to the triangle BCD are: 


(y’, 2’, ’) ~ cot (c’ — B), cot (b’ — y), cot (a — 4). 


Because of symmetry, the perpendicular from A’ to the plane A’;A’2,A’s 
meets the plane of face BCD in the same point A’’. Hence the theorem. 

If the tetrahedron T is orthocentric and A’=H, the feet of the altitudes of 
triangle BCD being Bi, Ci, Di, then 


U 


y 4 
tan B,AB tan CyAC tan 


But in the triangle DAD, 
DD, y’ 2! 


and 
tan D,AD BB, CC, DD 


Hence A’’(y’, 2’, t’) coincides with the centroid of BCD. This provides the solu- 
tion for problem 4150 [1945, 102]. 

Our theorem also answers the following question which we have raised 
(Annales de la Société Scientifique de Bruxelles, 1925, p. 302) and which has here- 
tofore remained unanswered. Two given tetrahedorns T=ABCD and T"’ 
=A’B’C’D’ are homological. The perpendiculars at A to the faces CDA, DAB, 
ABC of T meet the planes of faces A’C’D’, A’D’B’, A’B’C’ of T’ in Ai, As, As, 
and the perpendiculars at A’ to the faces A’C’D’, A’D’B’, A’B’'C’ of T’ meet 
the planes of faces CDA, DAB, ABC of T in A’;, A’z, A’3. In the same way we 
obtain, relative to the vertices B, C, D of T and B’, C’, D’ of T’, the triads of 


* Translation from the French by J. B. Lennes, University of Oklahoma. 
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points (B,, Bs, B;) and B’;, B’s, B’s), - - -. The perpendiculars from A, B, C, D 
to the planes A1A2A3, BiB2Bs, - - - , meet in a point O, and the perpendiculars 
from A’, B’, C’, D’ to the planes A’;A’2A’3, B’:B’2B’s, - - + , meet in a point O’. 
Are O and O’ corresponding points in the homology (ABCD, A'B'C’D’)? 

The answer is affirmative, for it is sufficient to show that lines AO and A’O’, 
BO and B’0O’, CO and C’0’, DO and D'O’ meet the plane of homology (P) inthe 
same points, as shown by the preceding theorem. It can be added that the point 
O is on the perpendicular to plane (P) through the isogonal conjugate O, of the 
center of homology with respect to the tetrahedron 7, and that O’ is on the per- 
pendicular to (P) through the isogonal conjugate O’ of the center of homology 
with respect to T’. 

In the general case first treated above, the perpendicular AX from a point 
P=A’' whose normal coérdinates in T are (x1, y1, 21, 41) to the plane AiA2A; has 
the following equations: 


y z t 


cos c’ 1 cos 1 cos a 1 


and meets the plane of face BCD in a point P’ whose coordinates are: 


y 3V 


A, B, C, D being the areas of the faces and V the volume of T. 
(a) If P is on Cayley’s cubic surface 


A B Cc D 
—+—+—+—=0, 
t 


the line AX =A P’ is parallel to plane BCD which is normal to plane A,A2A3. 
In particular, this is the case when P coincides with the midpoint of one of the 
28 line-segments joining the centers of the spheres tangent to the four planes of 
the faces of T. 

(b) The equations of the perpendiculars BY, CZ, DT from B, C, D to the 
planes B,B2B3, CiC2C3, D:D2D;3 (paragraph 2) are similar in form to that of 
AX. The following relations are necessary and sufficient in order that the four 
lines AX, BY, CZ, DT shall be generators of a hyperboloid: 


cos c’ 1 cos c’ i 


+ — —_-+— 
cos c’ 1 A B Cc D i 
—4+—4+—4+— 
oY 1 1 
——+—=——+-, 
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cos a 1 cos a 1 


ty ty v1 


Hence =h. 

The point P coincides with the center I of the sphere inscribed in T and the 
hyperboloid is the one formed by the lines joining the vertices to the points of 
tangency of the opposite faces with the inscribed sphere. These lines are con- 
current if 


cos 4a cos 3a’ = cos 3b cos 3b’ = cos 3¢ cos 3c’, 


that is, if the tetrahedron is isogonic.* 

As a consequence we have the theorem: In any tetrahedron T=ABCD let the 
perpendiculars at A to the faces CDA, DAB, ABC, meet the bisecting planes 
ICD,'IDB, IBC,---, in Ai, As, Then the perpendiculars from 
A, B, C, D, to the planes AiA2A3, BiB2B3, - + -, pass through the points of con- 
tact of the inscribed sphere with the opposite faces. 

All the properties apply to the triangle. Given a triangle ABC and a point 
P in the same plane. The perpendiculars at A to AB and AC meet BP and CP 
in A, and A; we have the analogous points B; and Be, C; and C2. The perpendicu- 
lars from A, B, C to AiA2, BiBo, C:\C2 meet BC, CA, AB in X, Y, Z. Let (AB, 
AP) =a, (AP, AC) =a’, (BC, BP) =8, (BP, BA) =’, (CA, CP) =v, (CP, CB) 
We have 

BX/XC=tan y/tan CY/YA =tan a/tan y’, AZ/ZB =tan B/tan a’, 
and 


Now, if A’, B’, C’ are the orthogonal projections of P on BC, CA, AB, we 
have also 


(2) 


A'C B’A C’B tana’ tan tan 7’ 


When point P describes the cubic of Darboux circumscribed about the tri- 
angle ABC, the lines AA’, BB’, CC’ are concurrent; the two members of equa- 
tion (2) have the value —1, and, by (1), the lines AX, BY, CZ meet in a point 
Q. Conversely, if lines AX, BY, CZ are concurrent, point P describes the cubic 
of Darboux. 

The locus of point Q is the cubic of Lucas, and the perpendiculars to BC, CA, 


AB from X, Y, Z meet in the isogonal conjugate of P with respect to triangle 
ABC. 


* See Nathan Altshiller-Court, Modern Pure Solid Geometry, New York, 1935, p. 290, art. 
879. 


+ We do not know if this method of generating the cubics of Darboux and Lucas has been 
noticed heretofore. (V.T.) 
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In particular, if in a triangle ABC, P coincides with the orthocenter, points 
X, Y, Z are the midpoints of BC, CA, AB, and Q coincides with the centroid. 

When P describes the circumscribed circle, AX, BY, CZ are parallel to BC, 
CA, AB. 

Finally, if P coincides with the center of one of the tritangent circles, X, Y, 
Z are the points of contact of this circle with BC, CA, AB, and Q is the Ger- 
gonne point of the triangle, or one of its associated points. 


INSTRUCTION AND RESEARCH IN APPLIED MATHEMATICS* 


The history of mathematics reveals a continuous interplay of ideas between 
the mathematician and the experimental or practical scientist. In the last half- 
century this interplay has decreased, to the detriment, it may be believed, of 
both parties. The war caused a lively cooperation, but unless this continues 
under peacetime conditions, the prospect for the future is serious and warrants 
earnest consideration. 

The widening gap between the mathematician and other scientists may be 
traced to several causes. The most obvious of these is the great activity in all 
branches of science, which has necessitated a high degree of specialization in 
each branch, and consequently less opportunity for roaming interests. However, 
side by side with the expansion of science, increased understanding of basic 
theory acts as a simplifying compensation. Such simplification is achieved 
largely through mathematical methods, of increasing depth and generality, 
and it would be indeed unfortunate if mathematicians should not be in a posi- 
tion to contribute towards this simplification. 

As our educational system is at present organized, there are relatively few 
opportunities for a student of mathematics, undergraduate or graduate, to be- 
come acquainted with the mathematical structures of the theories underlying 
other sciences. His needs in this connection are not identical with those of stu- 
dents pursuing other branches of science as their major interest; many practical 
details, of importance to them, are not essential for his purpose, and he can go 
further and faster by concentration on the mathematical aspects of the subject. 
It is suggested that the study of some branches of other sciences from the mathe- 
matical standpoint (and that is a meaning commonly attached to the words 
“applied mathematics”) should be regarded, wherever feasible, as an essential 
part of the training of mathematics majors, undergraduate or graduate. Apart 
from the general advantage to science which we might hope to see as aresult 
of this procedure, it seems only fair to offer a young mathematician the oppor- 
tunity of pursuing a career in applied mathematics if his natural inclinations 
are so directed. Unless he has the opportunity of viewing the field of applied 


* Prepared by the special committee of the American Mathematical Society on applied mathe- 
matics, and approved by the Council of the Society and the Board of Governors of the Association. 
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mathematics from the mathematical standpoint, it is unlikely that he will be 
attracted to it if his interests are basically those of a mathematician. 

Consequently it is suggested that departments of mathematics throughout 
the country should consider the feasibility of enlarging their offerings in the 
direction of applied mathematics, both on the undergraduate and graduate 
levels, insofar as the departments are qualified to offer such instruction. Many 
branches of applied mathematics have not only well established axiomatic 
mathematical structures, but are also fruitful fields of research, involving mathe- 
matical ideas and techniques of the highest order. Attention may in particular 
be directed to mathematical statistics, theoretical mechanics (including elas- 
ticity and fluid dynamics), statistical mechanics and thermodynamics, heat 
conduction, electromagnetic theory, relativity, quantum mechanics, genetics, 
and the theory of high polymers. 

In some institutions the problem is already solved in part by the activities 
of other departments, containing members well qualified in mathematics. 
Nevertheless, since each of the subjects listed above has not only its mathe- 
matical side, but also an even greater experimental or practical side, it appears 
likely that full justice can be done to the mathematical theories only by a 
specialization on that aspect. Needless to say, the task of maintaining and in- 
creasing an interplay between mathematics and the other sciences is one that 
can be dealt with only by full cooperaton between all the parties concerned. It 
is to be expected that scientists in other fields would welcome the desire of 
young mathematicians to increase their knowledge of the mathematical struc- 
tures of those fields, for the sake of a wider mathematical perspective and in 
the case of those qualified to pursue research, with a view to solving outstand- 
ing problems and simplifying basic theory by mathematical methods of in- 
creased generality and power. 

A similar gap in the background of our students lies in the field of computa- 
tion. Even if we carefully distinguish mathematical computation and com- 
putation engineering, there is much of importance concerning numerical 
computation and its organization for hand computer, IBM equipment, relay 
calculator or electronic computer, which will be of value to the student whether 
he is to become a topologist, a flutter analyst, or a mathematical logician. Of 
equal importance is the field of algebraic and analytic computation and its 
possibilities of mechanization. With the exception of a few institutions where 
the engineers have begun to work on computation engineering, this whole field 
is at the disposal of the departments of mathematics. 

Consequently, it is suggested that departments of mathematics throughout 
the country should consider the feasibility of enlarging their offerings in com- 
putation, algebraic, analytic, or numerical, at both the undergraduate and 
graduate levels. 
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THE WILLIAM LOWELL PUTNAM MATHEMATICAL 
COMPETITION 


G. W. MACKEY, Harvard University 


The following results of the seventh annual William Lowell Putnam Mathe- 
matical competition held May 24, 1947, have been determined in accordance 
with the rules of the Competition agreed to by the representatives of the Mathe- 
matical Association and the trustees of the William Lowell Putnam Intercol- 
legiate Memorial Fund. The contestants were known to the reader only by 
number. 

The first prize, four hundred dollars, is awarded to the Department of 
Mathematics of Harvard University, Cambridge, Massachusetts. The mem- 
bers of the team were C. W. Hewlett, Jr., J. J. Newman, W. F. Stinespring; to 
each of these, a prize of forty dollars is awarded. 

The second prize, three hundred dollars, is awarded to the Department of 
Mathematics of Yale University, New Haven, Connecticut. The members of 
the team were Murray Gell-Mann, Murray Gerstenhaber, and Henry Otto 
Pollak; to each of these a prize of thirty dollars is awarded. 

The third prize, two hundred dollars, is awarded to the Department of 
Mathematics of Columbia University, New York, New York. The members of 
the team were Alex Heller, Harold Lehrer, and Maxwell Rosenlicht; to each of 
these, a prize of twenty dollars is awarded. 

The fourth prize, one hundred dollars, is awarded to the Department of 
Mathematics of the University of Pennsylvania, Philadelphia, Pennsylvania. 
The members of the team were Joachim Ehrman, Donald F. Hunt, and William 
Turanski; to each of these, a prize of ten dollars is awarded. 

The five persons ranking highest in the examination, named in alphabetical 
order, were: C. W. Hewlett, Jr., Harvard University; Maxwell Rosenlicht, Co- 
lumbia University; W. F. Stinespring, Harvard University; William Turanski, 
University of Pennsylvania; Eoin L. Whitney, University of Alberta. Each of 
these will receive a prize of forty dollars. 

The five succeeding persons ranking highest in the examination, named in 
alphabetical order, were: Leonard Geller, Brooklyn College; Harry Gonshor, 
McGill University; Julian Kielson, Brooklyn College; John F. Nash, Jr., Car- 
negie Institute of Technology; Henry Otto Pollak, Yale University. Each of 
these will receive a prize of twenty dollars. 

The following teams, named in alphabetical order, won honorable mention: 
Brooklyn College, Brooklyn, New York, the members of the team being Melvin 
Hausner, Robert Margolies, and George Shapiro; University of California, Ber- 
keley, California, the members of the team being Marvin P. Epstein, George E. 
Gourrich, and Roger A. Stafford; Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania, the members of the team being George W. Hinman, John F. 
Nash, Jr., and Robert I. Van Nice; Northwestern University, Evanston, IlIli- 
nois, the members of the team being Donald B. MacMillan, James E. Murrin, 
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and Ernest Tilden Parker; Queen’s University, Kingston, Ontario, the mem- 
bers of the team being Robert W. Butcher, Thomas G. Donnelly, and Donald W. 
Dunn; University of Toronto, Toronto, Canada, the members of the team being 
W. T. Sharp, Miss M. J. Straus, and W. J. D. Lewis. 

Fifteen individuals were given honorable mention. The names are listed in 
alphabetical order. G. F. D. Duff, Toronto University; Joachim Ehrman, Uni- 
versity of Pennsylvania; Murray Gell-Mann, Yale University; Murray Gersten- 
haber, Yale University; George E. Gourrich, University of California; Melvin 
Hausner, Brooklyn College; Donald B. MacMillan, Northwestern University; 
J. J. Newman, Harvard University; John F. Riordan, Massachusetts Institute 
of Technology; William Riordan, Massachusetts Institute of Technology; George 
Shapiro, Brooklyn College; W. T. Sharp, University of Toronto; Roger A. 
Stafford, University of California; James E. Storer, Cornell University. 

The following is a list of all colleges and universities which entered teams in 
the Competition. The list, in alphabetical order, is: University of Alberta, Bos- 
ton College, Brooklyn College, University of California (Berkeley), University 
of California at Los Angeles, Carleton College, Carnegie Institute of Technol- 
ogy, College of St. Thomas, University of Colorado, Columbia University, 
Cornell University, Harvard University, Haverford College, Loyola College 
(Montreal), Loyola College (New Orleans), Massachusetts Institute of Technol- 
ogy, McGill University, Northwestern University, Oklahoma Agricultural and 
Mechanical College, University of Oklahoma, University of Pennsylvania, 
Queen’s College, Queen’s University, Rutgers University, Swarthmore College, 
Texas Technological College, University of Toronto, Ursinus College, United 
States Naval Academy, Wayne University, University of Washington, and 
Yale University. 

The following additional colleges and universities entered individual con- 
testants only: Case School of Applied Science, Purdue University, University 
of Saskatchewan, and Washington University. 

A total of 145 undergraduates representing 36 institutions took part in the 
competition. 

Participants in the competition were given the following lists of problems. 


Part I. THREE Hours 


(Answer the questions in any order and by any method. Show all your work in logical 
sequence, and indicate your answers clearly. No tables or other books may be used.) 


1. If {a,} is a sequence of numbers such that for »=1 
(2 Gn) An41 = 


prove that lim a,, as n— ©, exists and is equal to one. 
2. A real valued continuous function satisfies for all real x and y the func- 
tional equation 


+ y*) = f(x)f(y). 
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Prove that 


f(x) = 


3. Given this figure and any two points Q,, Q2 in the plane not lying on any 
of the segments si, se, - - - , 5s, Show that there does not exist a polygonal line 
P joining Qi, and Q2 such that: 


(1) P crosses each s;,i=1, 2, - - + , 6, exactly once; 
(2) P does not intersect itself; 
(3) P does not pass through any vertex V, Vi, Ve, Vs. 


V2 


V3 Se Vi 


4. A coast artillery gun can fire at any angle of elevation between 0° and 
90° in a fixed vertical plane. If air resistance is neglected and the muzzle velocity 
is constant (=v), determine the set H of points in the plane and above the 
horizontal which can be hit. 

5. a1, by, ¢; are positive numbers whose sum is 1, and for n=1, 2, - + + we 
define Show that aa, Dn, Cn 
approach limits as n— © and find these limits, 

6. A three by three matrix has determinant zero, and has the further prop- 
erty that the cofactor of any element is equal to the square of that element. 
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(The cofactor of a;; is (—1)* multiplied by the determinant obtained by strik- 
ing out the ith row and jth column.) Show that every element in the matrix is 


zero. 
Part II. THREE Hours 


(Answer the questions in any order and by any method. Show all your work in logical 
sequence, and indicate your answers clearly. No tables or other books may be used.) 


7. Let f(x) be a function such that f(1) =1 and for x21 
1 
f(x) = 2+ f(x) 
Prove that 
lim f(2) 


exists and is less than 1+7/4. 
8. Let f(x) be a differentiable function defined in the closed interval (0, 1) 


and such that 


| f(x) | O<x<1. 
Prove that 
1 12 k M 
0 N kel n n 


9. Let x, y be Cartesian codrdinates in the plane. I denotes the line segment 
1<x<3, y=1. For every point P on J, let P* denote that point that lies on the 
segment joining the origin to P and such that the distance PP* is equal to 
1/100. As P describes J, the corresponding point P* describes a certain curve 
C*. Let 1(I), 1(C*) be the lengths of I and C* respectively. Which one of /(J), 
1(C*) is greater? Prove your answer. 

10. Given P(z) =z*+-az+8, a quadratic polynomial of the complex variable 
z with complex coefficients a, 6. Suppose that | P(z)| =1 for every z such that 
|z| =1. Prove that a=b=0. 

11. a, b, c, d are distinct integers such that 


(x — a)(x — b)(x — c)(x —d) —4=0 


has an integral root r. Show that 4r7=a+b+c+d. 

12. Cis a fixed point on OZ and U, V are variable points on OX, OY respec- 
tively, where OX, OY, OZ are mutually orthogonal lines. Find the locus ofa 
point P such that PU, PV, PC are mutually orthogonal. 
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MATHEMATICAL NOTES 


EpiTep By E. F. BECKENBACH, University of California 


Material for this department should be sent directly to E. F. Beckenbach, University of 
California, Los Angeles 24, California. 


ON THE NUMBER OF PATHS IN A FINITE PARTIALLY ORDERED SET 


E. W. CHITTENDEN, University of Iowa 


Let P be a finite partially ordered system* with initial elements ay, d2, + + + , Gm 
and terminal elements y, bo, - - + , bn. Let N be a matrix (x;;), (@=1,2,+++,m; 
j=1, 2,--+,m), where x;; is the number of different paths from a; to b;, where 


a path is a maximal simply ordered subset of P which contains a; and b;. We 
show that the matrix N can be computed in terms of a series of incidence matri- 
cest which we proceed to define. The initial elements of P, ai, +--+, @m, are 
elements of rank zero. Anelement p of P is of rank 1 if it is not of rank zero and 
all its predecessors are of rank zero. In general an element p of P is of rank r+1 
if it has a predecessor of rank r and no predecessor of higher rank. If p’ is of rank 
r' and p”’ of rank r’’>r’+1, p’’ isa successor of p’, and there is no element p of 
P of rank r’<r<r’’ between p’ and p’’, we agree that p’’ may be counted 
among the ranks r’<r<r’’ and that in passing from rank 7 to rank r+1, p’’ 
is regarded as its own successor. Thus as there is actually only one path from p’ 
to p’’ in P, the convention just adopted permits us to admit the existence of a 
path from p’ in rank 7’ to p’’ in rank r+1 without actually increasing the total 
number of paths in the system. 

The incidence matrix A,=(e;;) is defined by the rule: If p;, is an element 
of P of rank r and is followed by p;,-4:, then e;;=1. Otherwise e;;=0. Let k be 


the maximal rank of any element of P. Then the 0; are of rank k and N is the 
matrix product 


r=k 
N 


r=0 


The proof is easily made by induction. Evidently Ao is the number of paths 
from elements of rank zero to their successors, and A, is the number of paths 
from an element of rank r to an element of rank r+1. Hence if N, is the matrix 
giving the number of paths from rank 0 to rank r, we see that the ith row of N, 
will give the number of paths from a; to each element p of rank r. The elements 
(e;;) of A, which are different from zero in a column give the incidences from an 
element of rank r to a particular element p; of rank r+1. Thus the element 
xi; of the product matrix N,-A, becomes the number of paths from a; to pj. 

In the case of a checkerboard, the number of paths from one king row to the 


* Garret Birkhoff, Lattice Theory, American Mathematical Society, 1940. 
ft O. Veblen, Analysis Situs, Colloquium Lectures, American Mathematical Society, 1916. 
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other is given by the matrix 


14 14 6 £1 

28 34:21 6 

20 35 14 14 

7 22 28 14 


Thus the number of paths, without jumps, from the second element of the start- 
ing king row to the third element of the ending king row is 21. 


NOTE ON CONJUGATE HARMONIC FUNCTIONS 


Epwarp Kasner, Columbia University 
Joun De Cicco, Illinois Institute of Technology 
In many textbooks on functions of a complex variable, the problem of finding 
a harmonic function (x, y) conjugate to a given harmonic function ¢(x, y) is 
solved in the following manner. Since the functions @ and y are conjugate- 
harmonic, they obey the Cauchy-Riemann equations 


(1) oz: = Vy dy = — vx, 

and also each of the functions satisfies the Laplace equation, that is 

(2) t+ dy = 0, Ver t = 0. 

Hence the function y is found by integrating the equation 

(3) dy = — o,dx + 

which is exact since the given function ¢ is harmonic. Moreover all such func- 
tions y differ from one another merely by a real constant. 

We shall show how to find all such functions y without resorting to integra- 
tion. The advantage of our method is that it enables one to prove certain 
theorems concerning the conjugates of rational, algebraic, or entire harmonic 
functions. 

In the first place, ¢(x, y) is harmonic and hence is an analytic function of 


x and y. Therefore by the obvious relations, u=x+iy and v=x—vty, we can sub- 
stitute x =(u+v)/2 and y=(u—v)/(27) into this function. The result is 


] = + no, 


2° 
where y(v) is an analytic function of v with coefficients the conjugates of the cor- 
responding ones of the analytic function \(x). 

Since we are searching for the second component (x, y) of the conformal 
transformation: X =¢(x, y), Y=y(x, y), we know that ¢+ip=2f(u) and 
@—iy =2g(v), where g(v) is an analytic function in v with coefficients the con- 
jugates of the corresponding ones of the analytic function f(u). Thus we have 


(5) d=etf, w=ilg—/f). 


(4) = 
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Comparing (4) and (5), we find 


(6) ic/2, 
where ¢ is a real constant. Therefore the solution for yp is 
(7) v= i(u—dA) +e 


Thus from (7), the function y may be deduced immediately. 
As an example, consider the problem of finding the conjugate of the har- 
monic function ¢=x/(x?+~?). Performing the obvious substitutions, we find 


= x/[(x + iy)(x — iy)] = (uw + = 1/(20) + 1/(2u), 
so that \=1/(2u) and »=1/(2v). Substituting into (7), we have 


= (i/2)(1/0 — 1/u) + ¢ = i(u — 0)/(2uv) +o = — y/(x? + +e, 


The related analytic function is @6+i)=1/u-+ic. 

As another example, let us find the conjugate of the harmonic function 
@=arce tan y/x. First recall that arc tan m=(1/2i) log [(1++-im)/(1—im) ]. 
Hence we have 


@ = arc tan y/x = (1/2i) log [(x + iy)/(x — iy)] = (1/22) (log u — log 2). 
Thus \=(1/27) log wu and p=(—1/22) log v. By (7), we find 
y = (— 1/2)(log u + log v) + ¢ = (— 1/2) log uv +c = — log (x? + y*)/? +c. 


The related analytic function is ¢+i~ = —i log u+ic. 

It is evident that the conjugate of a harmonic polynomial is also a harmonic 
polynomial. However it is not so obvious that the conjugate of a rational har- 
monic function is also rational. 

The conjugate of a rational harmonic function is rational. Thus if one com- 
ponent of an analytic function is a rational function of (x, y), then the analytic 
function is rational in u=x+iy. 

For if @(x, y) is a rational function of (x, y), obeying the Laplace equation, 
it follows from (4) that \(u) and y(v) are both rational functions. Hence by (7), 
it follows that y is rational also. Finally the function ¢+17y is a rational func- 
tion of u=x+iy. 

The conjugate of an algebraic harmonic function is algebraic. Therefore an 
analytic function of u=x-+1y is algebraic if only one component is given to be an 
algebraic function of (x, y). 

This result is also a consequence of (4) and (7). 

Similarly the conjugate of an entire harmonic function is entire. Hence an 
analytic function of u=x-+iy is entire if only one component is given to be an 
entire function of (x, y). 

For if @ (x, y) is an analytic function of (x, y) convergent for all finite values, 
then A(u) and y(v) converge for all finite values of u and v. Hence by (7), 
¥(x, y) is an analytic function of (x, y) convergent for all finite values. Therefore 
is an entire function of u=x++y. 
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CLASSROOM NOTES 


EpiTep By C. B. ALLENDOERFER, Haverford College 
All material for this department should be sent to C. B. Allendoerfer, Haverford College, 
Haverford, Pennsylvania. Contributions are invited on topics of immediate interest to teachers 
of undergraduate mathematics. 


END-POINT MAXIMA AND MINIMA 
C. O. OaKLEy, Haverford College 


It sometimes happens that a very elementary maximum-minimum problem 
will lead to seemingly impossible or even foolish results. The following examples 
are of this category and belong to the type commonly known as end-point 
maxima. They are of such simple nature that they may be readily assigned to a 
beginning class in the calculus.* Although not original, these two particular prob- 
lems do not seem to be widely known and should prove to be both interesting 
and instructive to teachers and students alike. 


PROBLEM I. Find the position of the point P(x, y), on an ellipse, such that the 
distance to a focus F is a maximum (minimum). 


Let the ellipse have the equation b?x?+a’y?=ab?. Then the distance PF is 
a continuous function of x in the interval (—a, a) and hence it has both a 
maximum and a minimum value. (Any function which is continuous in a closed 
finite interval possesses both a maximum value and a minimum value in that 
interval). Now the fortunate student, bright or dull, might begin to solve this 
problem as follows: 

SOLUTION 1. Let the ellipse be represented parametrically by x=a cos 8, 
y=b sin 6 and let F be the right hand focus (c, 0). Then 


L = PF? = (a cos — c)? + B? sin? 8, 


and 

dL 

—* 2(a cos 6 — c)(— asin 6) + 28? sin @ cos @ = 0, 
from which it follows that sin @=0, or cos @=a/c. The solutions of sin @=0 yield 
the answers while cos @ cannot equal a/c since a/c>1. 

But even an able student might attack the problem differently—and come to 
grief. 

SoLuTION 2. Now PF=a-—ex. It’s derivative is —e#0. And yet a—ex is 
obviously least when x =a and greatest when x = —a, if the condition -a Sx Sa 
is kept in mind. Even though the square of the distance is used, L=PF* 
=(a—ex)*, the unwary student is still in trouble for dL/dx = —2e(a—ex) =0 
yields x =a/e, a point outside of the ellipse (on the directrix, indeed). The func- 


*See J. L. Walsh, A Rigorous Treatment of the First Maximum Problem in the Calculus, 
this MonTuLy, vol. 54, (1947), pp. 35-36. 
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tion L does have a (relative) minimum at x=a/e (what about the maximum?) 
but there is no (real) y-codrdinate for the corresponding point P on the ellipse. 
Whether PF or PF? is used, the end-points of the interval determine the 
(absolute) maximum and the (absolute) minimum values wanted. 

Even the method that most beginning students are likely to employ will 
result in the same difficulty: 
SOLUTION 3. 


L = PF = (x — c)? + ¥?, 
2 

= (2-08 += (a 24); 
a? 


or again, 


x=a?/c = a/e. 


Pros_eM II. A man is in a boat at P one mile from the nearest point A on 
shore. He wishes to go to B which is farther down the shore M miles from A. If 
he can row r miles an hour and walk w miles an hour, toward what point C should 
he row in order to reach B in least time? 


Usually this problem is given with simple numerical values that present no 
difficulty. The shore line is assumed to be straight; let x=AC. Then the total 
time ¢ that it takes to go from A to B is readily determined as 


1 1 
t= —(i+ x’)'/?4+ —(M — 2), 
r w 


dt x 1 
dx r w 


— 


The first astonishing thing about this “answer” is that it is independent of 
the distance M. If r<w and if r<M+/w*—?’, routine testing of the derivative 
will show that a minimum time is attained. But if r<w and r>M+/w?—?, then 
the derivative vanishes for a value of x that lies outside of the interval AB and 
the problem has an end-point minimum, that is, the man should row directly 
toward B. Finally if r2w, the man should again row directly to B. 


3 dL b? 
dx a? 
whence 
0, 
and 
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Plotting the graphs of the functions involved in these problems is a great aid 
in clarifying the points at issue. 

DiscussIon. Problems involving end-point maxima and minima usually 
arise only when the independent variable is restricted to a certain interval in 
order to insure reality of the function or for certain physical reasons such as 
the avoidance of negative values of distance or time. In Problem I we found 
that the answers were either relative maxima and minima or end-point maxima 
and minima according to the particular choice of the independent variable. This 
is indeed a general situation: extreme values of the one type can be converted 
into those of the other type by an appropriate choice of a parameter or, what 
amounts to the same thing, by an appropriate transformation on the inde- 
pendent variable. Indeed let the function to be maximized be y=f(x) and let 
the range of x be (a, b). Then put x=a cos? 0+5 sin? 9. This transformation 
restricts x to the interval (a, b) but permits the new independent variable @ to 
vary from — © to +o, It is readily seen that now dy/d0@ = (dy/dx)2(b—a) sin 
cos 9 which vanishes at @=0 and 7/2, namely, at the points x =a and x=), as 
well as at any relative maxima and minima inside the interval (a, b). Applying 
this transformation to the above exercises will convert them into the usual 
relative maxima and minima problems. For example in Problem II if we set 
x=M sin? 0, we find dt/d@ vanishes when sin @=0, cos 6=0 and when 
sin‘ 9 =r?/M?(w?—r?). An analysis of these values gives the results stated above. 

If the interval is an infinite one, say (a, ©), the transformation x =a+tan?@ 
will change an end-point maximum (minimum) at x=a into a relative one. 

Editorial Note. Another example of an end-point maximum has been sub- 
mitted by V. L. Klee, Jr. of the University of Virginia. It is as follows: 

PROBLEM. We are given a straight fence 100 feet long, and wish by adding 200 
feet more to form a rectangular enclosure whose boundary contains the original 
fence. How shall this be done so as to enclose the greatest possible area? 

Let x denote the length of new fence which is aligned with the original 100 
feet of fence. Then the dimensions of the enclosure are 100+ and 50—x, and 
its area A(x) is —x?—50x+5000. The student may reason as follows: D,A 
= —2(x+25); hence we would choose x = —25 in order to enclose the maximum 
area. He will then be perplexed to note that this solution obviously does not 
satisfy the conditions of the problem. The correct solution, of course, is x =0 
which is an end-point maximum. 


LETTERS TO THE EDITOR 
The Trial Integral Method 


The essential idea of the method explained in my note, The Trial Integral 
Method, this MonTHLY, vol. 34 (1947) pp. 159-160, is contained in a paper by 
T. H. Hildebrandt, Marginal Notes, this MONTHLY, vol. 36 (1929), pp. 216-221. 

M. F. SMILEY 
Northwestern University 
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The Equation of an Ellipse 


In this Montuty for April (vol. 54 (1947) pp. 219-220) there appears a note 
by Frank Hawthorne in which the Cartesian equation of the ellipse with respect 
to its axes is derived, without the use of radicals, from the constancy of the sum 
of the focal radii. May I point out that the derivation of the equation of the 
ellipse in this manner appeared in 1707 in the posthumous Traité analytique des 
sections coniques (pp. 22-25) of the Marquis de L’Hospital. A concise account 
of this and other portions of the work of L’Hospital is 3iven in J. L. Coolidge, 
A history of the conic sections and quadric surfaces, Oxford, 1945, p. 77. 


C. B. Boyer 
Brooklyn College 


The Remainder Theorem 


The basic ideas of the note by R. W. Wagner, An A pplication of the Remainder 
Theorem, this MONTHLY, vol. 54 (1947) p. 106, were contained in a lecture by L. 
C. Karpinski, “Mathematical Short Cuts for Engineers and Architects,” pub- 
lished in a lithographed form by the Multi-Color Company of Detroit, Michigan 
in February 1930. 

Karpinski extended this idea to make an interesting and pedagogically de- 
sirable connection with a common arithmetical concept ,the “check by nines.” 
The proof of this arithmetical trick, taught usually by rule, if at all, follows, of 
course, from the fact that if, as in Wagner’s note, p(x) represents an integer 
conceived of as a polynomial in 10, then (1) is both the remainder if the number 
is divided by x—1, and the sum of the digits. In this case, x—1 is 9. The “check 
by casting out nines” then follows from the fact that if p, g, r are integers, then 
(p+qt+r)/9=p/9+q/9+1r/9 and hence the remainder of the sum should equal 
the sum of the remainders (with nines cast out). The extension of the check to 
the other fundamental operations of arithmetic is obvious. 

Karpinski pointed out in his lecture that one may also derive checks by 11, 
99, 101, and in a recent conversation he added that if the number is written 
as a polynomial in 1000 (e.g. 1,484,365 = p(x) =x?+484x+4365), then p(—1) is 
the remainder after division by 1001 (in the example p(—1)=—118). Since 
1001 =7-11-13, the divisibility of the remainder by any of these factors is 
necessary and sufficient to show the divisibility of the original integer. 


P. S. Jones 
University of Michigan 


Linear Differential Equations 


The purpose of these remarks is to express the opinion that the note On the 
Solution of Linear Differential Equations in this MONTHLY, vol. 54 (1947) p. 160 
is in violation of (a) sound mathematics (b) sound pedagogy (c) sound common 
sense 

I will briefly restate the author’s thesis, using his first example, the equation 


yy" 24’ + y = e%, 
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He sets down the general solution, and observes that as it stands it is not valid 
when a=1, containing as it does a denominator (a—1)*. He then rigs a special 
form which when a approaches unity does not become infinite but is inde- 
terminate, and naively assumes that the limit of this expression is the solution 
when a=1. This is to assume the validity of a theorem which certainly has not 
been proved at this stage of the game, namely, that if a solution of 
f(x, y, y’, a) =0 is bounded, it is a continuous function of a. The result, in the 
case before us, happens to be correct, but the process is unjustifiable. 

What does the conscientious mathematician do when a general result be- 
comes meaningless in a special case? He does not try to patch it up by ingenuity 
but goes back to the original problem and works it through. In this instance 
whether or not we like to use the symbolic operator D, we have a well-defined 
process for solving linear equations with constant coefficients, namely the repeti- 
tion of the relationship: 


If y’ — ay = f(x), y= f + |. 


Direct application of this rule yields easily the desired solution of the equation 
y’’ —2ay’’ +a*y =e* in a neat and simple form. 

Pedagogically, the operations which I suggest are straightforward and direct. 
The author remarks that the usual methods seem to involve pulling something 
out of thin air, but it seems to me that on the contrary his method involves 
pulling a rabbit out of a previously prepared hat. What student could see the 
sense in setting up his peculiar form of the general solution, or would understand 
it well enough to be able to handle another similar situation? The process which 
I suggest, on the other hand, goes back to sound and well-established principles. 
This argument leads naturally and logically to my third point, sound mathe- 
matics and pedagogy are bound up with common sense. We attack any given 
problem on its own merits, not as an offshoot of some other problem. 


R. A. JOHNSON 
Brooklyn College 
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ELEMENTARY PROBLEMS AND SOLUTIONS 


EpiTep By Howarp Eves, Oregon State College 


Send all communications concerning Elementary Problems and Solutions to Howard 
Eves, Mathematics Department, Oregon State College, Corvallis, Oregon. This department 
welcomes problems believed to be new, and demanding no tools beyond those ordinarily fur- 
nished in the first two years of college mathematics. To facilitate their consideration, solutions 
should be submitted on separate, signed sheets, within three months after publication of 
problems. 

PROBLEMS FOR SOLUTION 


E 781. Proposed by P. D. Thomas, Navy Department, Washington, D.C. 

A heavy ball is gently dropped into a vase full of water. A section through 
the vertical axis of the vase is a semiellipse, the height being the semimajor 
axis, the diameter being then the minor axis. The size of the ball is such as to 
cause the maximum displacement. (1) Find the radius of the ball. (2) Show that 


the plane of the circle of tangency bisects the height of the submerged segment of 
the ball. (See E 687 [1946, 334].) 


E 782. Proposed by Joseph Rosenbaum, The Milford School, Connecticut 


Show that the product of two numbers each of the form x?+kxy+y?, where 
k, x, y are integers, k fixed, is also of that form. If x and y are relatively prime, are 
all factors of x?+kxy+y? also of the same form? 


E 783. Proposed by C. D. Olds, San Jose State College 


Given a parallelogram and its diagonals. Let each side of the parallelogram 
be divided into m equal parts and let lines be drawn through the points of di- 
vision, parallel to the sides and to the diagonals of the parallelogram. Find the 
total number of triangles in the resulting figure. 


E 784. Proposed by R. E. Gaines, University of Richmond 

Show that the locus of the intersection of two successive perpendicular 
tangents to a logarithmic spiral is another logarithmic spiral. 

E 785. Proposed by R. J. Walker, Cornell University 


Each of »—1 tanks, Ty, - + +, Tn-1, holds V gallons of water, and an nth tank, 
T,, holds V gallons of a salt solution containing M pounds of salt. Liquid is 
circulated at the rate of g gallons per minute from T, to Tn-1, Tn-1 tO Tn-2, ++ *, 
Tz, to T;, T; to T,. How much salt is in T, after ¢ minutes? 


SOLUTIONS 
E 743, Solution II [1947, 285]. Correction. In the third line of the solution 
the word “positive” should be replaced by “negative.” 
A Cryptarithm 
E 751 [1947, 38]. Proposed by Alan Wayne, Flushing, N. Y. 
Find the digits represented by the letters in the following addition, if no two 
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different letters represent the same digit: 
FORTY 
TEN 
TEN 


SIXTY 


Solution by A. Chulick, Hofstra College. Since the sum of column (1), (right 
to left), is Y and no digit can be carried, N is zero. Likewise, since the sum of 
column (2) is 7, or T plus 10, E is either zero or five. But N is zero; therefore E 
is five. A digit must be carried from column (4). That digit is 1, for the sum of 
column (4) can be at most 11. It follows that J is either zero or one. But N is zero; 
therefore J is one, and O is nine. The sum of column (3) must be at least 21, for 
X cannot be either zero or one, and at most 1 can be carried from column (2); 
hence T must be either seven or eight. If T is seven R must be eight and X three. 
But F and S are consecutive digits. The only possibilities are two and three, 
and three and four; therefore X cannot be three. T must be eight, R seven, and 
X four. The only consecutive digits available are two and three; hence F is two 
and S is three. Y is six. The solution is unique and we have 29786+850+850 
= 31486. 

Also solved by Richard Andree, Maurice Anthony, LeRoy Babcock, Murray 
Barbour, T. E. Berry, Daniel Block, R. A. Bradley, W. G. Brady, Paul Brock, 
W. E. Buker, M. I. Chernofsky, P. L. Chessin, W. H. Coulter, R. E. Crane, 
J. A. Cromelin, J. H. Cross, J. E. Darraugh, E. de la Garza, Monte Dernham, 
William Douglas, R. L. Duncan, Virginia Felder, J. H. Ferguson, Barkley Fritz, 
I. M. Gardoff, M. A. Geisler, E. St. John Gough, E. L. Harp, Jr., R. H. Hoskins, 
J. M. Kingston, Mrs. V. L. Klee, Jr., William Kruskal, Elmer Latshaw, H. N. 
Leifer, G. E. McAllister, Burnett Meyer, Leo Moser, Roger Osborn, S. T. 
Parker, C. R. Perisho, C. F. Pinzka, P. A. Pisa, A. O. Qualley, P. W. A. Raine, 
A. P. Rhodes, G. G. Roberts, E. D. Schell, David Sohn, G. W. Walker, Davis 
Wellinger, Hazel Schoonmaker Wilson, M. S. Constable and the proposer. 


Editorial Note. Problems of this nature are appropriately called cryptarithms, 
and constitute a special class of the more general problems known as arith- 
metical restorations. See Ball-Coxeter, Mathematical Recreations and Essays, pp. 
20-26 and Kraitchik, Mathematical Recreations, pp. 79-80. Although the above 
problem is easy, cryptarithms can be very difficult. Devotees, however, prefer 
charming cryptarithms, and to be charming a cryptarithm should (1) make sense 
in the given letters as well as the solved digits, (2) involve all the digits, (3) 
have a unique solution, and (4) be such that it can be broken by logic, without 
recourse to trial and error. Judged by these standards the above cryptarithm is 
very charming indeed. 

The proposer has composed a number of these cryptarithms, some of which 
have appeared in The Cryptogram, publication of the American Cryptogram 
Association. Thus he has offered SEVEN+SEVEN+SIX =TWENTY (unique 
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solution), SEVEN+THREE+TWO=TWELVE (unfortunately possessing 
two solutions, since N and O may be interchanged), TWENTY+FIFTY 
+NINE+ONE=EIGHTY (unique solution but requiring an undesirable 
amount of trial and error), SCAN+THESE=DIGITS (unique solution), and 
RODE+ MAIL+ADDED+EMBLEM =OFFERED (having, as a literary ac- 
cessory, the key word FORMIDABLE, whose successive letters correspond to 
the ordered set of digits 0,1, ---, 9). 
Dernham and McAllister both suggested that if anyone wants more, let him 
SEND+MORE= MONEY. 
Of the proposed problem W. H. Coulter wrote: 
“Two sums that oddly cleave as one 
Though far apart like Earth and Sun.” 


Two Six-Point Circles Associated with a Right Triangle 
E 752 [1947, 38]. Proposed by Victor Thébault, Tennie, Sarthe, France 


Show that in a right triangle the twelve points of contact of the inscribed 
and escribed circles form two groups of six points situated on two circles which 
cut each other orthogonally at the points of intersection of the circumcircle 
with the line joining the midpoints of the legs of the triangle. 


Solution by J. H. Butchart, Arizona State College. Let ABC be the right 
triangle, C the vertex of the right angle, and let X, Y, Z be the points of 
contact of the incircle with the sides a, b, c opposite the vertices A, B, C respec- 
tively. Let Xi, X2, Xs be the points of contact of the escribed circles opposite 
A, B, C respectively with side a, and similarly define Y; and Z;. Let P, Q be the 
points where the line bisecting a, b at A’, B’ meets the circumcircle. Then, 
using the relations like AZ =s—a, BX;=s—a, where s is the semi-perimeter of 
triangle ABC, we can show that the sets of points P, Y, Q, Y3; P, X Q, Xs; 
X, Xs, Ys, Y are concyclic. For instance, (PA’)(A’Q) =(CA’)(A’B) =a?/4 
=(XA’)(A’Xs3). Since XX3, YY3 meet at C, which is not on PQ, it is impossible 
for circles (PYQY3) and (PXQX;3) to be distinct. Similarly Z;, Z2 lie on circle 
(XX3Ys3Y) and the center is the intersection of the perpendicular bisectors of 
XY and 2Z,Zz, i.e., an end, K’, of the diameter of the circumcircle perpendicular 
to AB. The points Z, Z3, X1, Q, Yi, X2, P, Y2 are proved concyclic in a similar 
manner, and K, the other end of the diameter through K’, is proved the center 
of this circle by showing that KQ and KZ are equal. Since the radii KQ, K’Q 
of the two circles are perpendicular, the circles are orthogonal. 

Also solved by the proposer. 

Editorial Note. It is believed that the two six-point circles of this problem are 
new in the geometry of the right triangle. 


The Elliptical Race Track 
E 753 [1947, 38]. Proposed by L. M. Kelly, University of Missouri 
How can one convince a class in elementary analytics that if the inside of a 
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race track is a non-circular ellipse, and the track is of constant width, then the 
outside is not an ellipse? 


I. Solution by Murray Barbour, Michigan State College. Let the equation of 
the inside of the track be 


x?/a* + y?/b? = 1, ax b. 

If the track is of uniform width and the outside ¢s an ellipse, its equation would 
have to be 

+ #)? + y*/(b+ 8? = 1, t #0. 


At any point of these two ellipses where the slope of the tangent is m, the width 
of the track will be 


{(a + + (6 + — {a%m? + 1 
(m? + 1)1/? 
After expanding and simplifying, this equation becomes 
t(am? + b) = t{ (m? + 1)(a%m? + 5%) }1/2, 

Since 0, both sides may be divided by ¢. Further simplification then gives 

2abm? = m?(a? + 5). 
Since this must be an identity for all values of m, we find 
a? — 2ab + b? = 0, 


(1) 


or a=b, contrary to our original hypothesis. Thus the outside of the track is 
not an ellipse. 


II. Solution by Leo Moser, University of Toronto. Assuming the inside and 
outside curves of the track to be the ellipses of Solution I, we note that the point 
A(a/s/2, b/4/2) lies on the inner ellipse and the point B((a+#)/+/2, (b+#)/+/2)) 
lies on the outer ellipse. Now the distance between A and B is ¢t. However, there 
is a shorter distance between the curves, for the slope of the line joining A to B 
is 1, while the slope of the tangent to the first ellipse, at the point A, is —b/a 
#—1. Thus we have a contradiction which proves the theorem. 

Also solved by R. A. Bradley, M. A. Geisler, E.G. Goman, and L. S. Shively. 

Shively obtained relation (1) above and suggested that one way to convince 
a class that the left member does not always equal ¢ is to substitute the special 
values a=4, b=3, t=1, m=1. Goman remarked that a class may be convinced 
that the inside and outside of an oval track of constant width ¢ cannot always 
both be ellipses by taking the outside of the track as an ellipse, and then con- 
sidering the inside curve for ¢ greater than the semi-minor axis of the outside 
curve. The inside curve clearly intersects itself, and thus cannot be an ellipse. 


Editorial Note. A curve, which is at a constant distance along the normal 
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from another curve C, is called a parallel to C. A curve and any of its parallels 
have the same normals and the same evolute. The parallel curves of a non- 
circular ellipse are curves of the eighth degree. See, e.g., Ex. 3, art. 372, of 
Salmon’s Conic Sections. R. C. Yates has devised a linkage for describing curves 
parallel to an ellipse. See this MonTHLY [1938, 607]. 


Skew Ordered Sequences 
E 754 [1947, 39 and 1947, 163.] Proposed by S. T. Thompson, Tacoma, Wash. 


A finite sequence of positive integers will be said to be skew ordered if either 
each integer in an even position of the sequence is greater than or each such 
integer is less than its immediate neighbors. If the eight integers 1, - - - , 8 are 
placed in random order in a sequence, what is the probability that the sequence 
will be skew ordered? 


I. Solution by J. B. Kelly, Hampton, Va. Let Q, be the number of skew 
ordered arrangements of m distinct integers. There is no loss in generality in 
supposing that the integers are 1, 2, - - - , m. Suppose that the integer occurs 
in the kth position. If we have a skew ordered arrangement, the sequences on 
either side of m must be skew ordered. The sequence to the left of m must have its 
last element less than its next to last element. Once the integers in this sequence 
are chosen, its elements may be arranged in 4Q,;_,; different ways so as to satisfy 
this condition. The sequence to the right of m must have its first element less 
than its second element. Once the integers in this sequence are chosen, its 
elements may be arranged in 3Q,_; different ways so as to satisfy this condition. 


There are iat ways of choosing the integers in the sequence to the left of 


n and once these integers are chosen, the integers in the sequence to the right 
of m are determined. Thus the number of skew ordered permutations of m dis- 
tinct integers for which the greatest integer (here 7) occurs in the kth position is 


(7-1) Qr-1Qn-x. It follows that 


In applying this formula, it is necessary to make the convention that Qo =Q:=2. 
Let P,, be the probability that a given permutation of m distinct integers will be 
skew ordered. Evidently P,=Q,/n!, and relation (1) becomes 


(2) P,= n—ky 


N kal 


where again we make the convention that Py) =P, =2. Calculating Ps by succes- 
sively calculating P2, Ps, - - -, Pz by means of (2) we obtain Pg=277/4032. 


II. Solution by Frederick Mosteller, Harvard University. The probability 
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that a finite sequence of m unequal integers will be skew ordered when all 
permutations are equally likely can be translated directly into statistical ter- 
minology. It is the probability that the sequence will have no run up or run down 
of length greater than or equal to 2. This problem was treated by P. S. Olmstead, 
Distribution of sample arrangements for runs up and down. Annals of Mathe- 
matical Statistics, vol. 17 (1946), pp. 24-33. 

For n=3 to 14 we derive from Olmstead’s Table 2 


n Probability of skew ordering 
3 0 .66666667 
4 0 .41666667 
5 0. 26666667 
6 0.16944444 
7 0.10793651 
8 0 .06870040 
9 0 .04373898 
10 0.02784447 
11 0.01772647 
12 0.01128499 
13 0.00718426 
14 0.00457364 


Since for m even the number of arrangements with runs of length no more 
than one is just twice Euler’s number E,, and for m odd the number of arrange- 
ments is twice the tangent number, we might use the approximation 


2\nt1 
“(=)”. 
(See, e.g., Milne-Thompson, Calculus of Finite Differences, p. 147.) 


For n=4 this approximation agrees to 0.002, for n=6 to 0.0001, for n=8 
to 0.000003, and for »=9 to 0.000001. The exact answer for n =8 is 277/4032. 
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ADVANCED PROBLEMS AND SOLUTIONS 


EpiTep By E. P. STarKE, Rutgers University 


Send all communications concerning Advanced Problems and Solutions to E. P. Starke, 
Rutgers University, New Brunswick, New Jersey. All manuscripts should be typewritten, with 
double spacing and with margins at least one inch wide. Problems containing results believed 
to be new or extensions of old results are especially sought. Proposers of problems should also 
enclose any solutions or information that will assist the editor. In general, problems in well 
known text books or results found in readily accessible sources should not be proposed for this 
department. 


PROBLEMS FOR SOLUTION 
4259. Proposed by Richard Bellman, Princeton University 


If 
= (1+ 2”), <1, 
show that, except perhaps for order, 
nk = 2°, 


4260. Proposed by Victor Thébault, Tennie, Sarthe, France 


In a triangle ABC inscribe two triangles A;B,C; and A2B.C, whose sides are 
parallel to the medians. Show that (1) the triangles ABC, A1B,C,, A2B2C2, have 
the same centroid and the same Brocard angle; (2) the triangles A:BiCi, A2BeC, 
are inscribed in an ellipse concentric and homothetic to the ‘inscribed Steiner 
ellipse, the ratio of homothety being 1/+/3. 


4261. Proposed by F. J. Dyson, Trinity College, Cambridge, England 


The number of partitions of an integer into a sum of positive integral parts 
is denoted by p(n). The result of subtracting the number of parts in a partition 
from the largest part is a positive or negative integer called the rank of the 
partition. Ramanujan proved that p(5”+4) is always divisible by 5, and 
p(7n+5) by 7. Show that the number of partitions of 5%+4 whose ranks are 
congruent modulo 5 to a given residue is the same whichever of the five residues 
is chosen, and the number of partitions of 7%+5 whose ranks are congruent 


modulo 7 to a given residue is the same whichever of the seven residues is 
chosen. 


4262. Proposed by L. A. Santalé, Rosario, Argentina 


Let C be a rectifiable plane curve of length L, contained within a given circle 
of radius R. Prove that there is a circle of radius p= R which cuts C in m points, 
where 


(1) n = L/rR. 
In particular there is a line which cuts C in m points, where n satisfies (1). If 
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p<R, the inequality (1) must be replaced by 


4Ip 


"= 


See L. A. Santalé, A theorem and an inequality referring to rectifiable curves, 
American Journal of Mathematics, 1941, p. 635. 


4263. Proposed by Howard Eves, Oregon State College, and Paul Halmos, 
Syracuse University 


Criticize the following alleged proof of the continuum hypothesis. 

Let X be the set of all infinite sequences of 0’s and 1’s, and let E be an 
arbitrary uncountable subset of X. Corresponding to any finite sequence, 

ay,+**, a}, of 0’s and 1’s, write E(a:, - - +, ax) for the set of all sequences 

xn} which belong to E and begin with {as, sey ax} . Since E= E(0)+E(1), at 
least one of the two sets E(0) and E(1) is uncountable; write a;=0 or 1 accord- 
ing as E(0) is or is not uncountable. Then, in either case, E(a;) is uncountable. 
If a; has already been defined for 7=1,---,k, so that E(a;, +--+, a) is un- 
countable, then write ai4,;=0, or 1 according as E(qi, - + - , ax, 0) is or is not 
uncountable. The resulting infinite sequence °° } has the property 
that for any k it is true that E(a, - - - , ax) is uncountable. Write E* for the 
union of all E(a, +--+, ax), for R=1, 2, 3,---+;then E* is a subset (in fact 
an uncountable subset) of E. 

For certain positive integers k it is true that both E(a, - - - , a, 0) and 
E(a, +++, a, 1) are uncountable; in fact this must happen for an infinite 
number of k’s. (Otherwise, for a sufficiently large k, E(ai, - - - , ax) would not 
be uncountable, contrary to its construction.) Let ki, ke, ks, - - - be the integers 
for which this is true, and write, for any {21, Xe, X3,° °° } in E*, 


Vn = 


then {y1, ye, - + - } is an infinite sequence of 0’s and 1’s. From the way in which 
the k, are defined it follows that every possible sequence of 0’s and 1’s occurs 
as a y sequence, and that consequently the sequences {x1, x2, - - - } in E* cor- 
respond (in possibly a many to one manner) to a set (viz. the set of all y se- 
quences) having the power of the continuum. It follows that the cardinal number 
of E* (and hence of E) cannot be less, and since E is a subset of X it cannot be 
greater. In other words it has been proved that every uncountable subset of a set 
having the power of the continuum has also the power of the continuum. 


4248 [1947, 232], corrected. Proposed by Victor Thébault, Tennie, Sarthe, 
France 


Having given a tetrahedron ABCD, place a sphere (S) of given radius in 
such a manner that the volume of the polar tetrahedron of ABCD with respect 
to (S) will be a relative minimum. 
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SOLUTIONS 
A Configuration of n Spheres 
3895 [1938, 631]. Proposed by Victor Thébault, Tennie, Sarthe, France 


Four spheres (w:), (we), (ws), (ws) pass respectively through the vertices 
Ai, Ao, As, A, of a tetrahedron and intersect in pairs on the corresponding edges. 
Straight lines AsAi, AsA2, AsA3, AsA4 which join an arbitrary point As to the 
vertices of the tetrahedron cut the respective spheres again in aj, a2, a3, a4 and 
these four points with As lie on a sphere (ws). Continuing in this way, the lines 
AnA1, AnA2, +++, AnAn—1 which join an arbitrary point A, to the points of a 
preceding set A2, +++, Ant cut the respective spheres (w1), (we), , (Wn—1) 
in the points a, a, - + + , @n-1 Which with A, lie on a sphere (w,). (1) Show that 
the spheres (w:), (we), + ++, (wn) meet in a point P. (2) Show that the points 
diametrically opposite to Ai, As, - - +, A, on the corresponding spheres lie in 
a plane through P. 


Note by Proposer. The part (2) is an extension of a theorem by Bouvaist- 
Delens, Annales de la Société Scientifique de Bruxelles, 1937, p. 159. 


Solution by the Proposer.* (1). A solution of this first part has already been 
given for n=5 (this MONTHLY, 3819 [1939, 239]). The same reasoning extended 
by adding one more point at a time shows that the spheres (w:), (we), - + - , (wn) 
meet in a point P common to the first four spheres. 

(2) To solve the second part it is convenient to establish a preliminary 


THEOREM. [f three spheres (w:), (we), (ws) pass respectively through the vertices 
A, Ao, As of a triangle A\A2As3, intersect in pairs on A1A2, A2A3, AsA1, and have 
in common a point P not in the plane of triangle A1A2A3, then the point P and the 
points Ai, Ad, Ag diametrically opposite Az, As on the given spheres are 
coplanar. 

The plane A;A2A; cuts the spheres (w:), (we), (ws) in three circles (Fi), (F2), 
(E3), which have a common point The points Az, Aj diametrically 
opposite Ai, As, As on (w;), (we), project orthogonally on the plane 
into points aj, az, a3 diametrically opposite A1, Az, A3 on (F:), (£2), (Es). The 
triangle aj a/aj has its sides perpendicular to the corresponding sides of triangle 
A,A2A3. The sides A? Aj, Aj Ai, Ai A? of triangle Ay Aj Aj are therefore per- 
pendicular to the homologous sides A2A3, A3A1, A1A2 of triangle A;A2A3. 

Likewise the orthogonal projections aj’, az’, a3’ of Ai, Az, A3 on the plane 
Aj A? Ag coincide with the points diametrically opposite A{, Az, Aj on the 
circles (Ej), (Ez), (Ej) cut by this plane from the spheres (w1), (w2), (ws). The 
sides of ai’ az’ az’, perpendicular to the corresponding sides of Aj A? Aj, inter- 
sect (E?/), (EJ) on AJ Aj, Aj Ai, A{ Ad. These three circles meet in the 
point P, which is the second point of intersection of the spheres (w1), (ws), (ws) 


* Translation by W. E. Byrne, Virginia Military Institute. 
¢ R. A. Johnson, Modern Geometry, 1929, Art. 184. 
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since P is not in the plane A;A2A3. The theorem is thus proved. 

The solution of (2) is merely a corollary. By the Miquel-Roberts theorem* 
the four spheres (w:), (we), (ws), (ws) of the present problem have in common 
a point P. The points Aj As, Aj diametrically opposite Ai, As, As, A, on 
(w1), (we), (ws), (ws) are in the same plane (m) passing through P, since the 
points (Aj, Ad, As, P), (Az, As, Ad, P), (As, Ad, Ai, P), (Ad, Ai, Az, P) 
are coplanar by our preliminary theorem. 

The lines AsA;, AsA2, AsA3, intersect (w:), (we), (ws), (ws) respectively 
in a1, @2, @3, as. By the Miquel-Roberts theorem applied to the tetrahedron 
A;A2A3Au4, the spheres Asc2aza4, (we), (ws), (ws) meet at P. Hence the point Ag 
diametrically opposite As on the sphere (ws) =Asa2a30%4 is in the plane (7). Step 
by step it is clear that the points Ag,---,A, diametrically opposite to 
As, On (ws), +++, (Wn) are also in the plane (7). 


Note. Part (2) generalizes an interesting theorem due to R. Bouvaist, who 
obtained it by analytical means. It was proved also in various ways by P. 
Delens.t The reasoning may be extended step by step to the polygon Aj Az -: - 
A, to prove part (2). The theorem is: 

Given a tetrahedron A,;A2A3A,4 and four spheres (a), (we), (ws), (ws) passing 
respectively through A, Ao, Az, As and intersecting in pairs on the corresponding 
adjacent edges , these spheres have a point P in common (S. Roberts), and the points 
Ai, Az, Az, Ad, diametrically opposite to Ai, Az, As, As on the spheres (w), 
(we), (ws), (ws) are in a plane passing through P. 


A Theorem on the Tetrahedron 
4150 [1945, 102]. Proposed by Victor Thébault, Tennie, Sarthe, France 


In a tetrahedron ABCD with the orthocenter H, the perpendiculars at A 
to the faces ACD, ABD, ABC meet respectively the planes HCD, HBD, HBC 
in (A;, Ae, As), and similarly for the points (B,, Bz, B;), and so on. Show that 
the planes A,A2A3, B,B2B3, --- , are perpendicular to the medians of ABCD. 


A solution will be found in the Proposer’s article, Tetrahedrons having a com- 
mon face, p. 395, this issue of the MONTHLY. 


Barycentric Cotrdinates 
4192 [1946, 103]. Proposed by Victor Thébault, Tennie, Sarthe, France 


For a given tetrahedron A1A2A3A,, if P is the centroid of its antipedal tetra- 
hedron ABCD, its barycentric coérdinates are inversely proportional to the 
squares of its distances to Ai, Az, A3, Aa; and conversely. 


* J. L. Coolidge, A Treatise on the Circle and the Sphere, 1916, p. 240. The priority of the 
theorem usually attributed to S. Roberts, seems to belong to Miquel, who gave an elementary 
demonstration of it in 1838 (Journal de Liouville, t. III, p. 522, theorem VI). Miquel obtained it 
by a stereographic projection. 

t P. Delens, Ann. de la Soc. Scient. de Bruxelles, 1937, p. 159. 
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Solution by R. Goormaghtigh, Bruges, Belgium. The following generalization 
will be proved: Let A1A2A3A, be a tetrahedron and B,B,B3B, the antipedal tetra- 
hedron of an arbitrary point P; if a1, a2, a3, a, and Bi, Be, Bs, Bs are the barycentric 
coérdinates of P as to and as to B\B,B;B,, then the ratios B2/a2, 
B3/as, B4/as are proportional to the squares of the distances PA;, PA2, PA3, PAg. 

The volume of the tetrahedron PA2A3A, is 


sin (P A2A3A4). 


Further, the areas o1,--- of the faces B,B;By,-+-+ are proportional to the 
sines of the supplements to the solid angles (B1— 22B;B,), - - - . But the solid 
angle (P—A2A3A,4) is the supplement to (B:—B,B3B,), and so on. Hence 
a1, +++ are proportional to o:i/PA:, and to - - - ; and this proves 
the theorem. The Proposer’s theorem is an immediate corollary. 

It follows also that the products a161, a282, a383, a484 are proportional to the 
squares of the faces of the tetrahedron B,B.B;By,. 


Collinear Orthopoles 
4195 [1946, 160]. Proposed by R. Goormaghtigh, Bruges, Belgium 


There are ten ways to divide six points on a circle into two groups of three 
so as to form pairs of triangles having no common vertex. The midpoints of the 
segments joining the orthopoles of a given straight line with respect to each pair 
of triangles are ten collinear points. 


Solution by Ou Li, Yenching University, Peiping, China. Let the given circle 
be the unit circle and let ¢; (¢=1, 2, - - - , 6) be the complex coérdinates of the 
six points A; on it. Denote the elementary symmetric functions of, say, ti, fe, fs 
by sj)” (j7=1, 2, 3); and those of th, f2, - + +, t by ox (k=1, 2, +++, 6). The given 
straight line has an equation in the form 


—+—-=1 
a a 


where a is constant. The orthopoles of this line with respect to the triangles 
A1A2A3, are, respectively ,* 


1 a 1 a 
2 a 2 a 
The midpoint of these two orthopoles will be 


1 a 
M = — [20 — 


Now let M’ be the analogous midpoint for the triangles A1A42A4 and A3;A5A¢. 


*See for instance R. Goormaghtigh, A Study of a Quadrilateral Inscribed in a Circle, this 
MontTaiy 1942, p. 178. 
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The clinant of the line MM’ becomes 


M — M’ (: 
=—= 
M — M' a 
Since this is symmetric in the ¢;, it is evident that all the ten midpoints are col- 


linear. 
Solved also by H. E. Fettis and the Proposer. 


Isosceles n-Points 
4254 [1947, 345]. Proposed by Paul Erdés, Syracuse University 


We have seven points in the plane. Prove that we can always select three 
which do not form an isosceles triangle. For six points this does not necessarily 
hold. (If A, B, C are on a line we can define that they do not form an isosceles 
triangle if AB#BC.) 


Editorial Note. This is essentially the same as the Proposer’s problem 
E 735 [1947, 227-229], which see for solution and comments. 


RECENT PUBLICATIONS 


EpiTep By H. P. Evans, University of Wisconsin 


All books for review should be sent directly to the editor of this department, American 
Mathematical Monthly, 531 West 116th Street, New York 27, N. Y., and not to any of the 
other editors or officers of the Association. 

An Introduction to Mathematical Genetics. By Lancelot Hogben. New York, 
W. W. Norton and Co., 1946. 12+260 pages. $5.00 


Genetical theories and various breeding practices lead to pleasing mathe- 
matical problems. They have intrigued abstract algebraists to construct special 
“genetic algebras,” and some of the recurrence relations and functional equa- 
tions describing the genetical structure of succeeding generations have not yet 
been solved. Readers of this MONTHLY will be more interested in the fact that 
genetical theories provide perhaps the easiest and best accessible examples il- 
lustrating how even very elementary mathematical reasoning leads to interesting 
insights and applications. It is a fairly common practice to use the Mendelian 
laws as a first application for the elementary rules of combining probabilities. 
The equally elementary and intuitive applications of simple difference and 
other recurrence equations, of matrix computation, and so forth, are less known. 
Mathematicians, and in particular mathematics teachers, will, therefore, greet 
a book on mathematical genetics with great interest. 

Unfortunately, the book under review will disappoint them inasmuch as it 
is clearly not intended for them. It is written for students of biology and practical 
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breeders who have difficulties with simple fractions but wish to understand the 
results of a quantitative analysis of breeding systems, of mutations, selection, 
and so forth. Accordingly, the choice of material, though natural and justified, 
is not ideal from the peculiar point of view of the mathematician. The mathe- 
matical tools and ideas are hidden so far as possible and the access to the book is 
made difficult by the standard practice in genetics of using as many technical 
terms as possible and not defining a single one. It is possible that students of 
biology would be acquainted with these terms long before they desire to con- 
sult Hogben’s book. Nevertheless, one page of explanations would have opened 
the book to many other readers. 

As has been indicated, the book relies on most elementary methods. An 
integral, a differential, and a statistical test of significance occur once or twice, 
but in a way calculated not to disturb a reader familiar only with simple algebra. 
All probability considerations are translated into the language of mixing and 
selecting cards. The notion of a sequence u, is carefully analyzed, but it may 
mislead some readers that u, usually refers simultaneously to any generation 
and to the initial or zero generation. In practice, uo is arbitrary and u,(m>0) is 
determined from a recurrence relation. Hence, identifying u, with uo is not 
legitimate. 

The content is as follows. Chapter I, Gene Frequencies, Genotypic Fre- 
quencies, and Systems of Mating. This chapter is mostly concerned with the 
card-pack model and the simplest rules. Chapter II, Basic Types of Algebraic 
Series in Genetical Theory, dedicated to recurrence relations, and limiting 
values. Chapter III, First Steps in the Calculus of Finite Differences, including 
approximate extrapolation. Chapter IV, Binomial Series (by which the author 
means the binomial probability distribution). Chapter V, Non-assortative Mat- 
ing in the Absence of Selection or Mutation. The term non-assortative stands 
for the usual “random” to which the author objects. This chapter treats the in- 
heritance of one or several characters, sex-linked or not. Chapter VI, Selection. 
Chapter VII, Assortative Mating and Consanguinity, covering various systems 
of breeding. Chapter VIII, Mutation Pressure and Isolate Effects. The Ap- 
pendices describe the simplest statistical techniques: I, Significance Tests for 
Mendelian Ratios, and II, The Estimation of Linkage and Determination of 
Variance Formulas for Gene-Frequency Analysis by the Method of Maximum 
Likelihood. 

W. FELLER 


Die Mathematische Denkweise. By Andreas Speiser. Basel, Birkhauser, 1945. 
132 pages. S. Fr. 14.50. 


In this book we have the philosophy of a mathematician. It is written with 
the enthusiasm of a distinguished mathematician who penetrates the arts and 
the world in his peculiar way. It will transmit, I imagine, this enthusiasm to 
every mathematician who is not only a craftsman but possessed by the sacred 
fire as the poet and philosopher ought to be. 
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The first section shows the réle of symmetry in ornamentics. Examples with 
beautiful pictures are given from Egypt and Crete, from Arabic and Renais- 
sance art. 

The following section gives a description of the complex of motives in a 
piece of music with characterization of the motive by the number of its meas- 
ures. The author gives examples from music by Mozart, Beethoven, Verdi and 
other composers. 

Very impressive is the third section, in which the author investigates Dante’s 
Divine Comedy in regard to natural philosophy which is essentially new 
Platonic. Here one finds a cosmic building which is the result of very scanty 
experience and a construction guided by the paramount desire for symmetry 
and harmony. If we compare this with the natural philosophy of today, we 
find the same tendency; however, the substructure of experience has grown to 
an immense extent and has so many interrelations that it represents in itself 
already a gigantic design. 

The fourth part is about the philosophy of mathematics by Proclus Dia- 
dochus. This philosophy is essentially based on Plato’s idea that mathematics, 
liberated from all the little devices of the craft, is a great ornament of divine 
origin. At the end of this section we find the following hymnic words of Proclus: 

“This, therefore, is Mathematics: she reminds you of the invisible forms of 
the soul; she gives life to her own discoveries, she awakes the mind and purifies 
the intellect; she brings to light our intrinsic ideas, she abolishes oblivion and 
ignorance which are ours by birth... . ” 

The average mathematical reader will, probably, appreciate only with diffi- 
culty the chapter on numbers and space in the Neoplatonic doctrine. As the 
author justly remarks, numbers were something much more alive, much more 
real with the Greeks than with us. One may, perhaps, add that in the develop- 
ment of human thought the concept of the general number is a rather late 
achievement and that the Greeks were close enough to the birth of this concept 
to feel the wonder of its nature. In these theories as in the cabala the numbers 
are not used to describe the order of the universe. They are felt as creating this 
order. In trying to penetrate this realm of thought a mathematician may have 
the feeling that he comes out of the clear day into an almost uncanny and 
mystical dark. 

The section on Goethe’s Theory of Colors follows. The author shows that 
Goethe tries to bring a system of symmetry into the realm of colors. In this way, 
through an apparent psychical order, he achieves a connection of the physical 
phenomena with the realm of beauty and art. 

The last section deals with astrology. Here we are introduced into the at- 
tempt to analogize the mathematical order of the stars with the somewhat 
chaotic fate of human beings and human society. The remarkable thing is that 
from the order of the heavenly bodies at the time of the birth, that is, from the 
horoscope, the great astrologers constructed human characters of a peculiar 
force. Sometimes, this construction formed the personality and life of the man 
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for whom the horoscope was cast. 

The book has an epilogue, consisting of an oration which Speiser gave at 
the tercentenary of Kepler’s death. Here, it is especially important to be re- 
minded of the fact that Kepler’s firm belief in the mathematical structure of the 
universe showed him the way and strengthened his amazing power of concen- 
tration and endurance to achieve his immortal work on the orbits of the planets. 

It is interesting to realize that Speiser, in his discussions, dwells mainly on 
long bygone times, when the physical world and the world of the human soul 
were felt, by the wise men as well as by the common people, as one. 

M. DEHN 


Concise Analytic Geometry. By C. H. Sisam, New York, Henry Holt and Co., 
1946. 9+155 pages. $2.00. 


This book comprises the section on analytic geometry contained in the au- 
thor’s more extensive College Mathematics, which was reviewed in the February, 
1947 issue of this MONTHLY (volume 54, number 2, page 119). 

The definitions and proofs of theorems have been formulated with care and 
accuracy and no essential topic has been omitted. In spite of the concise treat- 
ment the explanations, illustrative examples, and exercises are quite adequate. 
Polar coérdinates are introduced early and the polar equations of straight lines, 
circles, and conics are derived. The application of polar coédrdinates to more 
complicated loci is deferred to a later chapter. The last three chapters provide a 
brief introduction to solid analytic geometry. Answers to the odd numbered 
exercises are given in the back of the book. 

In the opinion of the reviewer the book contains ample material for a three 
semester hour course in Analytic Geometry and it is believed that instructors 
who have experienced difficulty in finding a suitable text for such a course would 
do well to examine this book carefully. 

H. P. Evans 


New Books Received 


Antennae: An Introduction to Their Theory. By J. Aharoni. New York, 
Oxford University Press, 1946. $8.50. 

College Algebra. Third Edition. By W. L. Hart. Boston, D. C. Heath and 
Co., 1947. 8+416 pages. $3.00. 

College Algebra. Alternate Edition. By Paul Rider. New York, The Mac- 
millan Company, 1947. 15+407 pages. $2.50. 

Engineering Problems Manual. Fourth Edition. By F. C. Dana and L. R. 
Hillyard. New York and London, McGraw-Hill Book Co., 1947. 14+419 pages. 
$3.25 

Mathematics of Accounting. Third Edition, By A. B. Curtis and J. H. 
Cooper. New York, Prentice-Hall, Inc., 1947. 10+550 pages. $6.00 trade edi- 
tion, $4.50 text edition. 

Preparatory Business Mathematics. By L. L. Smail. New York, Ronald 
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Press Co., 1947. 10+244 pages. $2.75. 

Multiple-factor Analysis. By L. L. Thurstone. University of Chicago Press, 
1947. 19+535 pages. $7.50. 

Plane Trigonometry. By E. B. Mode. New York, Prentice-Hall, Inc., 1947. 
10+216 pages. 

Relativity. The Special and General Theory. By Albert Einstein. New York, 
Henry Holt and Co., 1930. Reprint, New York, Peter Smith, 1947. 13+168 
pages. $2.50. 

Lecciones de Analisis Infinitesimal. By F. J. Duarte. Caracas, Tipografia 
Americana, 1943. 605 pages. 

An Introduction to Business Statistics. Second Edition. By J. R. Stockton. 
Boston, D. C. Heath and Co., 1947. 7+478 pages. $4.00. 

Six-Place Tables. Seventh Edition. By E. S. Allen. New York, McGraw-Hill 
Book Co., 1947. 23+232 pages. $2.50. 


CLUBS AND ALLIED ACTIVITIES 


EpitTep sy L. F. OLLMANN, Hofstra College 


Send reports of all activities, such as club reports, special features, topics with references, 
student papers, and other material of interest to L. F. Ollmann, Hofstra College, Hempstead, 
New York. 


CLUB REPORTS, 1946-47 
Mathematics Club, Berea College 


Professor W. R. Hutcherson of Berea College, was in charge of a mathematics 
and astronomy exhibit in connection with a science exhibit held at Berea, Ken- 
tucky on March 29, 1947. Many favorable remarks were made by visitors who 
classed the display in mathematics ahead of those of some of the sciences. Ex- 
tensive use of placards, maps, blackboards and models was made in the follow- 
ing exhibits: 

History of Mathematics: Development and trends 

Computing Machines: The mechanical brain 

Statistics: Applications and graphs (bar graphs, curves) 

Mechanical Drawing: Blueprinting and layout work 

Mathematics and Music 

Slide Rule: Uses and types 

Mathematics Tables and Applications of Calculus 

Surveying: Instruments, and maps 

Astronomy: Telescope, relative star sizes, and pictures 

Magic in Mathematics: Magic squares, perpetual calendar, and puzzles 

Modern Geometry: Law of duality and Desargues’ Theorem 

Calculus: Mechanics, and friction problems 
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Members of the club were divided into groups according to their interests 
and each group had charge of planning one of the exhibits and then demonstrated 
it during the period for exhibition. 


Mathematics Club, Rutgers University 


The Mathematics Club of Rutgers University, which had suspended its 
activity during the war, resumed operation on December 5, 1946, and regular 
semi-monthly meetings were held during the balance of the school year. 

The following papers were presented: 

Farey series, by Professor H. S. Grant 

Pythagorean triangles, by Professor C. R. Phelps 

Maxima and minima, by Dr. L. M. Court 

The actuarial profession, by Professor R. Walter 

Constructibility, by George Y. Cherlin 

Mathematics in industry, by Professor F. G. Fender 

The officers for 1946-47 were: President, George Y. Cherlin; Vice-President, 
Robert Coursen; Secretary-Treasurer, C. F. Pinzka. The Faculty Advisor was 
Professor C. R. Phelps. 


Kappa Mu Epsilon, William Jewell College 

Mu Sigma Alpha, the mathematics honorary fraternity at William Jewell 
College was organized by a group of interested students in January, 1943. It 
has assumed a place of distinction and respect among the honor societies on 
the campus. During the year 1946-47, the programs consisted of the following 
speeches made by members of the club and guest speakers: 

Contribution of Einstein toward the atomic bomb and his mass-velocity rela- 
tionship, by Dr. W. A. Hilton 

The invention and development of the concept of zero, by Professor L. O. Jones 

Fallacies in mathematics, by Nicholas Housley 

Short-cuts of mathematics, by Professor C. O. Van Dyke 

Geometric designs in lighting, by LeRoy Heaton 

Telling time by the stars, by Lloyd Elrod 

Improvement of high school mathematics courses, by Paul Curau 

The case for mathematics, by Woody Rixey. 

During the year 1946-47, several of the students became interested in be- 
coming affiliated with a national mathematical fraternity. Kappa Mu Epsilon 
was petitioned for membership, and on May 7, 1947, the Missouri Gamma 
Chapter of Kappa Mu Epsilon was duly installed at William Jewell College by 
the Missouri Beta Chapter from Warrensburg under the direction of Dr. Claude 
Brown. The officers for 1946-47 were: President, Mary Ruth Carney; Vice- 
President, Truett Neese; Treasurer, LeRoy Heaton; Corresponding Secretary, 
Professor L. O. Jones; Secretary, Edwin Watson. 

The officers elected for the year 1947-48 are: President, LeRoy Heaton; 
Vice-President, Maynard Cowan; Secretary, Walter Binns; Corresponding 
Secretary, Paul Swedburg; Treasurer, Lloyd Elrod. 
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Mathematics Club, Haverford College 


During the past academic year, the following papers were prepared and 
informally discussed : 

Elements of topology, by Dr. C. B. Allendoerfer 

Squaring the circle, by Geert Prins 

The four-color problem, by James Thorpe 

Fermat's last theorem, by Murray Freeman 

Fibonacci’s series, by Norman Brous 

Ham sandwich theorems, by Dr. John W. Tukey of Princeton University 

The cattle problem of Archimedes, by Daniel Wagner. 

The members of the club who distinguished themselves by presenting solu- 
tions for the monthly prize problems were: Murray Freeman, William Warner, 
Daniel Wagner, and Thomas Crolius. 

The officers of the club for the year 1946-47 were: President, Geert Prins; 
Secretary, Leon Robbins, Jr. 


Pi Mu Epsilon, St. Louis University 


The Missouri Gamma Chapter of Pi Mu Epsilon has evidenced a very 
strong interest in the promotion of scholarship. At its three regular meetings, 
members of the chapter contributed scholarly papers on the following subjects: 

History of logarithms, by Miss Margaret Willerding 

Transcendental numbers, by Mr. Edwin Karlowicz 

Mechanism and use of integrating machines, by Mr. Herman Plew, Jr. 

A summary of each of these papers will appear in the regular issue of the 
Chapter news magazine, The Missouri Gamma News. Each of the meetings was 
followed by a social hour. 

The principal and final function of the year took place on April 8 at Parks 
Air College. It consisted in the initiation of new members, the election of the 
new Director, Mr. Oliver F. Anderhalter, the Annual Lecture, and the Tenth 
Annual Banquet. The lecturer, Professor N. A. Court of the University of Okla- 
homa, spoke on Mathematical Asides. Rev. Patrick J. Holloran, S. J., President 
of Saint Louis University, addressed the group at the banquet. He expressed 
his personal satisfaction at the activity and growth of Pi Mu Epsilon. He gave 
a vote of thanks to Professor Regan for being the inspiration for much of this 
development. 

Missouri Gamma News is a Chapter bulletin which is now in its second year 
of publication. It reviews the principal activities of the Chapter during the year 
and the plans for the year ahead. A summary by the author of the paper read 
at the regular meetings is included in the publication. News items about mem- 
bers of the Chapter, the plans of the Mathematics Department and of the allied 
sciences can be found in the Chapter news. 

Officers for the year 1946-47 were: Director, Mr. Herbert Gebhart; Vice- 
Director, Miss Margaret Willerding; Secretarv, Miss Frances Higgins; Faculty 
Advisor and Corresponding Secretary, Professor Francis Regan. 
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NEWS AND NOTICES 


EpiTep By HARRY POLLARD, Cornell University 


Readers are invited to contribute to the general interest of this department by sending news 
items to Harry Pollard, White Hall, Cornell University, Ithaca, New York. 


STANFORD UNIVERSITY MATHEMATICS EXAMINATION 


The second Stanford University Mathematics Examination (see this 
Monta .y, vol. LIII, no. 7, pp. 406-409 (1946)) was held April 19, 1947, 2:00 
p.m. to 5:00 p.m., in thirty-nine high schools in California; 196 students took 
part. The following problems were proposed. 

1. To number the pages of a bulky volume the printer used 1890 digits. How many 
pages has the volume? 
2. Among grandfather’s papers a bill was found: 


72 turkeys $—67.9— 


The first and last digit of the number that obviously represented the total price of 
those fowls are replaced here by blanks, for they have faded and are now illegible. 

What are the two faded digits and what was the price of one turkey? 

3. Determine m so that the equation in x 


— Gn +0? 0 


has four real roots in arithmetic progression. 
4. Let a, Band vy denote the angles of a triangle. Show that 


sin a + sin + sin = 4 cos cos 


sin 2a + sin 28 + sin 2y = 4sin asin B sin y 
and 
sin 4a + sin 48 + sin 4y = — 4sin 2a sin 28 sin 2y. 


The writer of the best paper, Miss Maxine Steineke, Los Altos, California, 
a student at Palo Alto High School, Palo Alto, California, received a $500 
scholarship at Stanford University. The following contestants received “Hon- 
orable Mention”: 
Mr. Raymond M. Hendricks, Santa Barbara, California, a student at 
Santa Barbara High School. 
Mr. Jerold Bernard Weller, Los Angeles, California, student at Alexander 
Hamilton High School, Los Angeles. 


PERSONAL ITEMS 

Paul Alexandroff has been elected a foreign member of the National Acad- 
emy of Sciences. 

K. J. Arrow of Columbia University has been appointed research associate 
with the Cowles Commission for Research in Economics at the University of 
Chicago. 

Professor W. L. Ayres of Purdue University has been appointed Dean of the 


430 


| 
| 


1947] NEWS AND NOTICES 431 


School of Science. 

Dr. L. M. Brown of the University of Edinburgh, Dr. J. M. Hyslop of the 
University of Glasgow, Dr. J. M. Jackson of University College, Dundee, and 
Dr. A. J. Macintyre of the University of Aberdeen have been elected fellows of 
the Royal Society of Edinburgh. 

Dean R. D. Carmichael of the Graduate School, University of Illinois, has 
retired with the title emeritus. 

Miss Mary L. Cartwright of the University of Cambridge has been elected 
a fellow of the Royal Society. 

Professor T. G. Cowling of the University College of North Wales has been 
elected a fellow of the Royal Society. 

Associate Professor J. B. D. Derksen of the Netherlands School of Eco- 
nomics at Rotterdam has joined the Statistical Office of the United Nations at 
Lake Success. 

Dr. Mary P. Dolciani of Cornell University has been awarded the Sigma 
Delta Epsilon post doctoral fellowship for 1947-48 and will study at the In- 
stitute for Advanced Study. 

T. K. Glennan of the Ansco Division of General Aniline and Film Corpora- 
tion, Binghamton, New York, has been appointed President of the Case In- 
stitute of Technology. He was wartime director of the U. S. Navy Underwater 
Sound Laboratory. 

Professor C. C. MacDuffee of the University of Wisconsin has been awarded 
the degree of Doctor of Science at Colgate University. 

Professor O. E. Neugebauer of Brown University and Professor Hassler 
Whitney of Harvard University have been elected to membership in the Ameri- 
can Philosophical Society. 

Messrs. Paul Olum of Harvard University and E. H. Spanier of the Uni- 
versity of Michigan have been awarded Frank B. Jewett fellowships for 1947-— 
48 by the American Telephone and Telegraph Company. 

Professor P. A. Smith of Columbia University has been elected to member- 
ship in the National Academy of Sciences. 

Associate Professor Max Astrachan of Antioch College has been promoted to 
a professorship. 

Dr. W. D. Baten is now Chief, Operations Analysis Branch (A-5), Air De- 
fense Command, Mitchell Field, New York. 

Professor C. C. Bramble of the Naval Postgraduate School, Annapolis, 
Maryland, has accepted a position as Director of Computation and Ballistics 
at the Naval Proving Ground, Dahlgren, Virginia. 

Professor Talmon Bell of Sterling College has retired as chairman of the 
Mathematics Department with the title of Professor Emeritus. 

Assistant Professor Lipman Bers of Syracuse University has been promoted 
to an associate professorship. 

Professor Garrett Birkhoff of Harvard University served as Walker Ames 
Professor at the University of Washington during the spring quarter. 


| 
i 
a 
| 


432 NEWS AND NOTICES [September, 


Associate Professor I. S. Carroll of Syracuse University has been promoted 
to a professorship. 

Dr. Myrtle Collier, Chairman of the Department of Mathematics of Im- 
maculate Heart College, has retired with the title of Professor Emeritus. 

Dr. John Dyer-Bennet of Purdue University has been promoted to an as- 
sistant professorship. 

Professor Samuel Eilenberg of Indiana University has been appointed to a 
professorship at Columbia University. 

Dr. Cleota G. Fry of Purdue University has been promoted to an assistant 
professorship. 

Assistant Professor Abe Gelbart of Syracuse University has been promoted 
to an associate professorship. He will be on leave for the academic year 1947-48 
to attend the Institute for Advanced Study. 

Dr. W. H. Gottschalk of the University of Pennsylvania has been promoted 
to an assistant professorship. 

Associate Professor S. H. Gould ofthe University of Toronto has been ap- 
pointed to an assistant professorship at Purdue University. 

Assistant Professor May N. Harwood of Syracuse University has been pro- 
moted to an associate professorship. 

Associate Professor M. H. Heins of Brown University has been promoted 
to a professorship. 

Associate Professor Magnus R. Hestenes of the University of Chicago has 
been appointed to a professorship at the University of California at Los Angeles. 

Associate Professor H. K. Hughes of Purdue University has been promoted 
to a professorship. 

Assistant Professor Mark Kac of Cornell University has been promoted to a 
professorship. 

Dr. Irving Kaplansky of the University of Chicago has been promoted to 
an assistant professorship. 

Assistant Professor M. W. Keller of Purdue University has been promoted 
to an associate professorship. 

Assistant Professor J. L. Kelley of the University of Chicago has been 
appointed to an associate professorship at the University of California at 
Berkeley. 

Dr. Carl Kossack of the Navy Department has been appointed to an asso- 
ciate professorship at Purdue University. 

R. R. Kuebler of Dickinson College, Carlisle, Pennsylvania, me been pro- 
moted to an assistant professorship. 

Professor C. G. Latimer of the University of Kentucky has sii apponted 
to a professorship at Emory University. 

Professor L. L. Lowenstein of Alfred University has been appointed professor 
and head of the department at Kent State University. 

Professor Saunders MacLane of Harvard University has been appointed to 
a professorship at the University of Chicago. He will be on leave as a Guggen- 
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heim Fellow until the summer quarter of 1948. 

G. E. Markle of the University of Detroit has been promoted to an assistant 
professorship. 

R. D. Mindlin of Columbia University has been promoted to a professorship. 

Associate Professor Harriet F. Montague of the University of Buffalo has 
been promoted to a professorship. 

M. G. Moore of Bradley University, Peoria, Illinois, has been promoted to 
an associate professorship. 

Assistant Professor Andrewa R. Noble of Montana State University has 
been appointed to an associate professorship at Pacific University, Forest 
Grove, Oregon. 

Nilan Norris of Hunter College has been promoted to an assistant pro- 
fessorship. 

J. I. Northam of the University of Wisconsin has been appointed to an 
assistant professorship at Kansas State College, Manhattan, Kansas. 

Edward Paulson has been appointed lecturer in mathematics at the Univer- 
sity of Washington. 

Assistant Professor A. M. Peiser of Rutgers University has resigned to ac- 
cept a position as mathematician with Hydrocarbon Research, Inc., of New 
York. 

Dr. Clarence Perisho of McCook Junior College has been appointed assistant 
professor of chemistry and mathematics at Nebraska Wesleyan University, 
Lincoln, Nebraska. 

Dr. Sam Perlis of Purdue University has been promoted to an assistant pro- 
fessorship. 

Professor H. B. Phillips of the Massachusetts Institute of Technology has 
retired. 

Assistant Professor R. S. Phillips of New York University has been ap- 
pointed to an associate professorship at the University of Southern California. 

E. K. Ritter of the Postgraduate School, United States Naval Academy, has 
been promoted to an associate professorship. 

Assistant Professor C. K. Robbins of Purdue University has been promoted 
to an associate professorship. 

Associate Professor Arthur Rosenthal of the University of New Mexico has 
been appointed to a professorship at Purdue University. 

Professor A. C. Schaeffer of Stanford University has been appointed to a 
professorship at Purdue University. 

Dr. Edith R. Schneckenburger of Michigan State Normal College has been 
appointed to an assistant professorship at the University of Buffalo. 

Dr. I. E. Segal has been appointed to an assistant professorship at the 
University of Chicago. 

Assistant Professor Ruth G. Simon of Berea College has been appointed to 
an assistant professorship at Morningside College, Sioux City, lowa. 

Dr. Jerome C. Smith of Lafayette College, Easton, Pennsylvania, has been 
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promoted to an assistant professorship. 

Assistant Professor B. M. Stewart of Michigan State College has been pro- 
moted to an associate professorship. 

Professor A. H. Taub of the University of Washington has been granteda 
leave of absence to accept a Guggenheim Fellowship. He will spend the year at 
Princeton University. 

Bradford Tye of Bethany College has been promoted to an assistant pro- 
fessorship. 

Assistant Professor G. B. Van Schaack of Union College has been appointed 
to an assistant professorship at Washington University. 

Assistant Professor G. L. Walker of Temple University has been ap- 
pointed to an assistant professorship at Purdue University 

Assistant Professor A. D. Wallace of the University of Pennsylvania has 
been appointed to a professorship at Tulane University. 

Professor Andre Weil of the University of Sao Paulo has been appointed 
visiting professor at the University of Chicago. 

Associate Professor E. L. Welker of the University of Illinois has been ap- 
pointed associate in mathematics in the Bureau of Medical Economic Research 
of the American Medical Association. 

D. W. Western of Brown University has been promoted to an assistant 
professorship. 

Dr. G. W. Whitehead of Princeton University has been appointed to an 
assistant professorship at Brown University. 

Professor R. M. Winger of the University of Washington has been appointed 
Executive Officer of the Department of Mathematics. 

Dr. Hyman J. Zimmerberg of Rutgers University has been promoted to an 
assistant professorship. 

Assistant Professor H. S. Zuckerman of the University of Washington has 
been promoted to an associate professorship. 

Professor Antoni Zygmund of the University of Pennsylvania has been ap- 
pointed to a professorship at the University of Chicago. 


The following appointments to instructorships are announced: 
Brown University: Dr. A. B. Blakers, Mrs. R. C. Buck, Dr. Kathleen 
Butcher. 
Cornell University: Dr. Christine Williams, Dr. Bertram Yood, Dr. Bryant 
Tuckerman, Jr. 
General Motors Institute, Flint, Michigan: J. A. Straw 
Rutgers University: C. D. Firestone 
United States Naval Academy: J. R. Gorman 
University of California: Dr. E. W. Barankin 
University of Chicago: Edwin J. Akutowicz, J. Bruce Crabtree 
University of Washington: Dr. Fumio Yagi 
University of Wisconsin: Dr. A. M. Mark 
Wellesley College: Mrs. Jerry Cowen 
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C. C. Carter, Bluffs, Illinois, honorary life member of the Association, lost 
his life in an accident May 5, 1947. 

Professor P. J. Daniells of the University of Sheffield died May 25, 1946. 

Professor Emeritus C. J. Keyser of Columbia University died May 8, 1947 
at the age of eighty-six years. 

The death of M. A. D. Kinsman has been reported. 

Professor L. G. Owens, formerly of the University of Rangoon, died Febru- 
ary 23, 1947. 


THE MATHEMATICAL ASSOCIATION OF AMERICA 


Official Reports and Communications 
MEETING OF THE NORTHERN CALIFORNIA SECTION 


The ninth annual meeting of the Northern California Section was held at the 
University of San Francisco on Saturday, January 25, 1947. Professor W. H. 
Myers, Chairman of the Section, presided at the morning and afternoon ses- 
sions. 

The attendance was eighty-one including the following twenty-three mem- 
bers of the Association: H. M. Bacon; G. A. Baker; T. J. Bass, Jr.; Mable L. 
Beckwith; W. G. Brady; E. L. Fitzgerald, S.J.; S. A. Francis; Alice Graeber; 
W. R. Hanson; Emma V. Hesse; D. H. Lehmer; Sophia L. McDonald; E. D. 
Miller; F. R. Morris; W. H. Myers; C. D. Olds; George Pélya; Edris P. Rahn; 
E. B. Roessler; Pauline Sperry; Gabor Szegé; K. J. Waider; L. A. Walker. 

At the business meeting the following officers were elected for the coming 
year: Chairman, Professor George Pélya, Stanford University; Vice-Chairman, 
Professor G. C. Evans, University of California; Secretary-Treasurer, Professor 
E. B. Roessler, University of California at Davis; Representative on the Cali- 
fornia Journal of Secondary Education, Mrs. Ruth G. Sumner, Oakland High 
School. 

By invitation of the Section, Professor L. E. Reukema of the Department of 
Electrical Engineering, University of California, gave an hour’s address during 
the morning session. 

The following papers were presented: 

1. Vibration modes of tapered beams, by Dr. Edmund Pinney, University of 
California, introduced by the Secretary. 

The vibration modes of thin beams whose bending stiffness and linear densi- 
ties may be expressed in simple powers of the distance along the beam may, in 
certain cases, be obtained in terms of Bessel functions. Four one-parameter fami- 
lies of beams are shown to have this property. These include the wedge and the 
cone. The results may be applied to finding bending modes for airplane wings 
while still in the preliminary design stage. 
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2. Estimating electrostatic capacity from geometric data, by Professor George 
Pélya, Stanford University. 

Professor Pélya’s paper appeared in this MONTHLY, April, 1947. 

3. Electronics in our modern world, by Professor L. E. Reukema, University 
of California, introduced by the Chairman. 

The speaker recalled that electronics includes the science and industry con- 
nected with vacuum tubes and their circuits. The tubes may be the tiny ones 
found in radio receivers, or tubes capable of passing many thousands of amperes 
or of generating pulses of millions of watts of radio-frequency power. They 
may act as rectifiers, amplifiers, oscillators, modulators, or detectors. Without 
vacuum tubes, long-distance telephones, world-wide radio communication, 
sound motion pictures, television, and radar would all be impossible. Vacuum 
tubes also guide our planes through the skies, our ships over the sea, and our 
missiles to enemy targets. They allow our planes to land through the densest fog. 
They automatically aim and fire our guns. In medicine and surgery electronic 
equipment has made notable contributions. 

It is truly amazing what powers the vacuum tube gives to mankind. It makes 
possible the accurate measurement of distances as much smaller than the diame- 
ter of a raindrop as the diameter of the raindrop is smaller than that of a sphere 
50 miles across. It enables one to weigh particles too small to be seen with the 
finest microscope, to detect energy in space in such tiny quantities that it would 
have to be received continuously for over a year to amount to as much energy 
as a fly expends in raising itself one inch, to measure currents so small that they 
would have to flow for over a billion years to pass as much electricity as flows 
through the filament of an ordinary electric light in a second. It enables a single 
television transmitter to hold a million or more receivers in synchronism to 
within one sixteen millionth part of a second. In the next few years it will proba- 
bly supply jobs to millions of men in completely new industries, and will in 
thousands of ways increase our safety, comfort, and material wealth, and hap- 
piness. 

4. The impact of high speed computing on mathematics, by Professor D. H. 
Lehmer, University of California. 

This paper dealt with the several projects now under way in high speed 
computing devices, and the influences they are likely to have on mathematical 
problems and mathematical points of view. 

5. Refresher courses for returning students, by S. A. Francis, San Mateo Junior 
College. 

Mr. Francis discussed the topic of refresher courses as related to the field of 
mathematics. The speaker based his discussion on actual experiences with teach- 
ing such courses in classes under the auspices of E.S.M.W.T. during the war 
period from 1943 to 1945, and for returning students majoring in engineering 
and science at Stanford University and at San Mateo Junior College. Many 
sections of refresher courses were organized at San Jose State College, San Mateo 
Junior College, San Francisco Junior College, and at the University of California 
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to accommodate returning students who asked for such courses in the fields of 
algebra, trigonometry, analytic geometry, and calculus. The consensus of 
opinion of those engaged in the teaching of such classes is that the courses were 
of great value to the students as well as an encouraging and challenging experi- 
ence for teachers of mathematics. 

6. Some problems in aerodynamics, by Professor J. G. Herriot, Stanford Uni- 
versity, introduced by the Chairman. 

Professor Herriot gave a brief description of a wind tunnel, noting that the 
speed which can be attained is limited by the size of the test model even when 
sufficient power is available. In order to explain this he gave a simple mathe- 
matical proof of the fact that in subsonic flow a decrease in the cross-sectional 
area of the flow produces an increase in velocity whereas in supersonic flow the 
reverse is true. It was also shown how the non-linear differential equation of com- 
pressible flow can be linearized under the assumption of small perturbations of 
the velocity. It is then possible to deduce the properties of a compressible flow 
field from those of a known incompressible flow field. A number of other prob- 
lems of current interest were mentioned briefly. These included pressure dis- 
tributions over bodies of revolution in subsonic and supersonic flow, and mixed 
subsonic and supersonic flows. 

7. Some remarks on mathematical education, by Professor Sophia L. Mc- 
Donald, University of California. 

E. B. Secretary 


MEETING OF THE OKLAHOMA SECTION 


The annual meeting of the Oklahoma Section of the Mathematical Associa- 
tion of America was held in connection with the convention of the Oklahoma 
Education Association in Oklahoma City on Friday morning, February 14, 1947. 
Professor O. H. Hamilton, Chairman of the Section, presided. 

Sixty-two persons attended the meeting, including the following twenty-one 
members of the Association: E. F. Allen, Arthur Bernhart, J. C. Brixey, H. N. 
Carter, N. A. Court, A. H. Diamond, O. H. Hamilton, J. O. Hassler, W. N. Huff, 
H. V. Huneke, J. E. LaFon, P. E. Lewis, G. E. Meador, H. A. Palmer, W. L. 
Shepherd, W. T. Short, D. R. Shreve, H. W. Smith, C. E. Springer, R. W. 
Veatch, J. H. Zant. 

At the business session the following officers were elected: Chairman, D. R. 
Shreve, University of Tulsa; Vice-Chairman, W. T. Short, Oklahoma Baptist 
University; Secretary, J. C. Brixey, University of Oklahoma. 

The program consisted of the following ten papers: 

1. On the representation problem for projective algebras, by Professor J. C. C. 
McKinsey, Oklahoma A. and M. College, introduced by Professor A. H. 
Diamond. 

C. J. Everett and S. Ulam have recently given postulates for projective alge- 
bras, and have solved the representation problem for all complete atomic projec- 
tive algebras. The author showed that every projective algebra is isomorphic to 
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a sub-algebra of a complete atomic projective algebra; this result, combined 
with that of Everett and Ulam, provides a solution of the representation problem 
for arbitrary projective algebras. The method of proof is similar to the method 
used by M. H. Stone to establish a representation theorem for Boolean algebras. 
If A is a prime ideal, then A, is defined to be the prime ideal which contains all 
elements of the form a,z, for a in A. If A and B are prime ideals, then ADB 
is defined to be the class of all prime ideals C which contain all elements of the 
form ab, where a is in A and d is in B. These definitions were extended to per- 
mit operations on classes of prime ideals. 

2. Some properties of algebras and their subalgebras, by Professors A. H. 
Diamond and J. C. C. McKinsey, Oklahoma A. and M. College. 

The problem considered was that of sub-algebras, S, generated by n elements 
of an algebra A for a given value of n, such that every S is a Boolean algebra 
but A is not Boolean. A proof was given that is not less than 2. This result was 
used to prove a theorem of Wajsberg that every set of axioms for the sentential 
calculus must contain at least one axiom in at least three variables. 

3. Distance in the complex domain, by Professor W. T. Short, Oklahoma 
Baptist University. 

A point P in the complex domain was represented by the coérdinates 
(x+ut, y+vi). The distance between two such points Pi and P: was defined by 
d? = (%2—%1)? + (y2— 1)? + (U2 1)? + (¥2—01)?, where the real part r?=(x2—4%1)? 
+(y2—y1)?, and the normal part n?= (u2—11)?+(ve—1)?. The vector joining the 
two points is then r-+-n1, where 7? = —1. The vector is in the plane determined by 
the perpendiculars at P; and P;. The vector in the complex domain consists of 
two parts. The real part 7 is the resultant of the two real vectors. To find the 
normal vector n, we take the resultant of the two normal vectors in the real 
plane. This length normal to the real plane is the normal vector. 

4. The inverse geometry of Clifford numbers, by Professor E. F. Allen, Okla- 
homa A. and M. College. 

The number z=x-+ry, r?= +k’, k real, was used to solve geometric problems. 
In particular it was used to demonstrate a generalization of the concept of the 
perpendicularity of two lines. 

5. A contribution to the four-color problem, by Professor Arthur Bernhart, 
University of Oklahoma. 

Though a map may not be 4-colorable, yet each hemisphere of a minimal 
map must be colorable, subject to certain boundary conditions on the equatorial 
ring. For an n-ring there are [3"-'+3(—1)"+2]/8 color schemes. In the case 
of an 8-ring this leads to homogeneous equations in 274 unknowns. The alge- 
braic solution is facilitated by a magic square arrangement. Results include a 
new reducible configuration: two hexagons guarded by two pentagons. 

6. A comment on the approximate solution of certain boundary value problems 
by submatric inversion, by R. R. Reynolds, Oklahoma A. and M. College, intro- 
duced by Professor O. H. Hamilton. 

Matric analysis was applied to the approximate solution and the improve- 
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ment of the approximate solution of certain boundary value problems. 

7. Mathematical analysis of artificial pulse lines, by P. J. Miller, Oklahoma 
A. and M. College, introduced by Professor H. S. Smith. 

This was a preliminary report on work which is still in progress. 

8. A conic in the tangent plane to a surface, by W. W. Dolan, University of 
Oklahoma, introduced by Professor C. E. Springer. 

The polar reciprocal of any line / at a point P of a surface, with respect toa 
Darboux quadric at the point, is a line in the tangent plane. If / varies so as to 
generate a quadric cone, the reciprocal line envelops a conic. Two theorems ap- 
pear: (1) the reciprocal conic will be a hyperbola, parabola, or ellipse according 
as the line of centers of the Darboux quadrics is outside, on, or inside the cone 
generated by /; (2) the point P will be outside, on, or inside the reciprocal conic 
according as the cone generated by / has two lines, one line, or no line in common 
with the tangent plane. 

Analytic treatment in a properly restricted case reveals that the reciprocal 
conic is an ellipse whose axes are simple functions of the principal radii of normal 
curvature of the surface at the point. 

9. On a differential equation satisfied by a certain set of polynomials, by Pro- 
fessor W. N. Huff, University of Oklahoma. 

The set {yn(x) } of polynomials given by the generating function 


g(x, t) = f(xt)d(t) = 
with 


f(«t) = > 0. and ¢(t) = >> i=, (dn ¥ 0, bo ¥ 0), 


n=0 n=0 


was briefly considered. A relation for which these polynomials are of A type k 
in according with the definition of Sheffer was given. Finally a differential equa- 
tion of infinite order satisfied by the set { yn(x) | was obtained along with a neces- 
sary and sufficient condition that the differential equation be of finite order. 

10. Diagnostic testing and remedial instruction in freshman mathematics for 
engineers, by Professor D. R. Shreve, University of Tulsa. 

A report on a three-year program of diagnostic testing and remedial instruc- 
tion in freshman mathematics for engineers as carried out at Purdue University 
(1938-1942) was discussed. 

Joun C. Brixey, Secretary 


MARCH MEETING OF THE MICHIGAN SECTION 


The spring meeting of the Michigan Section of the Mathematical Associa- 
tion of America was held at the University of Michigan at Ann Arbor on Sat- 
day, March 22, 1947. This meeting also constituted the meeting of the Mathe- 
matics Section of the Michigan Academy of Science, Arts and Letters. Morning 
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and afternoon sessions and a luncheon-business meeting were held, at all of 
which the Chairman, Professor E. E. Ingalls, presided. 

About ninety persons attended the meeting including the following forty- 
seven members of the Association: N. H. Anning, J. W. Baldwin, J. M. Barbour, 
F. A. Beeler, W. M. Borgman, J. W. Bradshaw, J. B. Brandeberry, R. V. 
Churchill, Nathaniel Coburn, A. H. Copeland, C. C. Craig, Wayne Dancer, 
P. S. Dwyer, C. M. Erikson, C. H. Fischer, J. S. Frame, V. G. Grove, G. E. Hay, 
Fritz Herzog, T. H. Hildebrandt, E. E. Ingalls, L. G. Johnson, L. S. Johnston, 
P. S. Jones, Wilfred Kaplan, G. P. Loweke, E. D. McCarthy, L. E. Mehlen- 
bacher, A. L. Nelson, C. V. Newsom, George Piranian, G. Y. Rainich, E. D. 
Rainville, M. O. Reade, C. C. Richtmeyer, E. H. Rothe, L. J. Rouse, T. R. 
Running, Hans Samelson, W. F. Smith, T. H. Southard, B. M. Stewart, C. G. 
Stripe, I. L. Stright, D. R. Sudborough, R. M. Thrall, C. P. Wells, and R. L. 
Wilder. 

The following officers were elected for the coming year: Chairman, Professor 
G. C. Bartoo, Western Michigan College of Education; Secretary-Treasurer, 
L. J. Rouse, University of Michigan. 

At the morning and afternoon sessions the following program of eight 
papers was presented: 

1. The geometrical representation of group C2, by Dr. Mary H. Payne, Uni- 
versity of Detroit, introduced by Professor L. S. Johnston. 

It is well known that the unitary matrices form a double-valued representa- 
tion of rotations in three-space, and that the entire group C2 is a double-valued 
representation of the true Lorentz group, consisting of rotations in three-space 
and ordinary Lorentz transformations. The functions on which the matrices of 
the group C, act are two component vectors of the forms 


0 0 
cos eilety) 3, sin ety) /2 


each multiplied by a constant. Dr. Payne showed that these two components 
may be interpreted as the relativistic distances of some point on the light cone 
from each of two reference points. 

2. “Kasner” monogenic functions, by Professor Emeritus V. C. Poor, Uni- 
versity of Michigan, introduced by the Secretary. 

This is a restricted class of polygenic functions dependent upon the existence 
of the Cioranescue derivative, a two-dimensional directed space derivative. 
This note is an analogue to the R. J. Haskell paper on areolar monogenic func- 
tions in the Bulletin of the American Mathematical Society, vol. 52, 1946, pp. 
332-337. 

3. Music and ternary continued fractions, by Professor J. M. Barbour, Michi- 
gan State College. 

The problem of dividing the octave (2:1) into a larger number of equal 
parts in terms of which the major third (5:4) and the perfect fifth (3:2) might 
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be expressed with nearly integral coefficients may be attacked by approximating 
the ratios of log 2:log 3/2:log 5/4 simultaneously by ternary continued frac- 
tions. The common division into twelve semitones gives the ratios 12:7:4. Some 
other known approximations are 19:11:6, 31:18:10, and 53:31:17. Slow con- 
vergence of the ternary continued fractions is desired in order to obtain ratios 
involving relatively small numbers. Hence in dividing to obtain the partial 
denominators for the left and right hand fractions it is well to divide by the 
larger remainder. For any right term S,=gS,-1+5S,-2, whereas for the left 
term Sp =QSn-1+Sn—-z-3 where k is the number of intervening right terms since 
the last left term. 

4. Transformations of sequences into regions, by Dr. George Piranian, Uni- 
versity of Michigan. 

The core of a sequence S, =x,+1y, is a certain non-empty set in the closed 
complex plane; in the case of a bounded sequence it reduces to the least convex 
set containing all limit points of the sequence (cf. Knopp, Zur Theorie der 
Lemilierungsverfahren, Math. Zeitschr., vol. 31, 1929, pp. 95-127). In the formal 
transformation ¢, = yor odnnse of an arbirtary sequence s, by an arbitrary 
matrix @,,4, some of the numbers ms — fail to exist; but the core c, of the 
SEQUENCE Wnp = GnpSp (P=0, 1, - - ) always exists. The sequence of cores 
is used to define a convex region I’ to be associated with the formal transforma- 
tion of the sequence s,. If the transformed sequence exists in the classical sense, 
the region T' is the core of the transformed sequence. 

5. Multinomial theorem, by Professor G. P. Loweke, Wayne University. 

By proper arrangement of terms within the summation, the expansion of a 
multinomial of m terms may be expressed as 


B=[(n—y—b— ++ +—y)/2] y=[n/3] 5=—[n/4] 


(A+ B+C+---+M)"= 


b=e 
c=[n/m—1]  p=[n/ m] n! 


M# 


where the brackets in the upper limit indicate the integral part of the frac- 
tion, and the number of the summation signs will be m when m <n and m when 
n<m. Expressed in this form the binomial theorem becomes, 


B=n/2 


= | 


This summation for the multinomial expansion can be used to derive the 
terms of the expansion directly, or if an arbitrary selection of terms is employed, 
the evaluation of the limits only will indicate the number of characteristic terms 
which arise in the solution. 

6. Undergraduate program in mathematics, by Professor C. V. Newsom, 
Oberlin College. 
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7. A machine shop problem, by Professor B. M. Stewart, Michigan State 
College. 

A geometric solution is provided for the angular adjustments, A, B, of a drill 
vise which is to hold a rectangular parallelepiped so that a vertical drill will bore 
a hole whose projections on a pair of intersecting faces of the parallelepiped will 
make specified angles, a, b, with respect to the line of the intersection of these 
faces. The solution is found as follows: A =90°—a; B=arc tan (cos a tan db). A 
nomogram is provided for the latter relation. 

8. On graphical methods for elementary differential equations, by Professor 
Wilfred Kaplan, University of Michigan 

The purpose of this paper is to indicate how interest can be added to the 
course in elementary differential equations by a closer correlation between 
analytical and graphical methods. As illustrations, some significant properties 
of the graphs of linear, homogeneous, and exact equations, are pointed out. 
L. J. Rouse, Secretary 


CALENDAR OF FUTURE MEETINGS 
Thirty-first Annual Meeting, Athens, Georgia, January 1, 1948. 


The following is a list of the Sections of the Association with dates of future 
meetings so far as they have been reported to the Secretary. 


ALLEGHENY MOUNTAIN 

ILLINOIS OKLAHOMA 

InpIANA, Muncie, Oct. 17, 1947 Paciric Nortuwest, Eugene, Oregon, 
Iowa, Fairfield, April 1, 1948 March, 1948 

KANSAS 


PHILADELPHIA, Bryn Mawr, Pa., Novem- 


KENTUCKY ber 29, 1947 
LouIsIANA-MIssIssIPPI Rocxy MountTAIN 


MARYLAND-DIsTRICT OF COLUMBIA-VIR- 


SOUTHEASTERN 
GINIA 
MetTRopouitan New Yorx SOUTHERN CALIFORNIA, Redlands, March 
MICHIGAN 13, 1948 
MINNESOTA SOUTHWESTERN 
MIssourI TEXAS 
NEBRASKA Uprer NEw York STATE 
NORTHERN CALIFORNIA, Berkeley, Janu- Wisconsin, Beloit, May 8, 1948 


ary 24, 1948 


] 
] 
| 
q 
¥ 
¥ 
: 
é 
& 


FUNDAMENTALS OF STATISTICS 


J. B. SCARBOROUGH, U. S. Naval Academy 
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Here is a clear, concise text covering the fun- 
damental principles and techniques of statistical 
analysis for the lower college level. Notable 
features are the early introduction and consistent 
use of the computing scale, the restricting of the 
terminology of moments and of the “normal 
curve” to conform to student progress and thor- 
ough discussion of the function of probability in 
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TEXTBOOK NEWS 


SOLID ANALYTIC GEOMETRY 


By 
John M. H. Olmsted, Ph.D. 


University of Minnesota 


a 


HERE is a new textbook designed for courses on the junior-senior or early graduate level in 

analytic geometry of three dimensions. The book contains the material for a very rich 
and extensive course, so arranged as to be conveniently adapted to courses of various content, 
and presented in a lucid understandable manner. Other features of the book are the simplicity 
of the presentation, the emphasis on logical reasoning and method, the many illustrative 
examples, over one thousand exercises, the systematic elementary treatment of Matrix Algebra, 
and the clarity of the fifty line drawings, representing surfaces and other spatial objects. 8vo. 
257 pages. $4.00. 


D. APPLETON-CENTURY COMPANY 
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ANALYTIC GEOMETRY 
By David S. Nathan and Olaf Helmer 
@ Careful training in identifying and drawing loci from their equations and in 


setting up the equations of loci—a unique, simple and rapid treatment of conic 
sections—these are but two of many distinctive features in this superior text. 


PLANE TRIGONOMETRY 
Revised Edition 


By Fred W. Sparks and Paul K. Rees 


@ An amazingly simple, efficient text showing the principles of Plane Trigonometry 
and their application. During the past decade the authors worked daily to increase 


clarity for the student. Every page of this revision has been proved successful in 
the classroom. 


INTRODUCTION TO COLLEGE MATHEMATICS 


By Carroll V. Newsom 


@ Here is the ideal mathematics study—a special two-semester terminal course— 
for that majority of college students who do not major in mathematics. The result 


of 12 years’ research, and pre-tested at two colleges, this text is the first ever 
written for this important group. 


MATHEMATICS IN HUMAN AFFAIRS 


By Franklin W. Kokomoor 


@ An unusual and absorbing text combining the history and the study of mathe- 
matics from algebra to calculus. Presented as exciting narrative, with human-interest 
applications and explanations simple enough for students with little previous training. 


SEND FOR YOUR APPROVAL COPIES 


a PRENTICE-HALL, INC. 70 FIFTH AVENUE, NEW YORK 11 
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Just published! 
COLLEGE ALGEBRA 


By Moses Richardson, author of FUNDAMENTALS OF MATHEMATICS, 
Former Member of Institute for Advanced Study at 
Princeton 


This new, remarkably lucid exposition of College Algebra 

® gives unusual insight into sound mathematics without excessive rigor. 

© is far more complete than most texts. 
Beginning with a full review of elementary algebra on a mature 
level, it covers thoroughly all conventional subjects, and adds op- 
tional material now needed by science students and students pre- 
paring to specialize in mathematics. 


® stresses fundamental concepts—but not at expense of technical 
achievement. 


® motivates and justifies techniques at the student's level. 
© includes ample exercises, many with practical application. 
First printing of 20,000 copies! 


Coming this fall! 


CALCULUS AND ITS APPLICATIONS 


By Raymond Douglas and Samuel Demitry Zeldin, Massachusetts 
Institute of Technology 


This beginning course for students of mathematics, physics, chemistry 
and engineering is designed to cover a complete calculus program in 
an unusually short time. The text is particularly suited to engineering 
and technical schools and universities where the emphasis is on tech- 
nique. 


Note these significant qualities: 


© Presents material with maximum of simplicity and order. 
® Combines brevity with sufficient rigor. 

© Emphasizes essential topics, stresses applications. 

® Includes many topics not found in one book. 


In Prentice-Hall Mathematics Series, Dr. Albert A. Bennett, Editor 


SEND FOR YOUR APPROVAL COPIES 


EY) PRENTICE-HALL, INC. 70 FIFTH AVENUE, NEW YORK 11 
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Intermediate Algebra 
By BRITTON and SNIVELY 352 pp., $2.00 
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Plane Trigonometry 
By WILLIAM K. MORRILL 245 pp., $2.50 


An important revision of a successful text designed to provide the student with a 
complete study of the fundamentals of trigonometry. Among many outstanding 
features, it offers a study of the fundamentals of trigonometry, a chapter on the 
problems of elementary physics, engineering and navigation, and an enlarged sec- 
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Mathematics of Finance 
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Differential and Integral Calculus 

By Harotp M. Bacon, Stanford University. 771 pages, $4.00 
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fundamentals of the calculus available. Presents the material so that the student can 


master it with a minimum of explanation by the instructor. Illustrative examples 
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Elementary Differential Equations. New third edition 
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A thorough revision of a standard textbook especially suitable for students of 
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optional chapters and supplementary section to facilitate using 

the essential course content; (2) fresh exercises except for the 
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alterations suggested by classrom use of the preceding edition; 

and (4) pleasing new format, clear type, wider pages. 362 pages | 
of text. 
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A rigorous course in plane analytic geometry with chapters on 
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By E. RICHARD HEINEMAN 
Associate Professor of Mathematics, Texas Technological College 
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calculus. The explanations of factoring and of operations with zero are 
unusually clear and explicit, as is the explanation of the troublesome question 
of the “sign of the fraction.” The book contains extensive sets of carefully 
graded problems. To be published in the fall. $3.00 (probable) 


The Mathematical Theory 
of Finance 


By KENNETH P. WILLIAMS 


Professor of Mathematics, Indiana University 


This text places emphasis on clear thinking, on thorough analysis of problems 
to find the simplest method of computation, and on thinking out solutions 
rather than searching for special formulae. Algebraic problems as well as 
numerical problems are given, and investment functions are included in the 
tables of logarithms. Revised Edition to be published in the summer. $4.50 
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Elements of 
Symbolic Logic 


By HANS REICHENBACH 
Professor of Philosophy, University of California at Los Angeles 
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this text then takes up the calculus of propositions, for which the truth table 
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of the higher calculus of functions, conversational language, and connective 
operations and modalities. $5.00 
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THE HERBERT ELLSWORTH SLAUGHT MEMORIAL PAPERS 


The editorial committee is happy to present this, the first of the Slaught 
Memorial Papers. In view of the very great interest which Professor R. E. Langer 
has had in these Papers from their very inception, it is particularly fitting that 
he should be the author of the first of the series. 

At this time it may be appropriate to give a brief account of the origin of 
these Papers, and an indication of their purpose. This can perhaps best be done 
by quoting from the reports of two different committees of the Association. In 
the Montuy for February, 1940, there appears the report of a “committee to 
review the activities of the Mathematical Association of America.” Along with 
other recommendations, this report, which was prepared by Professor Langer 
as chairman, suggested that the Association establish a series of expository 
pamphlets which might well take the form of a memorial to Professor Slaught. 
Subsequently, another committee was appointed, with Professor C. V. Newsom 
as chairman, to make specific recommendations for the establishment of such 
a series. The report of this committee appeared in the MONTHLY for February, 
1941. 

From the report of the “Langer Committee” we may quote as follows: 

“The encouragement and sponsorship of expository and critical writing is 
one of the objectives of the Association which enjoys the unanimous support of 
the members. There is a ready welcome and a general demand for more readable 
scholarly papers on all kinds of mathematical subjects from the classical to the 
modern, from the elementary to the advanced, on theory, on applications, on 
history, or on philosophy. In the past there have, of course, been the Carus 
monographs, and from time to time excellent papers in the MONTHLY. There 
seems, however, to be at the present little or no means for the ready publication 
of writings which in length are intermediate betwen the relatively few pages of 
a journal paper, and the relatively many pages of a complete monograph. Such 
papers, say in length between twenty and a hundred pages, could be profitably 
written on subjects in many categories, including among others, elementary 
introductory expositions of theories and their applications, more advanced 
expositions and interpretations of modern viewpoints and theories, philosophical 
essays and criticisms, broad historical accounts of important schools, or bio- 
graphical accounts of individuals.” 

The report of the “Newsom Committee” said in part: 

“The principal conclusion reached by this investigation is that there is a 
widespread interest in additional expository writing of the type discussed in the 
report of the ‘Langer Committee,’ and that those who are now sponsoring series 
of expository monographs would welcome the creation of additional opportuni- 
ties for the publication of studies pertaining to mathematical subjects. In truth, 
the members of the committee have been impressed with the enthusiasm which 
has been displayed by those who have given opinions relative to the possibility 
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of a new publication program sponsored by the Association. Syntheses of mod- 
ern investigations in many fields of mathematics seem to be wanted by college 
men who do not have an opportunity to follow developments in the mathe- 
matical literature. Instructors in our junior colleges and secondary schools who 
may have a limited preparation in mathematics are seeking easily accessible 
accounts of some of the older theories. Some correspondents have expressed the 
belief that there is an amazing dearth of readable mathematical material for 
college students who have studied little beyond the calculus. And finally, some 
have emphasized that the interest of the American public in mathematical at- 
tainments and methods needs to be cultivated; this interest is attested to by 
the recent wide sale of a few popular books upon mathematics.” 

In accordance with the recommendations of the “Newsom Committee,” the 
Board of Governors authorized a series of expository pamphlets to be known 
as the “Herbert Ellsworth Slaught Memorial Papers.” The long delay in the 
actual appearance of the first of these Papers was largely caused by the demands 
of the war which left little or no time for the writing of mathematical exposi- 
tions. 

The Slaught Memorial Papers are to be published in the form of supple- 
ments to the MONTHLY and, at least for the present, are being sent free to all 
subscribers. The success of this project will depend on the interest of mathe- 
maticians generally and, more particularly, upon the co-operation of com- 
petent scholars who will be willing to devote sufficient effort to the difficult but 
worth-while task of writing elementary expositions of their respective fields of 
interest. 

The editorial committee through the undersigned will welcome suggestions 
from any interested persons and, in particular, will be glad to hear from pro- 
spective authors of expository articles which might be suitable for publication in 
this series. 


N. H. McCoy 
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FOURIER’S SERIES 
THE GENESIS AND EVOLUTION OF A THEORY 
R. E. LANGER 


PREFACE 


In choosing to present in this exposition some chapters from the theory of 
the representation of arbitrary functions in infinite series, I have done so in the 
belief that this subject has an unusually broad appeal. For in singular measure 
it serves both theoretical and practical ends. The pure analyst finds in it a wealth 
of structure and subtle inter-relationship, while the applied mathematician and 
the related scientist find in it, no less, a tool of almost endless flexibility and use. 

The simpler formal elements of the theory of trigonometrical infinite series 
are, it may be assumed, in some measure familiar to all who aspire to a level of 
mathematical attainment above the elementary one. Presentations of them in 
text-book form are common, and many of them eminently readable. It is not 
my purpose to duplicate any such expositions of fact and procedure, but rather 
to present here other matters less usually considered. These I have, in the main, 
centered about two focal theses, namely first, a sort of case history of the incep- 
tion of the theory and its development to the stage attained by Fourier, and 
second, a generalization of the theory in which the trigonometrical form of it is 
subordinated to the status of a mere special case. 

In its modern form the theory of Fourier’s series and its applications to prob- 
lems of physics admit of presentation in a direct and logical manner that is, on 
the whole, strikingly economical in design. The reasoning is straight-forward and 
to the point, and has at almost every turn an aspect of complete inevitability. 
The trigonometric formulas invariably appear to fit the needs at issue with such 
precision and neatness as could not have been more so had they been specifically 
tailored to the purpose. So completely is this true, that it seems no far cry to 
the suggestion that the whole structure might be the creation of some single 
master architect, who, in his genius, could draw to hand the exact and unerring 
means for an orderly and consummate unfolding of all the whole essential 
machinery of thought and analysis. Of erstwhile possible deficiencies, no trace 
is left revealed. 

Well developed mathematical theories are prone to seem like that, and in the 
deceptiveness of this there is weakness as well as strength. The craftsman, whose 
concern with the theory is motivated by the mere search for a tool, naturally has 
small interest in the cruder forms of it that are now obsolete. For the student 
whose concern is with ideas no less than with facts, on the other hand, the too 
finished result is often concealing as well as revealing. The confusion of germinal 
ideas, the labor and stumbling of the early advance, the frustrations in imprecise 
notions—all these are matters which for a speedy mastery of the facts are well 
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left aside. But precisely in these as in none other are to be discerned the creative 
imagination, the initial inductions and the logical strategy by which the final re- 
sult was shaped. When, as in the present instance, the ingenuity and the tech- 
nical exploits which gave the impetus and direction were those of such masters 
as the Bernoullis, Euler, d’Alembert, Lagrange and Fourier, it need hardly be 
feared that a review of the course of the developments will prove to be an unre- 
warding venture. In the growth of mathematics, much more than in its refined 
and polished results, is its living character evidenced. 

Generalization is the medium through which the mathematician constantly 
seeks the enlargement of his conceptions and understanding. The vista revealed 
by an extant theory may be broad, but it is broadened further by generalization. 
And from the more expansive viewpoint the scene may be revealed not only 
more amply but also more distinctly. A greater simplicity in the intrinsic plan 
may be discernible, for many features originally judged to be quite essential may 
be shown, on the contrary, to be in fact merely incidental or fortuitous. This is 
quite the case with the Fourier theory. Its dependence upon the trigonometric 
formulas of combination is so conspicuous as to seem to be the very essence of 
it. And if in the related theories of representations in series of Bessel functions or 
Legendre polynomials etc. other combination formulas are basic, these in turn 
generally seem, if anything, even more specialized. A true generalization, from 
which the Fourier theory may be drawn forth as a special case, is the theory of 
ordinary differential boundary problems in which the fundamental interval of 
the variable is one upon which the differential equation is without singular 
points. 

The discussion which I have given here is intended to serve these two pur- 
poses. In the first part the theme is historical. It centers about the incipience 
and the classical development of the theory, and is in fact a digest of some works 
by different masters through which conspicuous advances were made. This 
comes to its terminus with the discussion of Fourier’s deductions, and therewith 
the historical thread is definitively dropped. In the second part, which is de- 
voted to the generalization, the purpose is purely expository. The material is 
there set forth in as elementary a manner as I found possible, not with the gen- 
erality in which it exists in the literature, but with such generality as seemed to 
be adequate to the display of its essential character. The fact that this larger 
theory embraces that of Fourier, and the manner in which it does so, is shown at 
appropriate points by drawing the trigonometric formulas forth as specializa- 
tions obtainable without any peculiar implementations from the more general 
relationships derived. 

I believe that in the main the paper will be readable for students who in 
mathematics have gone but little beyond a good course in the calculus. The 
simplest facts about infinite series and differential equations, the formulas for 
the trigonometric functions in terms of exponentials, and such, have been as- 
sumed. Beyond that all pertinent deductions have been included until the clos- 
ing chapters are reached. Incidental material has, in part, been relegated to 
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appendices. In the final chapters the elementary theory of functions of a com- 
plex variable, and in particular the theory of residues is a requisite. 

In taking material from the literature, especially in the earlier parts of the 
paper, I have felt under no obligation to hold to the letter of the originals. The 
excerpts are, therefore, distinctly not to be regarded as facsimiles or verbatim 
reports. Although it was my intention to preserve the spirit, many formal 
changes were made, in part to bring the contributions from diverse sources under 
a consistent scheme of notation, but also in part to eliminate discursive material 
and to avail myself of such advantages as the presentation to modern readers 
might afford. The sources in the literature from which excerpts were made, or at 
which more extensive deductions may be found, have been indicated in the text, 
and are listed at the end of the paper. I have, however, made no attempt whatso- 
ever to be complete in this matter. 
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PART I 


CHAPTER 1 


Introduction. It is perhaps true that few mathematical doctrines are built 
around facts which on first acquaintance seem so surprising as those of the infi- 
nite trigonometric series. That some such series represent functions is obvious. 
Simple examples are easily constructible in any abundance. How broad or how 
deep the adaptability of these series for functional representation may be, is, on 
the other hand, far from easy to see. The terms of the series are but elementary 
functions of a very simple type—the sines and cosines of the multiples of an 
angle. The theory of power series does not promise much by analogy, since those 
series can represent only such functions as have all the regularity which unlim- 
ited differentiability assures. It might well be expected, therefore, that the prop- 
erty of trigonometric representability attaches only to the functions of a quite 
restricted class. Historically that was the opinion which originally held sway, 
and which was very generally maintained. It was, in fact, so long maintained 
and so tenaciously, even by the greatest of mathematical masters, and in the 
face of most insistent evidence to the contrary, that the final breach with it took 
on quite definitely the character of an emancipation. ; 

The concept of the function lies at the very heart of mathematical analysis. 
As it is now currently accepted it is a notion of very great breadth, covering very 
general interdependencies of variables upon each other. During the eighteenth 
century this concept was not only much more restricted, but precisely what its 
content and delimitations were had not yet been brought to any complete or 
clear formulation. Imprecise notions rarely fail to breed confusion, and in this 
respect the functional notion of that time was in no way an exception. It was 
differently conceived by different investigators. And these latter then disagreed 
among each other because unconsciously they talked at cross purposes. The 
written words flowing from different pens had different meanings. While, for 
instance, the function and the analytic formula were one to d’Alembert, the func- 
tion was thought of as a graph by Euler, and probably meant something else 
again to still another. 

The basis of the functional notion was originally drawn, of course, from ob- 
servations upon concrete examples—in the main from such functions as we 
now designate as of elementary type, such as present themselves in the simpler 
applications of mathematics to the problems of physics. Such functions are al- 
most invariably expressible by formulas. They generally have comparatively 
simple, orderly, and continuous graphs, and the identity of any two of them is 
restricted to isolated values of the variable. Inasmuch as this category includes 
no examples of distinct functions whose graphs have an entire arc in common, 
it was no more than natural then to consider that generally the course of a func- 
tion over any interval was determinative and completely identifying, so that the 
graph over its entire range of definition was to be thought of as unambiguously 
fixed. Functional relationships such as are now commonly dealt with, in which 
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the variables are related by different laws in different parts of an interval, were 
not thought of then as subsumed in any single function at all, but were regarded 
rather as a composite of a plurality of functional fragments. The possibility of 
representing such a conglomerate by a single formula was not even conceived of. 

The eighteenth century stands out in mathematical history as an era of great 
genius. Through the work of an astonishing array of masters the science was ex- 
tended and broadened by the opening of many new fields. Technical skill at- 
tained to extraordinarily high levels and new ideas were crowded one upon the 
other. And yet through this period the facts of the trigonometric series withheld 
themselves. Euler, d’Alembert, Lagrange and others walked upon the very edge 
of them without falling upon them. A more conspicuous example of the confining 
effects of preconceptions is hardly to be found. The break with all this remained 
to become the accomplishment of the next century, the personal achievement 
of Fourier. Once the step to a broader conception of the function had been made, 
the results of computations upon trigonometric series could be given a much 
more inclusive interpretation. As is now generally familiar, such series may rep- 
resent functions which are not only discontinuous, but which may be quite 
arbitrary in the sense that over different portions of the interval they may ac- 
cord with laws that need have no logical relation with each other. Computations 
upon even a few of the initial terms of the series often reveal these facts quite 
clearly. 

As Fourier announced his famous theorem it was to the following effect: 

Any single-valued function f(x) defined over an interval —1<x<l, is repre- 
sentable over this interval by a series of sines and cosines in the manner 


k k 
(1.1) f(x) = + cos be sin =|. 


kel 
In this representaion the coefficients are those which are computable from the func- 
tion f(x) by the formulas 


1 kas 
f(s) cos—ds, 
l 


(1.2) 


kas 
b =—] f(s) sin —ds. 
l 


If the interval over which the representation is to maintain is only 0<x<l, then 
either sines or cosines alone suffice, the series being in the one case 


k 
(1.3) 
k=l 
with the coefficients 
kas 
(1.4) b= —f f(s) sin — ds, 
Jo 
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and in the other case 


(1.5) fla) => + 


k=1 


with the coefficients 
a krs 
(1.6) a, = —f f(s) 
0 


This theorem, in this utter generality, is now known, of course, to be not 
strictly true, even with the modern notion of the term function. The restrictions 
which must be applied are, however, in no way gross. They are, on the contrary, 
so subtle that in the pursuit of them much further clarification of analytic 
notions was achieved. The réle which the trigonometric series have played in the 
development of precise conceptions is, on these many accounts, one of unusual 
interest to the student of the evolution of mathematical ideas. [1] 


CHAPTER 2 


Of mathematical applications to physics. By the time of the year 1725, a 
decade after the death of Leibniz and near the close of Newton’s life, familiarity 
with the formal processes of the calculus had become widely disseminated, and 
facility in the use of these new techniques had been developed to a very sub- 
stantial degree. In particular their effectiveness as instruments for the treatment 
of problems in mechanics had been generally recognized. This science was, 
therefore, under assiduous study. The more immediate of its problems—in the 
main those centering themselves around the motions of single mass particles— 
had already been pretty effectively brought to their solutions. The forefront of 
interest had, therefore, already been pushed beyond them to matters of greater 
complexity, such as presented themselves in connection with the motions of 
bodies with several or many degrees of freedom, or even with the reactions of 
flexible continuous mass distributions. Problems in the vibrations of elastic 
bodies in particular had begun to receive attention, and in the ensuing period 
these were to preoccupy an increasing number of investigators. The scientific 
literature of the following half century is, therefore, heavily interspersed with 
memoirs bearing upon this field. One may readily conjecture that the prob- 
lems to be found there must have seemed almost endless in their abundance 
and variety. They must, moreover, have exerted a very strong fascination, 
if for no other reason, then because of the evident suggestion that in them lay 
an important key to an analytic mastery over the manifestations of nature. 

Any review of the activities of the time show this to have been so. Important 
and difficult works were produced in great variety—among them, to name a few, 
investigations upon the oscillations of plates in vacuum or immersed in fluids, 
upon rods suspended from fixed or flexible mountings, upon jointed pendulums, 
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upon heavy dangling chains efc. etc. But among all these the researches upon the 
motions of tautly stretched elastic strings or wires were to be especially sig- 
nificant. Though these were perhaps of no greater mechanical importance than 
the discussion of many another problem, they did contemporaneously assume a 
quite disproportionate prominence, and they still do so in retrospect. They be- 
came, namely, a conspicuous point of impact for many divergent conceptions 
and opinions, and in the réle thereby thrust upon them they became crucial 
to the development not only of mechanics but of mathematics over a much 
wider range. Inasmuch as the ideas at issue are in large measure our essential 
concern in this discussion, we shall have to give a considerable modicum of our 
attention to this problem of the string in subsequent pages. 

From the very start of investigations upon them it seems to have been as- 
sumed that continuous material bodies could, for the purposes of analysis, be 
approximated by systems of discrete mass particles. The conception of an ex- 
tensive body as composed of particles is a very natural one. Quite apart from 
whether this viewpoint was intended in the end to be philosophically maintained 
or not, it was apparently seen to suggest more strongly than any other a prac- 
ticable mode of procedure. It requires no vivid imagination to picture the dis- 
crete system as merging into the continuous one as the size of the individual 
particles is diminished indefinitely while their number and density is correspond- 
ingly increased. If in the analysis of a finite approximating system the formulas 
can be so framed and dealt with that results are deducible from them without 
any actual specification of the number of the particles that are involved, then 
this number, retaining its generality, may be assumed to figure in the results in 
the way of a free parameter. This parameter may then logically be made the 
crux of limiting considerations, and through this means the physical transition 
from the discrete to the continuous configuration may be thought of as imple- 
mented mathematically by a passage from the finite to the infinite. 

It may well be recognized that success in the execution of any such subtle 
program as this would be contingent both upon superior insight and a high 
level of technical skill on the part of the investigator. Such is the case, no less, 
with almost any application of mathematics to a phenomenon of nature. The 
actual responses of physical bodies to ponderable influences are invariably of a 
discouraging degree of complexity, and this is generally due much less to the 
influences that are primarily under scrutiny than to the many others that are 
inevitable and yet really incidental and effectively irrelevant. Were complete 
recognition to be given to all these latter, the formulation of a natural problem 
would, almost without exception, be quite submerged in intricacies of detail. On 
the physical side the fortuitous distracting features might well obscure the 
salient ones, and mathematically they might well throw the problem far beyond 
the range of possible solution. 

At the very outset, therefore, it is usual and necessary to regard the physical 
configuration not as it actually is, but as it might be were it to be disencumbered 
of all but its primarily intrinsic features. The result of this is at once a simplifi- 
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cation and an idealization. It is this which is made the subject of analysis. The 
idealization is, of course, a departure from the physical reality, and this fact is 
certainly of no secondary importance. The original problem having been re- 
placed by another, no immediate guarantee exists, of course, that results de- 
rived for the latter have a sufficient relevancy or applicability to the former. It 
is evident that with too great a divergence from the original all practical purpose 
would thus be defeated. 

This is the consideration which checks the extent to which idealizing ab- 
stractions may be made. In determining upon these a fine sense of values and a 
depth of understanding are the indispensable guides. Whether in any case the 
permissible bounds were exceeded or not, must in the end ordinarily be decided 
after the fact by experiment. The applicability of the mathematically deduced 
theory stands or falls according as at strategic points its results agree or disagree 
with the data of observations. The decision, either way, implies no reflection 
whatever upon the theory’s logical soundness. Its implications bear only upon 
the legitimacy of the simplifications which were made in the determination of 
its basis, namely upon the insight with which suitability and adequacy in the 
initial idealizing approximations were sensed. 


CHAPTER 3 


The loaded string. A natural phenomenon that is universally familiar and 
frequently observable in our surroundings, and to which, despite this, the atten- 
tion is often sharply drawn, is that of the behavior of a stretched elastic string or 
wire in its response to a displacement from its state of equilibrium. In many cases 
this response is acoustically conspicuous, ranging from the hum of the heavy 
structural wire in the wind to the eloquent notes of the strings of a musical in- 
strument. And then again, not rarely, the response is visually noticeable, being 
sometimes marked by so curious a feature as the presence of nodal points that 
maintain the state of rest while the string between them is in violent agitation. 

The elastic string thus almost obtrudes itself upon the notice of the experi- 
menter, and, having drawn his attention, recommends itself in many ways. Its 
geometrical configuration is of the simplest sort, permitting the identification of 
any of its points by a single dimension. Its motions are, under many circum- 
stances, markedly regular, and respond promptly and prominently to adjust 
themselves to any quantitative modifications in the length, the tension, the 
weight or the initial state. It can hardly be looked upon as surprising, in virtue 
of all this, that the string should have been drawn under analysis at as early a 
time in the development of mechanics as that science became capable of dealing 
with continuous flexible bodies. 

The finite discrete system of particles that most naturally approximates the 
continuous material string is suggested by the string of beads. More precisely, 
it is to be conceived of as comprised of any number 2 of equal concentrated mass 
particles, that are mounted respectively at equally spaced points along a string 
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which is itself weightless, though strong, perfectly flexible, and elastic. Such 
“loaded strings” were used as the conceptual bases of a variety of significant 
investigations. Euler and Daniel Bernoulli based upon them, in 1732 or 1733, 
studies of the motions of heavy dangling chains, and Euler, among other things, 
regarding the particles as oscillating longitudinally, built upon them in 1746 
a theory of sound. [2] Mechanically simpler and chronologically even earlier than 
these researches, were certain investigations of John Bernoulli upon the trans- 
verse vibrations of a string with its end points fixed. We find in a consideration 
of this problem a convenient point of departure for our present discussion. 

Consider, therefore, a loaded string such as has been described, with loading 
particles having a total mass M. Let this string be thought of as held taut under 
a tension 7, the magnitude of this being so large that the ratio M/T is negligible. 
This last stipulation amounts, of course, to an idealization. Substantially its 
purpose is to discount the effects of gravity, and thus to concentrate the con- 
siderations upon those forces which spring from the tension alone. The diver- 
gence from physical actuality which this simplifying hypothesis sanctions is in 
many important instances of a very small amount. In the case of musical 
strings, for example, the ratio M/T which is to be ignored is quite commonly of 
a magnitude no greater than one one-thousandth. 

The initial state of the string in question is to be one of displacement from its 
equilibrium position, the forces which hold it in this state being coplanar and 
directed perpendicularly to the line through the string’s fixed end points. The 
particles of the string thus lie along some plane curve. We shall choose the plane 
of this curve as the (x, y) plane, and shall take the x-axis through the string’s end 
points with the origin at one of them. Under the hypotheses made, and in this 
system of reference, the x-axis then marks the string’s equilibrium position. 
With / designating the length of the string, with x»=0 its initial end, and with 


%1, X2, X, the equilibrium abscissas of the particles, the formulas for these 
are to be 
l 
k=0,1,2,--:,m. 
n 


It will be noted that in this assignment one of the particles is allocated to be 
mounted at the point x, =/ which is the terminal end of the string. The motiva- 
tion for adopting this arrangement, which may well seem a bit curious, is not 
by any means profound. It is merely one having some formal advantages, since, 
as can be shown, it leads to somewhat simpler formulas. The réle of the nth 
particle under these circumstances is; of course, an entirely passive one since the 
particle is constrained from all motion. 

In the accompanying figure three adjacent particles and the tensions operat- 
ing between them are schematically indicated. The position shown is one of 
displacement from the equilibrium, and the particle which in the latter state is 
located at the point (xz, 0) is here shown to have the coordinates (xj, y,). Simple 
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mechanical considerations based upon this figure readily lead to usable results, 
as we now propose to show. 


eet 


kel Xk 


The segment of the string which joins the (k—1)th and the kth particles is, 
as is clear from the figure, of the length (x¢ —x,_i) sec ax. In the state of equi- 
librium this segment would be of the length //n, and since, by the laws of elas- 
ticity, the respective lengths are to each other as the tensions, it must be con- 
cluded that 


l 
— sec ax:— = T;:T, 
n 
or 
(3.1) (x; — = IT; cos ax. 


Now the forces by which the particles are initially held displaced are by 
hypothesis in the direction of the y-axis. The tensions in the string segments 
which balance them must, therefore, have x-components which annul each 
other in pairs, namely they must be such that 


(3.2) Ty CoS = COS &=1,2,---, — 1). 


Since the right-hand members of the equations (3.1) are thus all equal, this 
equality must extend to the left-hand members as well. The differences 
(x¢ —x,-1) therefore have a common value, and since their sum is the length 
of the string this value is evidently //n. It follows that for each k the relation 
xe =x, maintains, namely that each particle has when displaced the same 
abscissa as in equilibrium. The measures of the displacements are thus simply 
the ordinates y;, and the relations (3.1) reduce to the forms 


(3.3) cos a, = T, k=1,2,--+,(n—1). 


From the initial position which has been described the string is now to be 
thought of as released, while in the state of rest, at an instant which is to be taken 
as the origin of time, =0. The motion into which the kth particle springs is, 
of course, then governed by Newton’s law 
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the force F; acting upon it being shown by the figure to have the value 
Fy = Tesi Sin — Ty sin ay. 
By virtue of the relations (3.3) this formula is alternatively expressible as 
F, = T(tan — tan ax), 


and in terms of the codrdinates this is 


(1/n) 
Hence if the constant a? is defined by the relation 
(3.4) a= 
M 
the equations of motion are 


In many respects the simplest modes of vibration which the string is capable 
of are those in which the ordinates y, maintain constant ratios to each other. 
In these so-called normal vibrations all particles of the string traverse their 
positions of equilibrium in synchronism, and their displacements are expressible 
as functions of the time by formulas of the type 


(3.6) = uxp(t), k= 0, 


in which the coefficients u, are constants, with u)=0, vu, =0, and ¢(¢) is common 
to them all. The fact that the motion originates from the state of rest is then ex- 
pressed by the relation ¢’(0) =0. 

Now the substitution of the forms (3.6) into the equations (3.5) gives to 
these latter the aspect 


2 


From this it is clear that the second derivative $’’(#) stands in a constant ratio 
to the function ¢(#) itself, namely that 


(3.7) = 


=— co), = 0, 


the constant c being one for which the relations 


: 
; | 
= 
M dy 
k 
n dt? 
be 
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cl\? 
na 


together with u)=0, u,=0 maintain. In terms of the coefficient qg given by the 
formula 


(3.8) = (<)- 2, 


the values u, must thus be solutions of the algebraic system of equations 


u = 0, 
(3.9) + que + = 0, k= — 1), 
u, = 0. 


This system is neatly solvable (c.f. appendix I) having a non-trivial solution 
when and only when the coefficient g has one of the set of characteristic values 
Qi, G2, * * * » Qn—1, Given by the formulas (1.5). The values of ¢ which respectively 
correspond to these under the relation (3.8) are those of the set 


2na 
(3.10) = — sin y=1,2,---,(n—1), 
l 2n 
and the solution u,, of the system (3.9) which exists for the value c, is obtain- 
able from the formulas (1.6). It is 
kur 
Uy. = A, sin —» 
n 
with A, designating any constant. Since when c =c, the equations (3.7) are solved 
by the function 


= cos (ct), 


or by a constant multiple of this, it may be drawn from the relations (3.6) that 
the ordinates in any normal vibration of the loaded string must accord with the 
formulas 


kor 2ant vO 
(3.11) = A, sin — cos (= sin b= 
n l 2n 


A loaded string carrying m particles is thus seen to be capable of sustaining 
(n—1) distinct motions of the normal type, these being given by the formulas 
(3.11) in correspondence with the indices y=1, 2, ---, (n—1). 

Under assumptions that were somewhat more restrictive than those which 
we have here imposed, these normal vibrations were considered by John Ber- 
noulli as early as the year 1728 in the cases of loaded strings in which particles 
up to eight in number were involved. [3] 
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CHAPTER 4 


The equations of motion for the continuous string. When the differential 
equations (3.5) for the motions of the loaded string of m particles have once been 
deduced, two alternative modes of procedure for utilizing them toward the 
ultimate purpose of an analysis of the vibrations of a continuous string sug- 
gest themselves. On the one hand the equations may be integrated, as has 
already been done in the preceding chapter, and the resulting finite equations 
(3.11) may then be subjected to the limiting process in which the number 
n is indefinitely increased. On the other hand this limiting process may be ap- 
plied directly to the system of equations (3.5) itself, and the integration of the 
result may then subsequently be undertaken. Both of these procedures were 
carried out in the first half of the eighteenth century. As we shall see, their 
results are of quite dissimilar aspects. Indeed they seemed to the men of the 
time to be no less than contradictory, to the extent that the proponents of 
either method saw no alternative but to reject the other. That no real dilemma 
was actually involved therein at all, came to its realization only half a century 
or more later. The clarifications of ideas by which the way out of the quandary 
was ultimately found are of especial interest to us here. We propose, therefore, 
to pursue the analysis of the two mentioned procedures to such points, at least, 
as afford some surveys of their conclusions. 

Returning then, to begin with, to the equations (3.11), let any positive 
integer v be chosen. Once chosen, pv is to be regarded as fixed. A loaded string 
with particles in number greater than v, (n>v), may then be thought of, and 
for the normal motions of such a string the equations (3.11) are derivable. Let 
the attention then be fixed upon any one of the particles of this string, and let 
its abscissa and ordinate be designated by x and y(t, x). If this particle, in the 
enumeration that was adopted, is the kth one, the equalities 


(4.1) x= y(t, x) = (2), 


evidently maintain. The respective kth equation of the set (3.11) may then be 
written in the manner 


(4.2) y(t, x) = A, sin cos =| 


with the significance of 8, given by the formula 
sin 
2n 
(= 
2n 


(4.3) = 
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Suppose now that the parameters m and & are increased, and indefinitely so, 
in any way such that the ratio k/n remains fixed. The point x then clearly re- 
mains invariant, and inasmuch as the formula (4.3) familiarly shows that the 
value 6, approaches the limit 1 it follows that the relation (4.2) passes, as 
n— ©, into the limiting form 


x at 
(4.4) y(t, x) = A, sin cos 


This result is now evidently to be accepted as a formula applicable to the con- 
tinuous string and representing a normal vibration of it. Inasmuch as the integer 
v was initially open to an arbitrary choice, the inference that infinitely many 
such normal motions are possible and that they are given by the formulas (4.4) 
in conjunction with the indices y=1, 2, 3, - - - is inevitable. 

The simplest of the normal vibrations, namely that described by the formula 
(4.4) with »y=1 was deduced by Brook Taylor at as early a date as 1713. [4] In 
this motion the string ‘vibrates without nodes and emits its fundamental tone. 
The existence of other normal motions, namely those associated by the formula 
(4.4) with other values of v and in which the string emits its various over-tones, 
were known later to Daniel Bernoulli. We shall have occasion to return to this 
matter again. 

The alternative procedure, to which we now turn, is associated most prom- 
inently with the names of d’Alembert and Euler. With the notational changes 
(4.1) and with the definition of Ax by the formula Ax =1/n, the kth one of the 
equations (3.5) may evidently be written in the form 


E x + Ax) — 2y(t, 2) + y(t, 


ot? (Ax)? 


Now whenever the function y(¢, x) is one which is twice differentiable as to x, 
its second partial derivative with respect to x is obtainable as the limit of the 
difference quotient within brackets on the right of the equality (4.5) as Ax—0, 
namely as n—. Basing himself upon this observation d’Alembert deduced in 
1747 the partial differential equation 


d*y(t, A? y(t, 2) 


4.6 
ot? Ox? 


for the motion of the continuous string. This result is of course still a standard. 
A solution of it, if it is to represent the ordinates of a string that is fastened 
at its end points upon the x-axis and that springs into motion at the time ¢=0 
from the state of rest and from the position of a curve y=f(x), must, of course, 
also fulfill the conditions 


= 
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y(t, 0) = 0, 
y(t, 1) = 0, 
(4.7) dy(t, | =0, 
ot t=0 
x) = f(x). 


The problem as thus formulated was solved by d’Alembert through the following 
ingenious use of familiar formulas from the calculus. [5] 
In terms of the abbreviations 


1 Op 
ax a? at 
whereas it is familiar that generally 
9p 
a ax 


By the use of these relations the standard formula 
dq = —dt+— dz, 
Ot Ox 


may, however, be written thus 
ap 1 op 
dq = —dti+— —dx. 
q Ox a* dt 


From this together with the companion formula 


op op 
dp = — dt+— dx, 
ot Ox 
it may be seen at once that 


1 10 
a dt 


a@ Ox 


p —1 0p 1 dp 
— = |— —+— — ld(at — x). 
‘) | a a? 4 


Consider the first one of these equations. The quantity (p/a+q) is a function 
of ¢ and x. These variables are in turn determinable from the combinations 


(4.9) 


the differential equation (4.6) is expressible in the form Bee 
= 

d d(at + x), 

U 
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(at+x) and (at—x). It follows from this that it is permissible to regard the 
quantity as a function of the variables (at+x) and (at—x), and that accordingly 
its differential is given by the formula 


| 


But by the evaluation (4.9) this differential includes no term in d(at—x). The 
coefficient of this term must, therefore, be zero, namely 


O(at — x) 


This, however, is in effect the assertion that the quantity (p/a+q) does not de- 
pend upon the variable (at—x) but is a function of the remaining variable 
(at+x) alone. A similar chain of reasoning shows that the quantity (p/a—q) is 
a function of the variable (at—x) alone, namely that with appropriate functions 
designated by ¢ and y, 


— x). 
a 
If these relations are now multiplied respectively by the factors (adt+dx)/2 
and (adt—dx)/2 and are then added, the result on the left of the equality is 
pdt+q dx, which is identified as dy. Thus 


dy = 3¢(at + x)d(at + x) + 4y(at — x)d(at — x). 


In this formula each term is an exact differential. An integration is, therefore, 
possible and shows that 


(4.10) y(t, x) = + x) + (at — x), 


the functions ® and W being indefinite integrals of ¢ and y respectively. With 
the attainment of the result (4.10) d’Alembert had deduced the fact that every 
solution of the partial differential equation (4.6) is of necessity expressible as the 
sum of a function of the variable (a#+.x) and a function of the variable (at—<). 
It is a simple matter to show conversely, by direct substitution, that if @ and V 
are any suitably differentiable functions, the formula (4.10) does in fact give a 
solution of the differential equation. 

To apply to the particular case of the vibrating string the solution (4.10) 
must furthermore conform to the conditions (4.7). Of these the first one, which 
when applied to the relation (4.10) assumes the form 


( a ) 

0. 

4 
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+ 7¥(at) = 0, 


shows that the function V must be identical with — ®. That being so the second 
condition takes the form 


(4.11) + 1) = — J), 


and since this is to be an identity in ¢ it shows that the function ® must be pe- 
riodic with the span 2/ as a period. The third condition (4.7) reduces to the re- 
lation 


= x). 
Upon an integration this becomes 
(4.12) @(x) = — &(— x), 


and thus characterizes ® to be an odd function. The last condition, which must 
maintain over the string’s length, reduces then to the relation 


(4.13) = f(x), Osxsl. 
In total it is to be concluded, therefore, that for the vibrating string 
(4.14) y(t, x) = $b(at + x) — — x), 


every motion of the string being so representable with an appropriate function 
®. In the instance of any particular motion, in which the curve from which the 
string is released is y=f(x), the function ® that is concerned is determined over 
the interval (0, /) by the relation (4.13) and is defined for all other arguments by 
its character of being odd and periodic. 


CHAPTER 5 


The d’Alembert-Euler-Bernoulli controversy. [6] The method of d’Alembert 
in his analysis of the problem of the vibrating string was also the method chosen 
by Euler. Superficially, therefore, the initial memoirs of these masters, written, 
as they were at short intervals of each other, differed mainly in their details. In 
their over-all aspect they resembled each other markedly, at least insofar as their 
formal features were concerned. Only below the surface did the lines of thought 
show themselves to be divergent, as sharply so, at points, as were the two men 
in the characters of their genius. Euler’s temperament was an imaginative one. 
He looked for guidance in large measure to practical considerations and physical 
intuition, and combined with a phenomenal ingenuity an almost naive faith in 
the infallibility of mathematical formulas and the results of manipulations upon 
them. D’Alembert was a more critical mind, much less susceptible to conviction 
by formalisms. A personality of impeccable scientific integrity, he was never 
inclined to minimize short-comings that he recognized, be they in his own work 
or in that of others. 


: 
. 
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To Euler the solution of the problem of the string seemed definitive. Since 
any and every motion originating from the state of rest would necessarily stem 
from some initial shape of the string, and since he was willing to accept as the 
function f(x) any distribution of values consistent with such a shape physically 
realizable, he maintained the solution of d’Alembert and himself to be the com- 
pletely general one. The values f(x) involved in the formula (4.13) he regarded as 
appropriately subject, if necessary, to graphical definition. To such interpreta- 
tions d’Alembert took exception. He regarded the functional symbol as standing 
for an expression which could be constructed by the ordinary processes of alge- 
bra and the calculus from the independent variables. In having taken the ordi- 
nates of the string to be denotable in the form y(¢, x), he believed that the results 
could apply only to such motions as might be characterized by the fact that in 
them the string shapes at any two instants #, and ¢, are obtainable from one and 
the same formal expression y(t, x) by giving to ¢ the respective values. He saw no 
reason to suppose that all possible motions conform to this. Furthermore, inas- 
much as the differential equation from which the solution emerges involves 
the derivative 0?y/0x?, he was unwilling to admit the applicability of the 
analysis to cases in which the function f(x) is not twice differentiable. Finally, 
because of the relation (4.13) he insisted upon restricting the solution to in- 
stances in which the function f(x) is periodic. That there might conceivably 
exist expressions ®(x) and f(x) yielding the same values over some specific 
interval but not persisting in this relationship for other values of the varia- 
ble, was believed by neither d’Alembert nor Euler nor by any of their contem- 
poraries. 

A difference between d’Alembert and Euler lay in the fact that whereas the 
former was inclined to look upon the concepts of the function and the analytic 
expression as synonymous, the latter would not hold to this. Euler saw no reason, 
for instance, to rule out the possibility of releasing a string from the position of 
a curve made up of circular arcs of different radii, provided these arcs joined 
with each other continuously and with a continuously turning tangent line. It is 
evident that Euler had advanced measurably to the conception of an arbitrary 
curve. It is understandable, however, that d’Alembert should have declined to 
acknowledge the legitimacy of admitting such curves into consideration where 
the operations of the calculus were to be employed. 

In 1755 a memoir of Daniel Bernoulli’s upon the motions of the string 
turned the entire disagreement into new channels. Bernoulli, who had inter- 
ested himself in acoustics, had recognized the relation between the several 
normal vibrations and the respective overtones which the string could be made 
to emit. It was a generally recognized fact at the time that a musical string 
ordinarily responds with a combination of its fundamental and overtones. Ber- 
noulli had discovered that the motions involved in this do, in a very definite 
sense, retain their individuality—that in the entire motion the several normal 
vibrations are simply superposed upon each other. It was a relatively moderate 
step from this to the conception that all possible motions of a string are but 
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linear combinations of the normal vibrations, variations in the relative intensi- 
ties of the overtone components producing the observable differences in the 
timbre of the tone. In terms of symbols, and with the use of the formulas (4.4), 
this comes to its formulation in the assertion that every motion of the string is 
expressible in the form 


vrat 


vax 
(5.1) y(t, x) = A,sin—— cos 
p= 1 
with appropriate constant coefficients A,. 
Neither Euler nor d’Alembert was inclined to accept this, and each made 
his rejoinder. Euler quickly recognized the fact that a motion representable in 
the form (5.1) would be one for which the initial ordinates have the values 


(5.2) 

To assume that all motions are here involved would, he pointed out, come to the 
assertion that an arbitrary function f(x) could be represented by a series of the 
type (5.2). Since among other things any expression (5.2) is odd and periodic, 
it seemed to Euler that he had reduced Bernoulli’s claim to a manifest absurdity. 
That Bernoulli’s result gave solutions—special ones—he did not deny. He 
had, in fact, made that discovery on his own account some years earlier. 

D’Alembert, on his part, not only endorsed all of Euler’s objections but went 
in fact well beyond them. He was unwilling to concede even that any and every 
odd and suitably periodic function could be represented by an expression (5.2), 
maintaining, in particular, that the function would need to be twice differentia- 
ble since that is so of all terms of the series. Bernoulli’s analysis, and especially 
his passage from the finite case of the loaded string to the continuous one, had 
been at best sketchy and fragmentary. His opponents found much that could 
properly be rejected in that. 

On the whole the objections left Bernoulli unshaken. In replying to Euler’s 
claim of absurdity he referred to the fact that any finite sum of m terms from 
the expression (5.2) could, by an appropriate determination of the coefficients, 
be made to coincide in value with any given function f(x) at any chosen set of 
points m in number. He saw no reason, therefore, for rejecting the possibility 
that the series (5.2), involving infinitely many coefficients as it does, might not 
coincide with an arbitrary function at an infinity of points. This viewpoint was 
indeed a worthy one. A development of it will concern us in the following chap- 
ter. 

The three cornered polemic spreads itself through the mathematical litera- 
ture over a period of more than a decade. Since no one of the contenders suc- 
ceeded in convincing another, the upshot of the matter at the time was negligible. 
Each of the disputants was in part right and in part wrong. Time has given the 
lion’s share of its endorsements to Bernoulli. 


: q 
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CHAPTER 6 


Lagrange’s solution of Bernoulli’s problem in curve fitting. [7] The fact 
invoked by Bernoulli, that an arbitrarily given curve is representable at any 
finite set of abscissas by a suitable segment of a series (5.2), is one of considera- 
ble importance both from the theoretical and practical standpoints. It devolves, 
of course, upon the possibility of determining the coefficients c, in a formula 


y= sin 


k=1 


so that the curve here represented may pass through a prescribed set of points 
(x,, F,), v=1, 2, - +--+, (m—1), the abscissas of which lie upon the interval (0, /). 
Alternatively stated it comes to the fact that with an arbitrary assignment of 
constants F, the system of equations 


n—1 k X» 
(6.1) cx sin 


k=l 


= F,, y= 1,2,---, — 1), 


is solvable for the values cx. 

The solution of any linear algebraic system of equations, and hence in par- 
ticular of this system, is, of course, possible by elementary procedures. Such a 
frontal attack upon it by the familiar method of determinants leads, however, 
through much tedious and protracted computation whenever the number of 
equations is large. The solution, moreover, is not likely to emerge from such 
manipulations in any neat or elegant form. At the very beginning of his career, 
while he was still in his early twenties, Lagrange concerned himself with this 
problem and gave solutions of it for both the cases in which the abscissas are 
equally and unequally spaced. It is the former of these which is of peculiar 
pertinence to our discussion, and although it appears in Lagrange’s work inci- 
dentally to the wider investigation with which we shall be concerned in the 
next chapter, a self-contained exposition of it is possible and is to be given 
here. It is a prime merit of this solution that it shows clearly how it depends 
upon the number of points involved, and that it is therefore excellently 
adapted to an investigation in which this number is ultimately to be varied 
and to be allowed to become infinite. 

Let the abscissas x, be identified again thus 


(6.2) 
n 


and let any one of the integers 1, 2, - - - , (2—1) be chosen and designated j. 
If the equations (6.1) are multiplied by the respective constant of an unde- 
termined set D;, and are then added, the result is the relation 
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n—1 n—1 n—-1 
pop asin = 
kel 


Since kx, =vx, this can be given the alternative form 


k=l 
the function ®; herein being defined by the formula 


n—1 


(6.4) ®,(x) = D;,, sin 


Now each function sin (vrx/) in this relation is expressible as the product of 
sin (rx/l) by a polynomial of the degree (v—1) in cos (x/l) (cf. appendix II). 
The complete function ®; may therefore be similarly expressed, namely thus 


x 
(6.5) @,(x) = sin (cos =), 


with P,_: designating a polynomial of the-degree (n —2). The coefficients of this 
polynomial depend, of course, upon the multipliers D;,, and these have not thus 
far been specified. It is proposed now to specify them so that the function 
@,(x) may be zero at each of the points x, with the specific exception of x;, 
namely so that 


(6.6) (xx) = 0, 


Assuming this to be possible, it is clear from the equation (6.5) that each one of 
the (7—2) values cos (1x,/l), v¥j must then be a root of P,_2, and that each 
corresponding difference (cos +x/l—cos mx,/l) must therefore be a factor. 
The factors are thus all accounted for, and with an appropriate constant a the 
formula (6.5) may accordingly be written 


n—1 
(6.7) @,(x) = asin Il (cos-** — cos =). 


As is shown by the formula (II.1), however, a relation 


sin —— = sin — Cos —— 
l 


also maintains, and since the left-hand member of this is zero at each point x, 
without exception, the function p,_:(cos rx/l) must admit as a factor each of 


( 
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the differences (cos rx/l—cos 1x,/l). The factors being thus again all accounted 
for, it follows that 


sin —— = sin — cos —— — COs 
l er l l 


with 6 standing for some constant that is not zero. From this together with the 
evaluation (6.7) it appears that 


(6.8) ( =) &,(2) a 
B 


namely, because of the formula (6.4), that 


> Dj,» sin —( cos — — cos — 
v=1 l l l 


This equation may be reduced by the use of the familiar relation 


sin —— cos —— = — sin ——————- + — sin —————_ 
l l 2 l 2 l 


and by a rearrangement of its terms, to appear in the form 
(6.9) [Dj.n41 qiDj,x ot Dj.x-1] Dj,a-1 0, 
k= B l 
the coefficient g; being specifically 
TX; 


(6.10) = 


and D;,o and D;,, being zero. 
The equation (6.9) is identically fulfilled if the multipliers D;, satisfy the 
system 


Dj,o = 0, 
+ + = 0, k=1,2,---,(#—1), 
Din = 0, 


and furthermore 
2a 


This system is precisely that which is discussed in the appendix I, the coefficient 
(6.10) being that one of the values (I. 5) for which the system admits the solu- 


i 
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(I. 6) with y=j. The free coefficient A; in this may be determined, more- 
over, to yield the value prescribed above for D;,,-1, and is thus found to have the 
value 2a/8 sin (n—1)j/l, or because of the relations (6.2) 


2a 


The evaluation of the multipliers which thus results is 
2a sin —— 

(6.11) D;,. = (— 


Wx; 


and with these the equation (6.3) reduces by virtue of the values (6.6) to the 
form 


n—1 


(6.12) = Dj, 


It only remains, therefore, to determine the value of ®,(x;), and this may be 
done as follows. The formula for ®,(x), as it is given by the equation (6.8), is 
indeterminate at x=x,;. By an application of l’Hospital’s rule, however, its 
limiting value is found to be 


namely because of the definitions (6.2), 


= (— 1) 
UX; 


This result, together with the evaluations (6.11), causes a final reduction of the 
relation (6.12) to the form 


2 n—1 X; 
(6.13) 
n y=1 l 


and therewith the coefficients in the equations (6.1) have been determined. 


@ 
A;= (- 1) 
B sin 
y= 1 
na 
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CHAPTER 7 


Lagrange and the vibrating string. [8] In the controversy over the problem 
of the vibrating string Lagrange was inclined on the whole to enlist himself 
upon the side of Euler. To support himself in this position he undertook to re- 
examine afresh the behavior of the weightless loaded string with an unspecified 
number of particles, his explicit purpose being to elicit from this a proof that in 
the case of the continuous string no restrictions upon the shape of the curve 
marking the initial position are requisite. His method in this has become a 
standard one. As do the deductions of the preceding chapter, it hinges primarily 
upon an introduction of undetermined multipliers. By this means he carried 
through, as we shall see, a general integration of the differential equations for 
the string’s motion, and thus displayed in terms of explicit formulas the de- 
pendence of the string’s position at any instant upon its initial shape. 

The differential equations for the particles of the loaded string of the length / 
under the tension T, with m particles of total mass M located respectively at the 
points 

kl 


(7.1) (ho 
n 


were deduced in chapter 3 and are given under (3.5). If the initial ordinates of 
the particles are denoted by f;, and if the particles spring at =0 from the state 
of rest in these positions, the boundary relations to which the differential equa- 
tions are to be subjected are 
= fas 
(7.2) 
dt 


k=1,2,--+,(#— 1). 


Let the equations (3.5) be multiplied respectively by unspecified constants M,. 
The addition of them then results in the single equation 


n—1 d* na 2 n-1 


k 
k=l dt? l kel 


and this, under a re-grouping of its terms, together with the evaluations 
M,=0, M,=0, takes on alternatively the form 


a? n—1 na 2 n—1 
(7.3) = [Misr — + ye. 


dt? k=l 


Consider now the possibility of so choosing the multipliers M;, as to make the 
corresponding terms of the two sums in the equation maintain a fixed ratio to 
each other. With a constant of proportionality y, the condition upon the multi- 


é 
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pliers is then this, that they may comprise a solution of the linear algebraic 
system 


Mo = 0, 
(7.4) — + Mir = k = 1,2,---, (mn — 1), 
M, = 0. 


This, however, is precisely the system (1. 1) of the appendix I, with (—y—2) 
in the réle of the coefficient g. The values of y for which the system is non- 
trivially solvable are thus found from the relations (I. 5) to be (n—1) in num- 
ber, namely y=¥,, v=1, 2, (a—1), with 


(7.5) = —4sin?—- 
2n 


Upon designating by the symbols M,, those multipliers M, which satisfy the 
system when ¥ has the value (7.5), we find from the formulas (I. 6) the respec- 
tive evaluations 


kor 
(7.6) = A,sin k=0,1,2,:--,m, 
n 


the coefficient A, being arbitrary. 
Let o,(t) be used now as an abbreviation in the sense 


n—1 


k=l 
The differential equation (7.3) with its boundary relations (7.2) may then be 
written in the form 
(“). 


n—1 


o(0) = Mo 


k=l 


= 0 
dt 
This is a differential system which is easily solvable by elementary means. Its 
solution is 


o,(t) = cos (v=7, 


namely, in terms of the evaluations (7.5), (7.6), and (7.7), 


| 
|__| 
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22 2nat vr 
(7.8) > ye(é) sin > fe sin cos ( sin — 
k=l l k=l 2n 
Inasmuch as the index yp is free to take the values 1, 2,---,(m—1), thisisa 


system of (n—1) equations. 

In structure the system (7.8) is evidently of the form (6.1) with the values 
y.(t) in the place of unknowns ¢. Its solution is therefore given by the formulas 
(6.13) to be 


2 x x 2nat 
(7.9) y() =— DD fa sin (= sin =), 
NM k=l l l 2n 
1). 


With this result Lagrange’s integration of the equations of motion is complete. 

It is suggestive for our purposes to consider the formalisms of a passage from 
the discretely loaded string to the continuous one upon the basis of Lagrange’s 
formulas. If in the equations (7.9) the notational changes indicated by the sub- 
stitutions of x and y(#, x) for x; and y,(t), of s, and f(sz) for x, and f;, and of 
As for 1/n, are made, the formula assumes the aspect 


; sin (=) 
y(t, x) = 7 ( f(se) sin as) sin cos 


\ 


As n is indefinitely increased the relation 


sin (=) 
2n 
lim —————- = 1, 
(=) 
maintains for each value of v and from the very definition of the definite integral 


n—1 


lim >> f(sx) sin = As = f f(s) sin - ds. 
0 


k=1 


These relations suggest as the limiting form of the solution (7.9) the formula 


t 
(7.10) y(t, x) = f f(s) sin as) sin cos . 
v=1 0 


It will be noted at once that this is precisely of the type of the solution (5.1) for 


| 
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which Bernoulli had insisted upon holding out. In it the initial position of the 
string is along the curve whose ordinates are f(x). At the instant t=0 the 
formula (7.10) thus reduces to the form 


(7.11) fa (— f as) sin 


The representability of an “arbitrary” function by a series of sines in precisely 
the manner (1.3), (1.4), would hereby seem to be unmistakably presaged. 

Although Lagrange carried out a limiting analysis upon his formula (7.9), it 
varied in some respects from that outlined above. His comparative result was 
the relation 


(7.12) y(t, x) = (x > sin sin cos 


which differs from (7.10) in having the order of the integration and summation 
reversed. This form has the disadvantage of involving a series which is obviously 
divergent, a matter which was readily seized upon by opposing critics. Beyond 
that, both limiting considerations are open to criticism upon a number of ac- 
counts, for they fail to distinguish between the results of analytical operations 
upon an infinite series as a whole and upon the terms of the series individually. 
These distinctions, so essential to rigor, were but imperfectly understood at the 
time. 

While in the manner shown the formula (7.9) could easily have led to the 
conclusion (7.11), it remains a fact that it did not do so. In deducing the form 
(7.12) Lagrange was bent upon a different purpose, wholly remote from that 
of proving any such a theory as would be implied by the relation (7.11). In- 
deed, when such a theory was announced by Fourier more than a half century 
later, the then aged Lagrange is said to have remained incredulous of it. 


CHAPTER 8 


Euler’s determination of the coefficients. In the latter half of the eighteenth 
century the properties of trigonometrical series were very much to the fore of 
mathematical interest, and numerous memoirs were written during that time 
upon one phase or another of the subject of the representability of functions by 
means of such series. In the main, however, these papers advanced the general 
theory but little. They may well be left aside in the present discussion as of 
only subordinate interest. A conspicuous exception to this, however, is a work of 
Euler’s which he appears to have written in the year 1777, although its publica- 
tion was deferred until 1793, some years after his death. Concerned with func- 
tions known upon some grounds or other to be representable in terms of a cosine 
series of the type (1.5), Euler deduced in this work the formula (1.6) for the 
coefficients [9]. 


- 
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If f(x) is any function which in terms of the variable &, under the relation 
t=cos 1x/l, is expansible for —1$£<1 in a convergent power series 
then f(x) clearly admits of representation by a cosine power series of the form 

(8.1) f(x) = ¢; cosi 
j=0 
over the interval (0, /). For each term in this, however, there maintains a re- 
spective trigonometric identity (cf. appendix III), namely 


2F* 


with [7/2] designating the greatest integer not exceeding j/2. The substitution 
of these evaluations into the equation (8.1) and the subsequent collection of 
terms of like character, give the equation formally the aspect 


(8.2) f(x) = = + > a, cos 


This is the type (1.5). It will be clear that a rather extensive class of functions 
f(x) fulfills the assumptions that are basic to this reasoning. 

The argument given, although it is adequate to permit an inference of the 
form of the representation (8.2), is readily seen to be quite far from being prac- 
tical. It yields neither an easily applicable nor a generally lucid method by 
which the coefficients a, therein may be quantitatively evaluated. Since the 
series converges, it may, of course, be inferred that 


(8.3) lim a,=0. 


ve 


From this slender source Euler succeeded in deducing an evaluation of the 
contants a, which is directly referable to the function f(x) in question. 

Let the symbol n,,, for any integer m and any indices ¢, 7, be defined to have 
the values 


n, if o = 7 (mod 2n) 
(8.4) Nor = 
0, if o #7 (mod 2n). 
There is, then, a trigonometric evaluation (cf. appendix IV) to the effect that 


n—1 —1 
(8.5) > [1 + cos or] + 
p=l n 2 


From this a certain related formula can be easily deduced. If the equation (8.5) is 
written successively with o replaced by (v-+) and by (v—&), an addition of the 
results yields directly the equality 
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n 


n 
= — [1+ 4 0s + + 3 cos (v — + + 
and this is contractible in an obvious fashion into the formula 


n—1 kr 
(8.6) >> 2cos cos = — [1+ cos kw cos vr] + + 


u=1 


Consider now the representation (8.2). If this is multiplied by the factor 
2 cos (ukw/n) and is then evaluated at x=yl/n, it yields the equalities 


pl pkr per pkr 
2f “) cos = ay cos + a,2 cos — cos — - 
n n val n n 
When these are summed with respect to u the constants a and a, appear with 


coefficients that are given by the formulas (8.5) and (8.6) respectively. It is 
thus found that 


a2 ~) ites — a,[1 + cos kr cos vr | 
n v=] 


+ aonzo + + my,-x]. 
vel 


This relation can be materially reduced. In the first place it will be found on the 
basis of the definitions (8.4) that for k=0, or for k>0 and n>k, 


veal 


A=1 


In the second place the formula (8.2) itself yields the equation 
#(0) + f(1) cos kr = [1 + cos Ae] + + cos cos 
vel 
With these evaluations, however, the formula (8.7) is simplified into the form 


(8.8) 25 ) cos naz — [f(0) + cos kr] + nd + 


Let the notational substitutions s, =yul/n, As =1/n, now be introduced. These, 
together with a division by n, give to the equation (8.8) the form 


cos As =a,-— ~ + f(l) cos kr] + (Gona—k + Genr+k)- 


l p=1 A=1 


(8.7) 
| 
fe 
4% 
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In this ” is now to be permitted to become infinite. Since each term den,44 in the 
final sum approaches zero by virtue of the relation (8.3), and the term in 1/n 
does likewise, the right-hand member of the equation has as its limit a,. The 
limit of the left-hand member being a definite integral, the conclusion is that 


krs 
(8.9) f T(s) ds = ay. 
0 


Herewith the problem was solved. 

Once in possession of the formula (8.9), Euler recognized the more direct 
manner in which he might have found it, and by means of which a verification 
of it might be made. This is, namely, the procedure now generally familiar, of 
multiplying the representation (8.2) by the factor cos (kmx/l), integrating it 
then term by term and applying the elementary evaluations 


0, if k, 

krx 

(8.10) f cos = 41/2, if »=k#¥0, 
if v=k=0. 


To this day the constants a; as given by the formulas (8.9) are still widely 
known as the “Euler coefficients” of the function f(x). Such an attribution 
seems, however, to be somewhat over-generous if representations (8.2) of arbi- 
trary functions are in question. Euler was consciously concerned only with such 
functions as were known upon other grounds to be representable in a cosine 
series. The crucial observation that the formulas (8.9) are significant for func- 
tions of a much wider class than those which, for instance, are representable in 
the manner (8.1) apparently escaped him. There is no evidence, either in this 
connection or in any other, that he ever receded from his opposition to Ber- 
noulli’s claim that arbitrary functions submit to trigonometric representation 


CHAPTER 9 


Fourier and the theory of heat. In the interior of a material body heat is in 
general distributed in a manner that is both non-uniform and fluctuating— 
that is to say with temperatures that vary from point to point and from time 
to time. The distribution of temperatures throughout a body is, therefore, 
naturally determined by a function of the codrdinates of position and time. 
What the precise form of this function is, in any particular case, depends in part 
upon the thermal properties of the material of which the body is constituted— 
its density, specific heat, and conductivity—but also in large part upon the 
instantaneous state in which the body finds itself at some specific time, and 
upon the conditions which thereafter maintain upon its surface. 


FOURIER AND THE THEORY OF HEAT 31 


In the early years of the nineteenth century Fourier devoted himself to an 
analysis of this temperature function 7, and deduced from physical funda- 
mentals [10] the fact that it must satisfy a partial differential equation of the 
form 


(9.1) = —- 


In this x? is a positive constant whose value is determined by the thermal 
properties of the material, whilc V°r is the so-called “Laplacian of r.” In terms 
of rectangular coérdinates this differential expression is 


dat 
if all three of the coordinates x, y, 2, are significant, or, more simply 


respectively, if the only space coérdinates are x, y, or merely x alone. A tempera- 
ture function must accordingly solve a partial differential equation such as (9.1). 
In any specific instance it must be that solution of this equation which takes on 
those values which apply at some specific instant ¢, and which furthermore ful- 
fills upon the body’s surface the thermal relations that maintain there. 

A simple and familiar physical formulation illustrates this and will serve also 
to show the relevancy of this subject of the flow of heat to the basic matter be- 
fore us, namely that of the representation of arbitrary functions by the means of 
trigonometric series. Consider a homogeneous material bar in the shape of a 
right cylinder of the length 7 and of any cross section. We may choose our co- 
ordinate system so that the direction of this bar is that of the x-axis with the 
end faces of the bar located at the points x =0 and x =/. Let it be supposed now 
that at some specific instant, which may be designated as :=0, the temperatures 
at all points within the bar having the abscissa x have the common value f(x). 
From this instant onward each end face of the bar is to be held constantly at 
the temperature zero, while the lateral surface is insulated againct the passage 
of heat. Except in the trivial case in which f(x) is everywhere zero, heat will flow 
within the bar, and the lines of flow will be parallel to the x-axis. The problem is 
to determine the temperature at any point of the bar at any instant subsequent 
to the initial one, namely to determine the function r(x, ¢) for 0<x </, and t>0. 
The relations from which this is to be done are in this case evidently the dif- 
ferential equation 


: 
or 
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Or Or 
(9.2) 

0x? ot 
the boundary relations 

7(0, = 0, 


(9.3) 
4) = 0, t> 0, 


and the initial condition 
(9.4) r(x,0) = f(x), O<x<l. 


Fourier’s method of attack upon a problem such as this is one that is still 
widely current in practice, namely the method of the “separation of variables.” 
Let us, to begin with, seek a function 7,(x, ¢) to fulfill the equations (9.2) and 
(9.3), while being of the form of a product of a function of the single variable x 
by a function of the single variable ¢, namely 


(9.5) T(x, t) = 


Upon substituting this form into the equations in question, it is found that these 
are satisfied if the function ¢,(x) fulfills the system of relations 


(i) + = 0, 
(9.6) 
(ii) (0) = 0, 
(iii) = 0, 


with X, designating a constant, provided that with this same constant the func- 
tion y,(¢) fulfills the equation 


d 
(9.7) + = 0. 


As a solution of the ordinary differential equation (9.6i) the function ¢,(x) must 
familiarly be of the form 


b, sin Kx + 


in which 6, and y, may be any constants. It will fulfill the condition (9.6ii) if 
,=0, and then also the condition (9.6iii) if 1/\, is any multiple of the constant 
a/xl, namely, when 


(9.8) (=), 


| 
i 
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with vy an integer. The associated solution of the equation (9.7) is then clearly 


Through the formula (9.5) it has thus been found that the equations (9.2) and 
(9.3) are satisfied by the function 


in fact by each of the infinite set of functions obtained from this formula by 
setting yv=1, 2,3,---. 

Now it is a characteristic property of linear equations or of any system of 
such, that the sum of any set of solutions is itself a solution. One is motivated 
thus to infer from the set (9.9) a formal solution having the structure 


x 


Aside from questions of convergence which are clearly to be raised in this con- 
nection, a primary matter still to be dealt with is the fulfillment of the condi- 
tion (9.4) whatever the function f(x) may be. Upon substituting the value ¢=0 
into the formula (9.10), this is seen to devolve into the relation 


(9.11) sin = f(z), O<x<I. 


The question of the representability of any function f(x) in a series of sines is 
thus clearly brought into issue. 

We propose to review in the following chapeere Fourier’s mode of coping with 
this problem. There are a number of reasons why a consideration of this may be 
regarded as worthwhile and of interest. It is, to begin with, ingenious and skill- 
ful. Aside from that it is a notable exemplar of work in the spirit of mathe- 
matics in the eighteenth century. Although this theory of Fourier’s was actually 
created in the next century, and was crowned by the prize of the Academy of 
Paris in 1811, its disregard for rigor was even then outmoded and seemed, in 
fact, to run counter to the better standards of which Fourier himself was 
conscious. It is a formalism—no more—a play upon symbols in accordance 
with accepted rules but without much or any regard for content or significance. 
As such it has, of course, no place in the mathematics of our time. Fourier’s 
work has had the profoundest effect both upon the development of pure mathe- 
matical concepts and upon the extension of the range of mathematical applica- 
tions to the sciences and technology. These deserts, however, sprung in the main 
from Fourier’s interpretations and not from his manipulations. It was, no doubt, 
partially because of his very disregard for rigor that he was able to take con- 
ceptual steps which were inherently impossible to men of more critical genius. 


~ 
y=1 
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CHAPTER 10 


Fourier’s formal solution of his problem. [11] In the study of a representa- 
tion (9.11) it comes to a mere matter of the choice of a unit of measurement to 
identify the length / with the value 7. We shall suppose this to have been done 
since an appreciable formal simplification results from it. To summarize the 
problem at issue, then, it is that of determining from any given function f(x) 
a set of constants b, such that 


(10.1) :¥ b, sin vx = f(x), 


for 0<x <7. 
Fourier began his considerations of this relation by substituting in it for each 
sine function its power series equivalent, namely 
(- 1) 


(2n — 1)! 


Upon interchanging the order of the summations, an operation which was at 
that time generally resorted to without question, the relation (10.1) was made 
to appear in the form 


(- ( ) 

Ke) (2n — 1)! 

The function f(x) has thus been related to a series in powers of x, and since such 

a series must, in fact, be its MacLaurin series 


2. 


in which f!!(x) stands for the kth derivation of f(x), it was to be concluded by a 
comparison of the coefficients of like powers of x that f!*!(0) is zero whenever k 
is even, and that otherwise 


(10.3) = (— n=1,2,3,-°°. 


The effect of this consideration has thus been to throw the constants b, into the 
réle of the infinitely many unknowns in a system of infinitely many linear 
equations. 

To deal with a system of this type Fourier had to invent his own method. 
He chose to base this upon the use of a chain of ordinary algebraic systems 


y=1 
| 
y=l 
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of arbitrarily large degree r, the relevancy of these to the infinite system (10.3) 
to be assured by an adoption of the relations 


(10.5) on( ©) 1, 2, 3, 


Inasmuch as the solution 8,(r),v=1, 2, ---,7, of the system (10.4) is unique, 
whatever the degree r may be, it was tacitly assumed that the solution of the 
infinite system could be inferred from those of the finite ones in the manner 


(10.6) B,(©) = b,, y= 1,2,3,--°. 


In the system (10.4) let the mth equation be multiplied by r? and let the next 
following equation then be subtracted from it. If this is done for each value of 
n from 1 to (r—1) the result is the system of equations 


(10.7) [r? — v?]8,(r) } = 7°,(r) — m= 1,2,°-°+,(r— 1). 


This suggests imposing upon the members ¢,(r) the relations 


for, since the coefficients of the system (10.7) will then be precisely those of the 
system given by the equations (10.4) when r is replaced by (r—1), it may be in- 
ferred that the unknowns (r?—v?)@,(r) in the one case and 8,(r—1) in the other 
case, are, in fact, the same, namely that 


B,(r) = Bor — 1) 


r? yp? 
By suitable iterations of this relation it is evidently found that 


(10.9) Br) = — r>», 


II 


whence it follows by a simple formal step from the finite r to the infinite, that 


B,(v) 
I (nt 


n=y+1 


(10. 10) b, = 


To any modern investigator this conclusion could, of course, be only mean- 
ingless, for the infinite product involved in it is manifestly divergent. Nor is the 
source of this unfortunate result difficult to trace. It lies, namely, in the naive 
adoption of the relations (10.8) in the face of the fact that these are quite in- 
consistent with the previously adopted assignments (10.5). It requires no keen 
critical faculty to observe at once that under the formula (10.9) the sequence 
of values 8,(r) inevitably converges to zero with 1/r and hence that no conclu- 


= 
| 
‘ 
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sion other than that of the vanishing of each coefficient 6, would be logically 

admissible. Fourier had no intention whatsoever of drawing that conclusion, 

and hence proceeded undismayed with the analysis of his formula. In so doing 

he found, quite naturally, that the divergence it already involved could be 

formally compensated for only by the introduction of still other divergencies. 
Let the determinant D(x) be defined by the formula 


i 2  +++(v—1) x 
13 23 (v x3 
(10.11) D(x) =| 15 25 «++ (v— 1)5 x5 


and let the cofactor of its element in the 7th row and jth column be denoted by 
D;,;. Then since the determinant of the system (10.4) when r =r is precisely D(v), 
it is seen at once that by Cramer’s rule 


( ) B,(v) Do) 


This formula can be made much more explicit. The determinant D(x) is, in the 
first place, found to admit (cf. appendix V) of the evaluation 


(10.13) D(x) = (- 1), I (n? — x?). 


n=1 


Upon expanding the left-hand member of this equality by the elements of its 
last column, and by agreeing to define the coefficients c,(r) for any value of r not 
less than v by the relations 


(10.14) II (n? — x’) = Dd ca(r) 
n=1 
it evidently follows further that 
> Da = (— 1)D,,, >. Ca(v) 


n=1 n=1 


Now in any identical equation between polynomials or power series, the 
coefficients of like powers of x in the two members of the equation must be the 
same. It will be seen at once that because of this the substitution of any quantity 
whatsoever in the place of any specific power of x will not destroy the equality. 
In the equation above, therefore, the replacement of x?*-! for each value of n 
by the respective quantity ¢,(v) is legitimate. It leads to the result 


4 
. 
r r 
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n=1 n=1 


and this, together with the value of D(v) that is obtained from the relation 
(10.13), yields upon substitution into the equation (10.12) the formula 


(10.15) B,(v) = 
“IT (n? — v?) 


Let it be observed now that the removal of the factor (r?—x*) from the left- 
hand member of the equation (10.14) has the mere effect of reducing the index r 
to (r—1). From that relation it is thus seen that 


Ca(r)x2*-? = (r? — > aly — 
n=1 


n=1 


namely that 


If in this each power x?’ is replaced by the respective value ¢;,:(r), it follows 
because of the relation (10.8) that 


n=1 n=1 


The sum on the left of this equality is thus independent of r, and this having 
been established it is a simple formal step to write 


n=1 n=1 


It was familiar in Fourier’s time (c.f. appendix VI) that 


= x sin rx 
n=l n? 


‘a relation which in a formal sense, if in no other, yields the formula 


II (nt 24) = ( Il (1-2) sin 


Upon replacing the left-hand member of this by the equivalent series (10.14) 
and substituting on the right the series equivalents 


n=1 n=1 
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sin 7x (— 1)" 


2k—-2 


it becomes 
Il > > (- 
n=1 q=0 k=1 = 1)! 


The procedure of replacing each power x?/-? by ¢,() leads from this to the 
equation 


or, by virtue of the relations (10.5) and (10.16), to the relation 


With this result the formula (10.15) for 8,(v) becomes wholly explicit, and 
through it the formula (10.10) assumes the form 


(2k — 1)! 
» TL 


n=1,n¥#v 


(10.17) b= 


Inasmuch as all quantities in the right-hand member of this are to be regarded 
as known when f(x) is known, this result amounts, at least in a formal sense, to 
an evaluation of the constant B,. 


CHAPTER 11 


The reduction and interpretation of the solution. [12] The result (10.17), 
although it formally accomplishes the task originally set, namely the construc- 
tion of a formula through which the coefficients }, are expressed in terms of the 
given function f(x), will nevertheless hardly be found completely satisfying. For 
the purposes of practical calculation, namely, some reduction of its intricacy 
would clearly be imperative, and Fourier, realizing this, turned his attention to 
its simplification. His effectiveness in achieving this is eloquent commentary 
upon his skill in analytical manipulations. 

Consider, to begin with, the product of the values (m*—v?) which the formula 
contains. The first (v—1) factors of this are negative and have as their product 


v q=0 v q 
a 
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(- @—») + 2). 


n=1 n=1 


The product of the remaining (positive) factors may be written in the form 


n=r+1 n=r+1 
and thus, by simple changes of the index in each of the partial products, the 
entire expression may be made to appear as 


v—1 


n=1 n=r+1 n=1 n=2r+1 


Since in this each natural integer except v and 2» occurs twice, the excepted 
ones occurring just once, an alternative form for the product is evidently 


It has been found thus that formally 


2 


Consider now the formula 


(11.2) f(x) = x 


which is the equivalent of the relation (10.2) by virtue of the fact that in the 
latter each coefficient fl*!(0) with an even index k is zero. A 2q-fold term by 
term differentiation leads from this to the companion formula 


», (2n — 2q— 1)! 


and if in this x is given the value 7 and the index of summation is suitably 
changed, the result is the equation 
2 


(11.3) = flal(z), 


The substitution of the evaluations (11.1) and (11.3) into the formula (10.17) 
causes a reduction of this latter to the form 


(11.4) 


L 
in 
v n=1 
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The infinite series which still appears in this result suggests the formal 
definition of a function u(x) by the relation 
(= 
u(x) = ‘ 
Upon differentiating this twice term by term and thereupon adjusting the index 
of summation, it is found that 


(x) p> 


and hence that 
+ v?u(x) = vf(x). 


This is a differential equation of an elementary type. Its reduced equation is 
solved by the function (c; sin yx+c, cos yx) and from this fact the method of 
“variation of parameters” [13] leads readily to the conclusion that the equation 
itself has a solution 


U(x) = Jf 1 sin v(x — s)ds. 


This is, of course, verifiable at once by direct substitution into the differential 
equation. Upon setting x =7 it is thus found that 


(11.5) U(r) = (- f sin ysds. 


Let the integral in this formula now be integrated by parts 2n times in 
succession, the trigonometric factor being each time the one to be integrated. 
Since, as may be seen from the formula (11.2), every even ordered derivative of 
f(s) is zero at s=0, while the function sin vs vanishes at both s=0 and s=7, the 
result of these integrations is the relation 


U(r) = > + f'?"+21(s) sin ysds. 


By the step from the finite » to the infinite the formal relation 


(= 
may be drawn therefrom, and this together with the formula (11.5) yields the 
equality 

— > = ( 


f f(s) sin vsds. 
0 


| 
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By virtue of this the formula (11.4) is now once more and finally reduced to the 
form 


(11.6) b= =f sin vsds. 


This is, of course, the familiar formula in virtue of which the coefficients }, are 
generally known as those of Fourier. It is the formula (1.4) for the special case 
in which /=7. 

Elegant though the conclusion (11.6) unquestionably is, the verdict of any 
critical appraisal of Fourier’s accomplishment to the point of its derivation must 
inevitably be profoundly disappointing. As to the result, in the first place, that 
was not new. It had been contained in the mathematical literature for over a 
decade—to be precise, since the publication of the memoir of Euler that was 
discussed in chapter 8. While, to be sure, Euler’s results applied only to func- 
tions of a certain class, that is no less true of Fourier’s, since his deductions 
were based upon such material restrictions as the representability of the func- 
tion f(x) in power series of the form (11.2) that converge when x=7. 

Nor could any advantage be claimed by Fourier in the matter of method. 
On the contrary—and even leaving aside the important fact that by its em- 
ployment of divergent processes it divested itself of all rigorous validity—the 
method of Fourier suffers in almost every respect by comparison with that of 
Euler. The device of referring the problem to a system of linear equations, 
ingenious though it is, is nevertheless quite foreign to the nature of the problem. 
The trigonometric functions are conspicuously endowed with many peculiar 
properties and fulfill a great many characteristic interrelationships. Of this im- 
portant fact Fourier’s approach in no way avails itself, while Euler’s, by contrast, 
exploits it to the utmost. In cutting directly to the heart of the matter Euler 
thus attained his result more perspicuously and incomparably more cheaply. In 
this respect the superiority is all his. 

With the priority and preference in manipulative matters thus denied him, 
Fourier’s claim to renown must be based upon other grounds, and these are, 
namely, those of interpretation. Approaching the formula (11.6) afresh, without 
regard for the manner of its derivation, it was observed by him that through it 
each coefficient 6, admits of interpretation as the area between the abscissas 
x =0 and x=7, and under the graph 


(11.7) y= See sin yx. 


Such an area is evidently conceivable, and retains its clear-cut significance, in 
association with functions f(x) that are in a very general sense quite arbitrary. 
Certainly these functions need not be assumed to be continuous or expansible 
by any simple analytical formulas. They might be graphically defined and 
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might well represent distributions of functional values that are extremely er- 
ratic. On the basis of such considerations Fourier concluded that any and every 
function f(x) had associated with it a set of constants B,. 

From this fact alone it would not follow, of course, that with such coefficients 
the representation of the function f(x) by a series (10.1) would result. As has 
been seen in chapters 1 and 5, the masters of the eighteenth century had re- 
jected such a possibility as manifestly absurd. In this matter, however, Fourier 
was willing to disregard opinions and precedents, however well established, and 
to look further for himself. From calculations of the coefficients 6, with small 
indices v in the cases of a great variety of functions f(x), and from subsequent 
plottings of the respective initial segments of the resulting trigonometric series, 
he came to convictions upon two salient points, namely: (i) that the series (10.1) 
always represents the function over the interval 0 <x <7, and (ii) that in general 
this representation does not persist for values of x outside that interval. 

Fourier’s announcement of these facts was quite generally met with in- 
credulity. Even the mass of his substantiating evidence won, in many cases, only 
grudging and reluctant acceptance. The implications behind the new assertions 
were too revolutionary to be easily assimilated. They called for no less than a 
fundamental revision of many concepts that were wholly traditional, some of 
them lying at the very basis of mathematical analysis. On the other hand 
Fourier’s new theory did now finally vindicate the half century old reasoning of 
Daniel Bernoulli by which he had convinced himself, if no others, that any curve 
from which a taut elastic string could spring into vibration could be represented 
by a trigonometric series. 


CHAPTER 12 


The Dirichlet integrals. Once Fourier had deduced the formula (11.6), he, 
like Euler before him, observed that in a schematic way the result is recoverable 
in a most direct and simple manner by the mere expedient of multiplying the 
relation (10.1) through by sin ux with any natural integer He, and then integrating 
term by term over the interval (0, +). The infinite series reduces under this 
process to a single term, because of the evaluations 


(12.1) f sin vx sin uxdx = 0, fory ¥ yu, 
0 


and the formula for b, thus emerges. The procedure is obviously adaptable also 
to the case of a cosine representation 


(12.2) fix) = 
y=l 


the evaluations 


i 
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(12.3) f cos vx cos uxdx = 0, for y ¥ n, 
0 
leading in this case to the formulas 
2 
(12.4) 6, = f(s) cos vsds, y=1,2,3,--:. 
0 


An extension of these results to the case of a function that is given over the 


larger interval (—7, m7) is of importance and is easily deduced. Through the 
relation 


f(x) = fo(x) + f.(x), 
with 

fo(x) = — f(— 2)], 

f(x) = + f(— 2)], 


the function f(x) is expressed as the sum of two components of which the first 
one is an odd function and the second one even. Now the relations (10.1) and 
(12.2) for these respective functions, namely 


fo(x) = > b, sin vx, 


a 
f(x) = + a, COS vx, 
obviously remain quite unchanged if x is replaced by —x. Any validity they may 
have over the interval (0, 7) therefore implies the same over the larger interval 


(—7, 7). It follows at once that for the originally given function the representa- 
tion in question is of the form 


(12.5) f(x) = > + > [a, cos vx + b, sin vx]. 


The formulas for the coefficients in this are, moreover, easily freed from 


reference to the components f(x) and f.(x), and may thus be brought to expres- 
sion directly in terms of the function f(x) itself. Since 


j fo(s) sin vsds = f sin vsds, 


as may easily be verified, and also 


fos (s) cos ysds = f(s) cos vsds 
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the formulas in question are, namely, seen to be 


1 
—f cos vsds, 
(12.6) 


1 
b, =— f f(s) sin vsds. 


In the instance that the interval (—7, 7) which has been taken to be basic is 
replaced by the more general interval (—/, /), the representation (12.5), (12.6), 
is correspondingly replaced by that given in the formulas (1.1), (1.2). 

Let the sum of the first 2N+1 terms of the series (12.5) be designated by 
Sy(x), thus 


N 
(12.7) Sy(x) = = + >> [a, cos vx + b, sin vx]. 


The substitution of the values of the coefficients (12.6) into this gives it the 
aspect 


Sy(x) = =f E + 00s v(s — nds], 


and this may be contracted by the use of elementary trigonometric relations 
(cf. appendix VII) into the wholly compact form 


sin [(V + 4)(s — 
sin [3(s — x) ] 


An equivalent manner of writing this, and one that has some analytic ad- 
vantages, is 


1 
(12.8) Sy(x) = f(s) 


(12.9) Sw(x) = -. © ds + Fox. x, N)ds, 
with the symbols F(s) and W(s, x, N) having the significance 
_ 1 


x, N) = sin [(N + 4)(s — a) 1. 
From the latter of these formulas it is easily inferred that 


| W(s, x, N)| $1, and 


f W(s, x, N)ds 


» for—-wSa<BSr. 
2N+1 
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The function (s, x, N) thus possesses the properties: 

(i) that it is bounded uniformly as to N and s; 

(ii) that its integral over any sub-interval of the range (—7, 7) converges 
to zero with 1/N uniformly as to the sub-interval. 
These properties are sufficient to insure [14] the relation 


(12.10) lim "F(s)U(s, x, N)ds = 0, 


N-o 


whenever the function F(s) is integrable (in the sense of Lebesgue) over the 
interval (—7, 7). 

Now from the definition of the function F(s) it may readily be seen that 
integrability is assured to it by that of the function f(s) provided the point x is 
in the interior of the interval (— 7, 7). Thus for every integrable function f(x) the 
final integral in the relation (12.9) converges to zero, and that relation thus 
implies that 


N + 4)(s — x)] 
AY 
(s — 2) 


whenever the right-hand limit involved in this exists. This permits us at once 
the conclusion: that any function f(x) which is integrable and which is further- 
more such that for it 


si 
(12.11) lim Sy(2) = lim — [( 


(12.12) 
(s — x) 


is a function which is representable by a trigonometric series in the manner 
(12.5), (12.6). 

Although Fourier made calculations upon many specific cases, he gave no 
general proof of his ultimate assertion of the trigonometric representability of 
arbitrary functions. Indeed no such proof could have been given, since in the 
omission of all qualifications upon the functions the assertion is too broad to 
be valid. The first proof that was both satisfactory in the matter of rigor and 
ample in the matter of generality was given by Dirichlet in the year 1829. In 
the manner that has been indicated above, this proof was based by Dirichlet 
upon an establishment of the relation (12.12). The integrals involved in that 
relation and in (12.8) are accordingly known generally as “Dirichlet integrals.” 
Dirichlet’s proof in its original form, or as it has been improved and refined, is 
to be found at many places in the mathematical literature [15]. We shall, 
therefore, go no further into it here but shall draw this part of the discussion 
to its close. In the following part the primary subject of study is to be a gen- 
eralization of the entire theory in which the representations of functions in 
trigonometric terms sink to the status of special cases. 


ds = f(x), 
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PART II 


CHAPTER 13 


The differential boundary problem. Among the most powerful of mathe- 
matical means for the formulation of natural laws are the linear differential 
equations of either the partial or the ordinary type. Varied and diverse as 
physical phenomena certainly are, they nevertheless submit quite generally to 
description by such equations. The flow of heat in material bodies and the 
vibrations of elastic strings under tension are instances of this that have already 
been noted. 

By its very nature as the formulation of a more or less general law, any par- 
ticular differential equation applies, of course, to the entire category of mani- 
festations which the law itself governs. To single out for description any indi- 
vidual phenomenon from such a category, therefore requires that some auxiliary 
means beyond the equation itself be resorted to. This ordinarily takes on the 
form of a set of one or more restricting relations that are expressive of the 
characterizing initial or boundary conditions. The differential equation to- 
gether with such relations is commonly designated as a differential system. It is 
also said to define a differential boundary problem. The system consisting of the 
equations (4.6) and (4.7), for instance, thus defines a partial differential 
boundary problem, namely the one which is descriptive of a certain stretched 
string vibrating with specified end points and initial position of release. The 
equations (9.2), (9.3) and (9.4), define a similar problem, one which describes 
a linear flow of heat from specific initial temperatures and under certain bound- 
ary conditions. 

The method that was employed in the reduction of the boundary problem 
of chapter 9 was designated there as that of the “separation of variables.” Its 
immediate effect was to refer the partial differential system to an ordinary 
system (9.6), this latter having the peculiarity of involving an unspecified con- 
stant or parameter, which was designated by A,. Following the solution of this 
ordinary boundary problem at characteristic values of this parameter, the 
theory led in a natural way to the further problem of representing an arbitrary 
function in terms of the respective solutions. This method was in no way espe- 
cially designed for the problem of chapter 9. It is on the contrary one that is 
highly flexible and of very wide applicability. In the particular instance there 
considered the ordinary differential equation which characterized the problem 
was one whose solutions were trigonometric functions, and it was because of 
that, that the representation of the function f(x) took the form (9.11), namely 
that of the Fourier theory. This feature was special, to the extent that it would 
not even have maintained for the equation (9.2) if the coefficient x? involved in 
it had been dependent upon the coérdinate x, rather than constant. In the fol- 
lowing a discussion is to be framed which is free from such peculiar specializa- 
tions. 
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Let the variable x be real, with the range 
(13.1) eszs 5, 


and on this interval let p(x), g(x), and r(x), be differentiable functions. The sym- 
bol L(¢, \) is to designate the differential expression 


(13.2) = + p(x)o’ + + r(x) 


In this \ is to play the réle of a parameter, the range of which is to be the 
entire complex plane. The differential equation 


(13.3) L(g, ) = 0, 


is, then, one that is regular, in the sense that it has no singular points upon the 
interval (13.1). As an equation of the second order it will, of course, not gen- 
erally be explicitly solvable. Certain facts concerning its solutions are, however, 
familiar. Of these the following ones will be especially relevant to the discussion 
proposed [16]. 

(i) The equation admits of solutions ¢(x, \) that have continuous second 
derivatives as to x, and that are analytic in \ over the entire complex \ plane. 

(ii) There is a pair of such solutions ¢1(x, A), d2(x, A), that are linearly inde- 
pendent as functions of x for all values of X. 

(iii) The Wronskian Q(x, A) of this pair, namely the determinant 


o1(x, d) $2(x, d) 


13.4 Q(x, = ’ 


is subject to the relation 
(13.5) Q(x, = Q(a, 42, 


with Q(a, A) an analytic function of \ that is different from zero for all X. 
(iv) The general solution of the equation (13.3) is expressible in terms of the 
pair $1(x, A), d2(x, A) by a form 


(13.6) hapi(x, A) + d), 


in which the coefficients h;, h, are constants as to x, though they may be func- 
tions of X. 

For any equation (13.3) there are known to be infinitely many pairs of solu- 
tions ¢1(x, A), d2(x, A), that have the properties enumerated. Of these pairs any 
one serves in every way as well as any other, and the choice that is made is 
accordingly immaterial. For the sake of avoiding gratuitous complications, how- 
ever, it will be supposed throughout the discussion that when a choice of such a 
pair in the instance of any specific equation has been made, it will be con- 
sistently adhered to. To that extent, then, the designations ¢:(x, A), d2(x, A) 
will be understood to apply not to random but to specific solutions. 

Let the symbols §;,;, ¥;,;,4=1, 2;7 =1, 2, 3, 4; now denote any constants which 
are such that when 
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(13.7) = Bid + 
then the matrix 

@2,1(A) a@2,2(A) 


is of the rank 2 for every value of \. In terms of the values of any function @ and 
its derivative 6’ at x =a and x=), the forms 


(13.9) A) = + ai,2(A)O(a) + + = 1,2, 


are then always linearly independent. Under these circumstances the dif- 
ferential system 


(13.8) 


L(u, \) = 0, 
(13.10) Ai(u, A) = 0, 
Ao(u, = 0, 


defines an ordinary boundary problem. It is this problem which is to be central 
to our discussion. It will be seen at once to include as a special case the problem 
(9.6), and to do that even if in the latter the coefficient x? varies with x. Many 
other boundary problems that stem from physical origins are also included, as 
will upon occasion be seen in the following. 


CHAPTER 14 


The characteristic values and solutions. Since any solution u(x) of the differ- 
ential system (13.10) must in particular solve the differential equation (13.3), 
it must have the form 
(14. 1) u(x) = hodi(x, + hid2(x, d). 


Except for the trivial solution u(x) =0, which we shall herewith specifically and 
permanently rule out of this discussion, the values fy, 2 will not both be zero. 
The substitution of this form into the boundary relations of the system give to 


the latter the aspect 

h hyA1,2(A) = 0, 

(14.2) 2A1,1(A) + 41A1,2(A) 
+ h1A2,2(A) = 0, 


in which the abbreviations 


(14.3) Aj, = Ai($;(x, d), d), i, 1, 2, 


have been resorted to. 
The equations (14.2) constitute an algebraic system in which the values 
hy, fe, function as the unknowns. Since this system is homogeneous its non- 


trivial solvability is contingent upon the vanishing of its determinant A(A), 
where 


; 
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_ | 41aQ) 4120) 
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A proper solution of the boundary problem thus exists if and only if \ is a root 
of the so-called characteristic equation 


(14.5) A(A) = 0. 


(14.4) A(A) 


These roots are called the characteristic values (Eigenwerte) of the boundary 
problem. The multiplicity with which such a value occurs as a root of the equa- 
tion (14.5) is also designated to be its multiplicity as a characteristic value. If 
the root is one at which the elements of the determinant (14.4) do not all 
vanish, namely at which the rank of the determinant is 1, it is said to be a char- 
acteristic value of the index 1. On the other hand a value at which the elements 
do all vanish is said to be of the index 2. It is not difficult to see that in this 
latter case its multiplicity must also be at least 2, and hence that in every case 


(14.6) (The index) S$ (The multiplicity). 


There is material advantage to be gained by formally regarding a value whose 
index is 2 as being, in fact, two coincident characteristic values. We shall here- 
with, once for all, adopt this convention. 

The determinant A(A) is an analytic function for all values of \. The number 
of its zeros in any finite region of the \ plane is, therefore, finite. Thus, in par- 
ticular, only a finite number of these zeros fulfill a relation |\| <N, whatever the 
constant N may be, a fact from which it follows that they—the characteristic 
values—may be sequentially ordered in a succession of non-decreasing absolute 
value. With the assignment of subscripts in such a succession, the characteristic 
values thus follow each other in the array 


(14.7) Ao, A1, As, 
In this, then, 


for every subscript 7. Each characteristic value of the index 2 occurs in the 
array twice, occupying, as we may and shall assume, two consecutive positions. 
Each characteristic value of the index 1 occurs just once. 

Consider now any characteristic value A, that is of the index 1. The pairs of 
values =(h%, which satisfy the system (14.2), namely for which 


he + hy = 0, 


(14.9) 
he + hy = 0, 


have in this case members that stand in a fixed ratio to each other. Any such 
pair yields through the formula 


(14.10) u(x) = he Ae) + hy oo(x, dr), 
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an associated solution for the boundary problem, and the solution so obtain- 
able from any other pair is a mere multiple of this. The function u,(x) is known 
as a characteristic solution (Eigenfunktion) of the boundary problem. 

If now, on the other hand, A, and A,41 are the designations of a characteristic 
value of the index 2, then at this value the system (14.9) is vacuous since each 
quantity A,,;(A,) is zero. There is, therefore, no restriction upon the choice of 
hy, he, and hence in particular two pairs that are not multiples of each other may 
be taken. These lead through the formula (14.1) to two linearly independent 
characteristic solutions u,(x), u,41(x), which we may regard as associated with 
the values A, and A,4: respectively. Thus in every instance each symbol A, has 
associated with it a function u,(x). 

There exist boundary problems that have only a finite number of character- 
istic values or even none at all. There also exists, however, a large class of such 
problems for which the characteristic values are infinitely numerous. It is only 
to problems of this latter category that the continuing discussion will, in all its 
phases, be relevant. 


CHAPTER 15 


The adjoint boundary problem. With the coefficient functions p(x), g(x), and 
r(x), that occur in the expression (13.2) let M(@, \) be defined by the formula 


(15.1) M(¢, = 6” — + [g(x)d + r(x) ]o. 


This differential form is said to be adjoint to the form L(@, \). It is customary 
also to refer to it as the adjoint of L(g, d). If it is completely written out, thus 


— p(x)d’ + + r(x) — 


and its adjoint is in turn constructed, this latter is found to be again the form 
L(g, d). The relationship of being adjoint is thus a reciprocal one, either of 
two forms so associated being the adjoint of the other. For the two adjoint 
forms the equality 


d 
(15.2) ¥L(¢, — A) = x), 
with 
(15.3) Q(d, ¥, x) = o'(x)W(x) — + p(x)o(x)¥(4), 


is easily verified for any two suitably differentiable functions ¢ and y. The rela- 
tion is thus an identity. It is generally known as the “Lagrange identity,” and is 
the source of many important analytical formulas. 

The notion of the adjoint relationship is extensible in the most direct and 
immediate way to differential equations such as the equation (13.3). This latter 
and the equation 


i 
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(15.4) M(¥, ) = 0, 


are thus likewise said to be adjoint. As the discussion proceeds there will be 
ample illustration of the manner in which the differential equations of an ad- 
joint pair, or their respective solutions, interplay in the development of a 
theory. Even here it may be observed that the solubility of either equation 
implies that of the other, since the solutions of either are simply expressible in 
terms of those of the other. Thus if the solutions ¢;(x, A), 2(x, A) of the equa- 
tion (13.3) and their Wronskian (13.4) are used to construct the functions 
), A) by the formulas 


_ d) 


these latter are found by direct substitution to satisfy the equation (15.4). They 
are similarly seen to fulfill the relations 


(15.6) (i 
Wi, = 0, if j ¥i, 


and to be linearly independent. In terms of them the general solution of the 
equation (15.4) is accordingly expressible in the form 


(15.7) v(x, A) = Ropi(x, A) + Rivo(x, d), 


with coefficients ki, ke that are free from x. The relations (15.6) when applied 
to the form (15.7) yield at once the evaluations 


QO(¢1, vy, x) ko, 
2, x) = hi. 


The extension of the notion of the adjoint relationship to the complete 
boundary problem (13.10), although it is not immediate is nevertheless possible, 
and may in fact be made in several formally different but essentially equivalent 
ways. We shall do this in the following manner. Consider the differential system 


(15.8) 


M(v, ») = 0, 
v(@) = + 
(15.9) — v'(a) + p(a)v(a) = + 


— 0(b) = + wiae,s(A), 
v'(b) — p(b)0(b) = poars(d) + (A), 


in which there occurs besides the function v(x) also a pair of “parameters” yu, ue 
that are independent of x. The coefficients a;,;(A) are to be those which were 
defined by the formulas (13.7). We shall show that this system in fact defines a 
boundary problem which is essentially of the type (13.10), and that the param- 


are 
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eters j41, M2 may be looked upon as standing for certain specific linear forms in 
the values v’(a), v(a), v’(b) and v(b), with coefficients that are functions of X. 

The condition that the set of boundary relations of the system (15.9) be 
consistent in the “unknowns” yy, pe, is that the matrix 


v(@) a1 

— (a) + plajo(a) a2 
— 06) 

v'(b) — p(b)v(b) aia 


be of the rank 2. Now by hypothesis there is, for every 4, some two rowed 
determinant from the last two columns of this matrix that is not zero. This 
occurs as a minor in two of the three rowed determinants of the matrix. The 
results of setting these latter equal to zero are two equations in the quantities 
v’(a), v(a), v’(b), v(b). That these equations are independent follows at once 
from the fact that each contains one of the four quantities which the other does 
not contain. Thus the differential equation of the system (15.9) is seen to have 
imposed upon it two linear boundary conditions. A boundary problem is thus 
defined. The boundary relations of this problem evidently have coefficients 
that are polynomials in \. To this significant extent the problem is accordingly 
similar in form to the problem (13.10). It is true that the coefficients of the 
boundary relations of the system (13.10) were taken to be linear polynomials, 
whereas those of the newly found system may be quadratic. That, however, is 
not truly important, for the assumption of the coefficients of the problem (13.10) 
to be of the first degree in \ was motivated only by the desire for simplicity, 
and is in no way essential. Finally some pair of the equations (15.9) is always 
solvable for uy; and pe. By that solution mw; and pe are expressed as linear forms in 
the values v’(a), v(a), v’(b), v(b) as was asserted above to be possible. 

It is a matter to be observed that if in the equations (15.9) the parameters 
1, Me were to be both zero, it would follow that v’(a) and v(a) would also neces- 
sarily vanish. Only the trivial solution of the differential equation M(v, \) con- 
forms to these values. Since this solution is to be ruled out of the discussion, it 
is evident that the simultaneous vanishing of wi and pe is likewise to be barred. 
It is to be understood henceforth, therefore, that of the value pair 41, we at least 
one member is in every case different from zero. 

If in any differential system the differential equation is replaced by its gen- 
eral solution, and the boundary relations are replaced by independent linear 
combinations of them, the content of the system clearly remains unchanged. In 
consonance with this let the first two of the boundary relations of the set (15.9) be 
multiplied respectively by the factors ¢/ (a, A) and ¢,;(a, A) with 7 =1, 2, and let 
them then be added. The left-hand members thus obtained are found to be 
Q(¢;, v, a), and thus by the relations (15.8) the resulting equalities are 


(15.10) = + j= 1,2, 
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with 
(15.11) (A) = (a, A) + d), i,j = 1,2. 


In a similar way the last two of the relations (15.9) may be combined with 
the multipliers —¢j/(b, 4) and —¢@;(b, to assume the forms 


(15.12) = — + 42,0], 7 = 1,2. 
In content the equations (15.7), (15.10) and (15.12) are thus equivalent to the 
set (15.9). 

Consider now the case in which J is any characteristic value of the boundary 
problem (15.9), namely a value for which a set of elements v(x), 1, me, fulfilling 
the equations (15.9) exists. From the fulfillment of the equations (15.10) and 
(15.12) it follows then at once that 


1,1(A) + H1A2,1(A) = 0, 
+ = 0. 


Since w; and pe are not both zero, the determinant of this system must vanish. 
This determinant is, however, precisely A(A), as that was defined by the formula 
(14.4). Thus A must be a root of the equation (14.5), namely a characteristic 
value of the boundary problem (13.10). Every characteristic value of the adjoint 
boundary problem is thus also a characteristic value of the given one. 

The converse of this may also be established. Thus let A, be any value from 
the set (14.7). With \ at this value the system (15.13) is non-trivially satisfiable 


by values pu”, wu”, which, of course, fulfill the relations 


(15.13) 


+ = 0, 


15.14 n 


Let v,(x) be the solution of the first three equations of the set (15.9) with A=A, 
and wj;=y\”, j7=1, 2. The familiar “existence theorem” for a linear ordinary 
differential equation [17] gives assurance of both the existence and the unique- 
ness of this function. Now with the values at hand the equations (15.10) are ful- 
filled. Because of this, however, and with the equations (15.14), the relations 
(15.12) are seen to be also fulfilled. This means that the entire system (15.9) 
admits this solution, namely that X, is also a characteristic value of the boundary 
problem (15.9). We may go even somewhat further. If the index of \, relative 
to the boundary problem (13.10) is 1, the system (15.14) determines the values 
ui”, ws, except for a common multiplicative factor which remains arbitrary. 
From the manner in which »,(x) is determined it is then seen that this function 
is also fixed except for a constant multiplier. If, on the other hand, X, is of the 
index 2 relative to the boundary problem (13.10), the system (15.14) admits of 
solution by two linearly independent pairs of values 1, we. Each of these leads in 
the manner described to an associated solution v(x), and the two of these so ob- 
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tained are also linearly independent. It has thus been shown that adjoint bound- 
ary problems have the same characteristic values, and, moreover, that each 
such value has the same index relative to each of the two problems. 

For the characteristic solutions v,(x) of the boundary problem (15.9) the 
formula (15.7) yields the form 
(15.15) = ha” Vala, dn) + 
The coefficients in this, as they may be drawn from the equations (15.10), have 
the evaluations 


In general the boundary problems (13.10) and (15.9) are distinct. In a re- 
stricted class of cases, however, they may be effectively the same. Boundary 
problems of this class are said to be self-adjoint. Many familiar physical phe- 
nomena admit of mathematical formulations in terms of self-adjoint differential 
systems. 


CHAPTER 16 


Generalized orthogonality. Of the methods for the determination of coeffi- 
cients in a trigonometric representation, that of Euler (Chapter 8) was seen to 
be by all odds the simpler one. It is more direct and very much shorter than 
Fourier’s (Chapters 10, 11). And this advantage was recognized in Chapter 11 
to be attributable essentially to the fact that from the very start it exploits 
strategically the peculiar properties of the functions in terms of which the 
representation is made. These properties are particularly those which enter into 
the so-called orthogonality of the trigonometric functions, namely, those which 
come to their most familiar expression, at least in part, in the relations (12.1), 
(12.3), (8.10) etc. A directive influence upon the present discussion is evidently 
to be discerned in this fact. The réle that is filled by the trigonometric functions 
in the classical theory is to be assigned in this generalization to the characteristic 
solutions of the boundary problem (13.10). The discovery of interrelations be- 
tween these solutions, and especially of such as reduce to orthogonality under 
suitable specialization, therefore appears as a significant issue at this turn. 

In terms of the coefficients h”, hY’, with which the relation (14.10) main- 
tains, let the function U,(x, \) be defined by the formula 


This is evidently a solution of the differential equation (13.3), specifically one 
which becomes a characteristic solution of the boundary problem (13.10) when 
d is given the value A,. If, then, as usual, v,(x) is used to designate the mth char- 
acteristic solution of the adjoint problem, it is clear that 
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L(U,, = 0, 
M (vn, A) = (A — An) g(x) 


In virtue of these relations the identity (15.2), with U,(x, A) and Yn(x) i in the 
place of ¢ and y, yields, upon integration, the equation 


(16.2) q(x) U(x, A)on(x)dx = Q(Us, Mm, — Q(U;, Mn, 5). 


Now when p™, »” and X, stand in the place of 1, we and \, the boundary 


relations of the system (15.9) are fulfilled by the function v,(x), while their right- 
hand members may be written as 


wa” — (A — + [on — = 1,2, 3,4. 


Upon multiplying them respectively by U; (a, \), U,(a, ), (b, d), d), 
and then adding them, it is accordingly found that 


Q(Ur, — Q(Ur, 6) = ws + Aa(Ue, d) 


the symbols B; having been introduced here in the sense 


(16.3) = + + + $= 1, 2. 


The equation (16.2) therefore assumes, after a division by (A—A,), the form 


B.U,)], 


b 
(16.4) f q(x) A)on(x)dx + B(U;) + ui B(U,) = = 4,,,-(A), 


in which the right-hand member is explicitly given by the formula 


(n) 
d) + 


We are to be specifically concerned with the form of the relation (16.4) and hence 
with the value of the expression (16.5) when \=),. 
From the relations (16.1) it is to be seen at once that 


(16.6) A,(U,, = he + Aj,2(d), j = 1, 2. 
The equations (14.9) thus assure the relations 
(16.7) = 0, j= 1,2, 


and from these it follows that the expression (16.5) vanishes at \,, whenever 
\,¥X,. For its evaluation at \ =A, we may observe that by virtue of the formulas 
(16.6) the numerator of the expression (16.5) may be written out explicitly as 
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and that the relations (15.14) thereupon show it to reduce to zero at A,. At this 
value of \ the expression (16.5) is, therefore, indeterminate. To fulfill the rela- 
tion (16.4) it must, however, be continuous. The value to be assigned it at An 
is, therefore, that which is obtainable by an application of the familiar “l’Hospi- 
tal’s rule,” namely, with the use of a superscribing dot to denote a derivative 
with respect to A, thus 


a 
16.8 F ’ —F ’ ’ 
(16.8) (2, |= (x 

the value 


For the further analysis of this formula we must distinguish between the case in 
which the index of \, is 1, and that in which it is 2. 

If X, is of the index 1 we are still concerned with the formula (16.9) only in 
the case that r=n. In this instance the relation 


(16. 10) Ag,1(An) 0, 


is, moreover, fulfilled for some choice of the subscripts o, 7, and this implies 
through the equations (14.9) and (15.14) that 


(16.11) 0, #0. 


Now the relations (14.4) and (16.1) assure the equation 
Ai(Un, A2(Un, d) 
A;,,(A) 
and a differentiation of this, together with the evaluations (16.7), shows that 

A1(Un, Xn) A2(Un; An) 
| Aayr(An) 


The determinant on the left of this equality is, however, that of the system 
which is comprised of the equation (16.9) with r=m, and the equation 


= (— 1) kh, A(A), 


(n)> 


= (— 1) h, AQ,). 


which is one of the pair (15.14). The eliminant of the coefficient of A.,,(An) in 
this system is, therefore, found to be the relation 
(n) (n) 


t Me 


(16. 12) = (— A(An)- 


Ag,1(An) 


E 
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The right-hand member of this is different from zero or equal to zero according 
as the multiplicity of \, as a root of the characteristic equation (14.5) is equal 
to or greater than 1. It has been shown thus for the case at hand, namely when 
the index of X, is 1, that the relations 


16.13) if rn, 
#0, if r=n, 


maintain if and only if the multiplicity of \, and its index are equal. 

Consider now the case of a characteristic value A, of the index 2, and let its 
designations in the array (14.7) be Am and Anyi. The set of equations obtained 
from the relation (16.9) when m and r are given values from the pair m, m+1, 
yield the determinant relation 


Pmtt.m(An) Ai(Um An) Xn) 
Now the formulas (16.6) assure the equation 

A,(Um, d) Adta® ) (m+1) 


and a two-fold differentiation of this together with the evaluation (16.7) leads 
to the equality 


Ai(Um, An) A(Um, An) he he 
i cm) (m+) | A(An). 
mi, An) A2(U m41, An) hy hy 


The combination of this result with that above evidently leads to the relation 


If the multiplicity of \, is greater than 2 the right-hand member of the 
equation (16.14) is zero. It is clear that in that case the equation is contradicted 
by the relations (16.13) and that these latter, therefore, do not maintain. If, 
on the other hand, the multiplicity of \, is 2, namely the same as its index, no 
such contradiction is involved. If in this case the parameter pairs py, wy and 
uit’, wt’, which have thus far been left unspecified except that they be 
linearly independent, are now specified to fulfill the relations 


16.15 
A\(Un, Xn) + A(Um, Xn) = 0, 


it follows at once through the equations (16.9) that Bm m4i(An) = Pm4i.m(An) =0, 
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and then through the equation (16.14) that ®,,,(A,) #0, for n=m, m+1. The 
relations (16.13), which by virtue of the equation (16.4) take the form 


(n) (n) =0, if r¥n, 

#0, if 
thus maintain whenever the characteristic value \, is one whose index and 
multiplicity are equal, and do not maintain in any other case. These relations 
are the ones that were sought. They are, namely, those which are expressive of 
the generalized orthogonality which subsists among the characteristic solutions 
of two boundary problems that are adjoint. 


CHAPTER 17 


The formal representation of an arbitrary function. In the instance of any 
boundary problem which admits of infinitely many characteristic solutions, the 
physical context in connection with which the problem arises leads in a natural 
way to the question of the representability of an arbitrary function in terms of 
these solutions in the manner 


(17.1) fla) = > 


This has already been observed in several instances of the Fourier theory, which 
is, of course, exemplary of the more general case. As in the trigonometric case, 
the crux of the formal problem thrown up in this way devolves upon a deter- 
mination of the coefficients c,. We shall consider this matter now, not in a rigor- 
ous way, but formally. The relation (17.1) will, therefore, be taken to be amen- 
able to all such operations as shall be made upon it, and no consideration will be 
given to matters of convergence. The deductions will, therefore, of course, be 
devoid of all power of proof. Their purpose is purely an exploratory one. 

To begin with, let the symbols f;, i=1, 2 be used as abbreviations for the 
expressions B,(f), as these latter are obtainable from the relations (16.3). The 
equation (17.1) thus has associated with it the pair of auxiliary relations 


(17.2) rBi( ur), i= 1,2, 
r=0 
and in terms of the elements f(x), fi and fz the formulas 


(17.3) L(f) = f q(x) f(x)dx + fit uy n=0,1,2,-°-, 


may be taken to define their left-hand members. Upon substituting into these 
formulas the infinite series evaluations (17.1) and (17.2), interchanging the order 
of the integrations and summations, and collecting the terms in any coefficient 
¢r, it is found that alternatively 


r=0 


FORMAL REPRESENTATION OF AN ARBITRARY FUNCTION 59 


I,(f) = 


The symbols ®,,,(A,) in this stand, as heretofore, for the left-hand members of 
the relations (16.16). Since these latter vanish whenever rn, the equation 
reduces to the simple form 


(17.4) = CnPn,n(An)- 


Let it be assumed now asa hypothesis which is to cover the entire remaining 
portion of these deductions, that the boundary problem (13.10) in question is one 
for which there are infinitely many characteristic values, and for which, moreover, 
each of these is of a multiplicity that is equal to its index. By the relations (16.13) 
the equations (17.4) yield, then, for each n, the evaluation 


1,(f) 
(17.5) 


It will be evident even from the most casual review of the procedure de- 
scribed that it is limited, insofar as actual applicability is concerned, to the 
functions of a materially restricted class. The expressions B;(f) are, for instance, 
significant only for functions that are differentiable at x =a and x =}, while other 
heavy restrictions are manifestly involved. Beyond that it is clear that the result 
can imply nothing of the representability of a function f(x), since that repre- 
sentability was at the very outset assumed. A theory of representation must, 
accordingly, be approached differently. 

Let it be supposed, therefore, that a function f(x) which is arbitrary except 
that it possesses certain requisite properties of integrability, and a pair of con- 
stants f1, fe are given. The constants are likewise to be regarded as arbitrary. In 
particular they need have no specific relation to the values of f(x). From these 
elements f(x), fi, and fz the values J,,(f) are constructible through the formulas 
(17.3). There is thus associated with them a sequence of constants (17.5), or, 
in other words, a series of the form (17.1). We shall indicate this association in 
the manner 


(17.6) f(x) 


There is no implication at this stage that the series here written down is con- 
vergent, or, even should that be the case, that its value is f(x). The continuing 
theory is to be shaped toward the investigation of those matters. As it stands, the 
series is to be approached afresh and with no regard for the manner in which it 
was deduced. Its convergence is to be studied. Any conditions, if such are found, 
under which the series does converge and to the value f(x), will be conditions 
under which the symbol of association in the relation (17.6) may be replaced 
by one of equality and a theory of representation maintains. The direction of 
our further theoretical developments is thus forecast. 


r=0 7 
: 
= 
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CHAPTER 18 


Some examples. Throughout the discussion thus far no references whatso- 
ever to illustrative examples have been made. A stage has now been reached, 
however, at which the consideration of some explicit cases may be both interest- 
ing and instructive. We shall, therefore, interrupt the theoretical developments 
at this point to review some particular examples with the special purpose of 
illustrating concretely some of the salient features of the theory as it has thus 
far been deduced. It will be seen in connection with one of these cases how the 
present more general theory wholly includes that of Fourier. 

Because of the breadth of the basis upon which the theory is built, it is a sim- 
ple matter to construct in great variety examples that are subsumed under it. 
In the choice that is set forth below, however, an effort has been made to avoid 
in as large a measure as possible all such complications as are not germane to 
the matters essentially at issue. Such could be only distracting and not illuminat- 
ing. For this reason, in particular, the boundary problems that are drawn upon 
have all been chosen to be such as involve differential equations that are ex- 
plicitly solvable. The theory in no way requires that. For convenience in the 
notation the symbol p has been introduced to stand for a square root of \, thus 


(18.1) 


This, however, is to be looked upon merely as an abbreviation. The effective 
parameter will continue to be \ even when an expression is written in terms of p. 
It will generally be found on this account that a form that is indeterminate as 
to p is quite specific in A, as is the case with the form sin ap/p, which near \=0 
is to be thought of as defined by the power series 


It may also be worth observing in conjunction with the notation defined through 
the relation (16.8), that 


(18.2) F=-— 


Example 1. The boundary problem 
2 2 
ul’ — = 0, 
x 


(18.3) u(=) = 0, 
2 


a ad 
; 
: u = = 0, 


SOME EXAMPLES 61 


The fundamental interval at the ends of which the boundary relations of this 
problem apply, namely 7/2 $x <37/2, is one upon which the differential equa- 
tion has continuous coefficients and no singular points. Beyond that the coeffi- 
cients a;,; of the boundary relations are the elements of the matrix 
0100 
0001 
and this, the matrix (13.8), is of the rank 2 for all values of \. The boundary 


problem is thus of the type (13.10). Its differential equation admits as linearly 
independent solutions analytic in A, the functions 


o1(x, = + e-7], 


’ 


p 


and as it is formed from these the determinant (14.4) is 


A(A) = 
2 2 
The evaluation 
— 
2 p 
shows that its zeros occur at the points p= + nt, with n=1, 2,3, - - - . At each of 


these points the determinant A(A) is reduced to the rank 1 and A(A) is not zero. 
The characteristic values, which are infinitely many, are thus all of the index 
and multiplicity 1. They are to be arranged after the fashion (14.7) in the order 


= — (n+ 1), n=0,1,2,---. 


When r is any even integer, say r =2n, the coefficient of h2” in the first one 
of the equations (14.9) is zero. This is seen to determine 4°” to be zero, while 
leaving hf” arbitrary (not zero). This latter coefficient may evidently be chosen, 
therefore, so that the respective characteristic solution, as it is given by the 
formula (14.10), is explicitly 


thon(x) = x cos (2n + 1)x. 


In a quite similar manner it is found that 42"~” =0, and that we may accordingly 
choose 


Uon—1(x) = x sin 2nx. 


By the formulas (15.9) the differential system adjoint to (18.3) is 


‘ 
q 
: 
| 
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Z 
v’ +—v' —dv=0, 
x 


(=) 4 4 (=) 


Of the equations in this, the first, the second, and the fourth, effectively define 
the adjoint boundary problem, while the remaining ones give equivalents of the 
parameters pw and pe. The formulas (15.5) yield as solutions of the differential 
equation in this system the functions 


—1 
vi(x, A) = —— [er= — 
4px 


1 
\) = — + 
4x 


The equations (15.16) show that for any n, k?""'=0 and k3"=0. The character- 
istic solutions v,(x) may accordingly be taken to be 


sin 2nx 


Ven—1(*) = 


cos (2n + 1)x 
x 


Von(x) = 


The boundary problem (18.3) does not involve \ in its boundary relations. 
Hence the constants f;,; that appear in the formulas (13.7) and (16.3) are all 
zero, and the forms (16.3) accordingly all vanish. The relation (16.4) thus shows 
that 


= — f Uy(x)0_(x)dx, 


and the relations of orthogonality (16.16) therefore assume a familiar, purely 
trigonometrical form. With the special choice of constants f;=0, f2=0, the 
formulas (17.3) and (17.6) yield, in association with an arbitrary function f(x), 
the series of characteristic solutions 


f(x) ~ cox cos x + sin 2% + cox cos 3% + cgxsin4dx+---, 


ve + 
2 T 2 
3x 
| 
: x 
| 
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with the coefficients 


2 34/2 sin 2ns 
= — f f(s) ds, 


2 2 
f 
x/2 


Example 2. The boundary problem 


+ 
x 


Au = 0, 


(18.4) u(t) = 0, 


2u'(1) — u’(2) = 0. 


This boundary problem is of the form (13.10) on the fundamental interval 
1 <x <2. Its differential equation admits the solutions 


¢:(x, A) = cos (= 


1 2: 
¢2(x, \) = — sin (= 
p 


and as formed from these functions 


3 
—|sin{—]}-—sin —|cos{—}—cos {— 
a a 3 3 3 3 
The evaluation 


8r 
A(A) = rs (1 — cos 2zp), 


shows that the zeros occur at the points p=0 +1, +2,---. There are thus 
infinitely many characteristic values, and since at each of them the determinant 
is reduced to the rank 1 they are all of the index 1. We may accordingly write 


An = =0,1,2,---. 


But it is now seen at once that A(A,) =0 for each m. The characteristic values 
are thus each of a multiplicity that is higher than its index and the boundary 
problem is therefore excluded by the hypothesis of chapter 17. 


2p 1 : 
— sin | — 
« 
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Example 3. The boundary problem 
+ dru = 0, 
(18.5) u(— 3) — u(x) = 0, 
u'(— — = 0. 


On the interval —7Sx Sm this boundary problem is of the form (13.10). 
We find for it 


sin px 
(18.6) (x, A) = Cos px, $2(x, A) = 
and accordingly 
2 
0 — — sin pr 

(18.7) A(A) = p 

2p sin pr 0 
Since from this 
(18.8) A(A) = 4 sin? pz, 
its zeros occur for p=0, +1, +2,---. At the first of these the rank of the 


determinant (18.7) is 1, at all others it is 0. The first characteristic value \) =0 
is thus of the index 1 and the remaining ones are of the index 2. They are ac- 
cordingly to be arranged thus 


(18.9) ho = 0, Non—1 = Arn = m= 1,2,3,--°. 


Each value is found to be of a multiplicity equal to its index. By the relations 
(14.9) it is shown that h© =0, so that the first characteristic solution is uo(x) =1, 
or any multiple of this. For the other characteristic values the equations (14.9) 
are vacuous and the coefficients h* and h¥ accordingly arbitrary. We may 
therefore choose and as Aen—1) and g2(x, Aen) respectively, 


U2n—1(x) = COs nx, 
(18.10) 
tho,(x) = sin nx. 


The boundary problem (18.5) is readily seen to be self-adjoint. It is, there- 
fore, admissible to choose v,(x) =un(x) for every n. Since the boundary relations 
are independent of \ the forms B,(u,) all vanish, and it is found accordingly that 


= f 


The relations (16.13) are thus, in this instance, merely expressive of the familiar 
property of orthogonality of the sines and cosines of multiples of x. With the 
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choice of constants f:=0, fe=0, the formulas (17.3) and (17.5) define the co- 
efficients ¢, to be those of Fourier, and the relation (17.6) is simply the associa- 
tion of a function f(x) with its Fourier series. 


CHAPTER 19 


Another example. Although the examples that have been cited in the chapter 
above were earmarked among themselves by some conspicuous qualitative dis- 
similarities, they fail, even when taken together, to illustrate adequately some 
of the more pronounced departures from classical theory which the generaliza- 
tion permits. The reason for this lies in the fact that the boundary problems 
that are basic to them all involve only such boundary relations as do not depend 
upon the parameter \. In the following a problem of the excepted class is to be 
briefly analyzed. 


Example 4. The boundary problem 
+ 4ru = 0, 
(19.1) u(0) = 0, 
u’(1) + ++ 1—A]u(1) = 0, 


with vy standing for any real constant. 


Such boundary problems as this present themselves in connection with a 
variety of physical investigations of which the following ones may be looked 
upon as in some measure typical. [18] 


Problem (i). A right cylindrical solid with a cross-section of any shape and 
size, and with plane terminal faces at x =0 and x =1, has its lateral surface insu- 
lated against the passage of heat and has an initial distribution of temperatures 
depending only upon the longitudinal coérdinate x. At the time t=0 this solid is 
placed in contact with a quantity of liquid at one of its terminal faces, and the 
liquid is thereupon kept well stirred to insure that the temperature is uniform 
throughout it at each instant. The temperature of the liquid and the temperature 
distribution in the solid at any subsequent time are to be calculated. 


Problem (ti). A solid metal sphere with an initial distribution of temperatures 
that is symmetrical about its center, is cooled by being plunged into a mass of 
liquid. The liquid is kept well stirred. The temperatures of the liquid and those 
within the solid during the cooling are to be determined. 


Problem (iii). A mass of material is uniformly distributed at the time ¢=0 
throughout a jelly in a cylindrical container. The jelly is covered with a liquid 
that is kept well stirred. From the experimental measurements of the concentra- 
tion of material in the liquid as a function of the time, the coefficient of diffusion 
of the material in the jelly is to be found. 


£ 
| 
is 


2 
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The boundary problem (19.1) is of the form (13.10). Its differential equation 
admits as solutions the functions 


X) = cos 2px, 


(19.2) sin 2px 


and the evaluation of the determinant 


1 0 


19.3) AQ) = in 2p |, 
( —2p sin 2p+(v?+1—A) cos 2p 2 cos 2p+(v?+1—A) 


gives as the characteristic equation 


sin 2 
Alternative forms of this are 
p 
and 


The characteristic values are thus the squares of the roots of a transcendental 
equation and it is readily seen that they do not admit of expression by any 
elementary formula. The essential facts concerning them are nevertheless deduci- 
ble. 

Thus if p is any complex value with a positive imaginary part, the point 
(p+7) in the complex plane is more distant than the point (p—7) from any point 
on the axis of reals, and hence in particular from the points y and —v. Hence 


and from these relations it follows that the right-hand member of the equation 
(19.4c) is greater than 1 in absolute value. But for a value of p such as this the 
left-hand member of that equation has an absolute value that is less than 1. 
This p is, therefore, not one for which the equation is satisfied. A similar argu- 
ment establishes that same fact for any value of p that has a negative imaginary 
part. The roots of the equation must thus all be real. 


2 
an 
Pp 
0 
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Now that being so, it is observable from the equation in the form (19.4b), 
that its roots are the abscissas of the points common to the graphs 


y = cot 2p, 
p ve+1 


in the Cartesian (p, y) plane. These graphs are easily seen to intersect infinitely 
often and without being tangent to each other at any intersection point. From 
this it follows at once that there are infinitely many characteristic values and 
that they are each of the multiplicity 1. From the general relation (14.6), or from 
the fact that the determinant (19.3) has a constant non-vanishing element, it is 
to be seen that every characteristic value is of the index 1. 

By the first one of the equations (14.9) it is shown that for each r, hf =0. 
The characteristic solutions may, therefore, be taken as the functions 


(19.5) un(x) = sin x, 


and since the boundary problem is self-adjoint we may also take v,(x) =u_(x). 
The formulas (16.3) and the boundary relations of the system (15.9) give in this 
case the evaluations 


B,(u,) = 0, Bo(ur) 


(n) (n) 


— u,(1), 
— v,(1). 


It is thus found from the equations (16.4) that 
1 
0 


Because of the simple form of the functions u,(x) the integration in this expres- 
sion is explicitly possible, and the relations of generalized orthogonality (16.16) 
are therefore verifiable upon the basis of the equation (19.4a). The formulas 
(17.3) and (17.5) yield as the coefficients in the association 


(19.6) Cy sin 
n=0 
the values 


af f(s)un(s)ds — fiug (0) — fomn(1) 
(19.7) 


4 f ree + u,(1) 
0 


} 
2 
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A scrutiny of the formulas (19.7) reveals one of the salient features in which 
this and the Fourier representations are effectively dissimilar. While the coeffi- 
cients of the latter are wholly determined by the function f(x), those given by the 
formulas (19.7) depend also upon the prescribed constants f; and fe. The same 
function f(x) may thus be associated here with many different representations. 
This is quite consonant with the nature of the physical problems from which 
boundary problems of the type (19.1) arise, as may easily be appreciated from 
a consideration of the problems (i) and (ii) that were cited above. In each of those 
cases the function f(x) represents the initial temperatures in the solid. The sub- 
sequent temperatures therein depend, however, not only upon these but also 
upon the initial temperature of the fluid into which the solid is immersed or with 
which it is placed in contact. It is through the constants f; and f. that this initial 
fluid temperature comes to account. [19] 


CHAPTER 20 


The Green’s function. A noteworthy feature of the Fourier’s series, and one 
which is almost invariably taken as the point of departure for studies of its 
convergence properties, is the fact that any initial segment of it may be ex- 
plicitly summed and hence expressed by a compact formula. This was observed 
in chapter 12, the summation of the segment Sy(x) as given by the relation 
(12.7) having been accomplished there by the formula (12.8). It is far from 
obvious that the advantages inherent in this are retainable in the generalization 
of the theory, for the derivation of the formula in question is directly and ex- 
plicitly based upon relationships that are peculiarly trigonometric. We shall 
show that this may nevertheless be done, the key to the requisite deductions 
residing in a certain function, the so-called Green’s function, which generally 
plays an important réle in the theory of boundary problems. 

Let \ be taken and retained throughout the deductions of this chapter as not 
a characteristic value of the boundary problem (13.10). This problem then ad- 
mits of no solution, which is to say that a function which fulfills its boundary 
relations cannot also satisfy the stipulations of its differential equation for all 
values of x upon the fundamental interval. Now it is significant that the conces- 
sion which must be made on the part of the differential equation is very slight, 
amounting to no more, in fact, than the partial relaxation of its stipulations at 
only a single point x =s of the interval. The function which fulfills the specifica- 
tions except at x =s may even be required to be continuous there. The failure 
will occur in its first derivative, which is subject at this point to an ordinary 
discontinuity. When this discontinuity is of the proper sign and of the unit 
magnitude, matters which are adjustable without otherwise affecting the issues, 
the function is called the Green’s function. Inasmuch as it depends upon the 
location of the point s as well as upon the variables x and X, it is to be denoted by 
the symbol G(x, s, A). Its properties, by way of summary, are, then: 
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(i) As a function of x it satisfies the differential equation of the system 
(13.10) over each of the intervals a<x<s, and s<x<b; 

(ii) It is continuous in x over the whole interval (a, 5), but its derivative 
is discontinuous at the point x=s to accord with the prescription 


0 
— —G(x, s, »)| = 1; 
Ox 


r=s— 


(20.1) 5, 
Ox 


z=8+ 


(iii) It fulfills the boundary relations of the differential system (13.10). 

By these properties the function G(x, s, \) is uniquely determined, as the fol- 
lowing derivation of its form will incidentally show. 

As a solution of the differential equation (13.3) the function is expressible 
upon each of the intervals (a, s) and (s, b) as a linear combination of the solutions 
$;(x, d), j=1, 2, namely 
¥1,261(%, + d), for a Ss, 


G(x, s,) = 
4) + A), for sd. 


The specifications (ii) impose upon these forms the relations 


2,21(S, A) + y2,162(s, A) = Y1,261(s, A) + ¥1,1¢2(5, A), 
2,201 (S, A) + ¥2,162 (S, A) — Y1,267 (8, A) — (s, A) = 1, 


in other words, a pair of equations which may be solved into the form 


o1(s, d) 

¢$2(s, d) 


Because of the formulas (15.5) these relations are alternatively 
V1.1 = — ¥2(s, A), 
71,2 = — ¥als, d). 
The formulas for G(x, s, \) are, therefore, more explicitly 
(20.2) G(x, 0) = + + 5, 9), 
in which the final term is that defined by the relations 


— o1(%, A)¥als, A) — a(x, A), for 


20.3 »5,A) = 
( )  gilx ) { 0, for « 6. 


4; 
7 
| 
7 
j 
7 
| : 
‘ 
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The requirement (iii), that G(x, s, A) fulfill the boundary relations, comes to 
its expression in terms of the symbols (13.9), (14.3) and (15.11), in the pair of 
equations 


+ 72,4520) — Air Als, — 4) = 0, = 1,2. 


Their solution determines for the coefficients 2,2 and 2, the evaluations 


i+1 (a) (a) 


(- 1) (A)yr(s, d) + A1,2 (A)v2(s, d) 1 2 
| Aas) ASLO») + A520) 
The substitution of these into the formula (20.2) leads to the result 
g2(x, 5, d) 
(20.4) G(x, 5, d) gi(x, 5, d) + AQ) 
in which 
A) A) 0 
(20.5) go(x, s, 4) = —| Arald) Arr A) + Ara a) |. 


The Green’s function has thus been completely determined. 

The method which has in this way been described in connection with the 
differential system (13.10), may now be applied equally well to the calculation 
of the Green’s function G,(x, s, \) of the adjoint system (15.9). If this is done it 
will be found that the two Green’s functions are simply related by the equation 


(20.6) G.(x, s, 4) = G(s, x, d). 


The function (20.4) thus operates as the Green’s function of the adjoint differen- 
tial system when its second argument is taken to be the variable and the first 
one is fixed. The set of its properties enumerated above may thus be extended to 
include the following ones: 


(iv) As a function of s it satisfies the differential equation of the system 
(15.9) over each of the intervals a<s<x and x<s<b; 

(v) It is continuous in s over the whole interval (a, 6), but its derivative as to 
s is discontinuous at the point s=x to accord with the prescription 


0 
(20.7) - »)| = 1; 


s= r+ 


(vi) As a function of s it fulfills the boundary relations of the differential 
system (15.9), with appropriate determinations of u(x, A) and pe(x, A). 


eae 
~ 
} 
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It is on the whole simpler to verify these facts than to deduce them. By the 
formulas (20.4), (20.3) and (20.5), the function G(x, s, A) is evidently expressed 
as a linear combination of the solutionsy,(s, A), 7=1, 2, on each of the intervals 
(a, x) and (x, b). The property (iv) is thus assured to it. The forms 

o1(x, d) + $2(x, d), 
o1(x, (s, d) + $2(x, (s, d), 
are readily seen from the relations (15.5) and (13.5) to take on the values 0 and 1 


respectively at the point s=x. By that the properties (v) are implied. Finally the 
formulas (15.5) and (15.11) may be drawn upon to supply the evaluations 


By virtue of them it is clear that 
g2(x,A) 0 
G(x, a, = AiaQd) 41,20) 
Az,2(A)  @2,1(A) 


This, however, is simply the statement that 
G(x, a, d) d)ar,1(A) + A)a@e,1(A), 


namely that the first boundary relation of the system (15.9) is fulfilled, with the 
parameter values 
(—1)* | A) A) 


By the same procedure it may equally well be shown that the remaining bound- 
ary relations are also fulfilled with the same values (20.8). Thus the property 
(vi) maintains. 

Although it will not be relevant to the discussion which follows, a conspicu- 
ous property of the Green’s function may still be mentioned. This is, namely, 
that the non-homogeneous differential system 


L(w, d) = F(x), 
Ai(w, = 0, 
A2(w, A) = 0, 


which is obviously related to the system (13.10), is solved, whatever the function 
F(x) may be, by the formula 


w(x, A) = f “G(x, Ss, A)F(s)ds. 


= 
l 
| 
| 
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CHAPTER 21 


The residues of the Green’s function. If the requirement that \ be fixed is 
now abandoned and this parameter is, on the contrary, regarded henceforth as a 
variable that is free to range over the complex plane; it is recognizable from the 
formulas (20.3) and (20.5) that the functions gi(x, s, A) and go(x, s, A) are analytic 
in \. The same is also true of the function A(A). The formula (20.4) therefore 
shows that the Green’s function is analytic except for poles, these latter occur- 
ring at and only at the zeros of A(A), namely at the characteristic values. It has 
been assumed as a part of our general hypothesis that each characteristic value 
is of a multiplicity equal to its index. As a consequence of that it will be found 
that for each m the product (A—A,)G(x, s, \) approaches a finite limit when 
A—A,. This limit is known as the residue of G(x, s, \) at An. We shall use the 
prefix “res,” to designate it, thus 


(21.1) res, G(x, s, A) = lim (A — A,n)G(x, s, A), 


and shall show that certain residues closely related to these are significantly 
associated with the terms of the series in a representation (17.6). 

Since the function g;(x, s, A) is analytic in \ when x and s have any specified 
values, this part of G(x, s,A) evidently makes no contribution to the right-hand 
member of the relation (21.1). It is, therefore, only the part go(x, s, A)/A(A) that 
needs consideration. Upon expanding the determinant (20.5) by the elements 
of its last column, therefore, it is found that the relation (21.1) may be more 
explicitly written as 


Amr, i=1 


the functions E;(x, \) being defined by the formulas 


(—1)' | d)  2(x, A) 
A(A) |  Ai,2(A) 


(21.2) ) = 


From the alternative form 
res, G(x, s, A) 


(21 . 3) (a) 


= {res, Ex(x, Xn) + As Xn) |, 


t=1 


it thus follows that the evaluation of the residues of G(x, s, 4) may be made to 
depend upon the determination of residues of the functions (21.2). To make 
these determinations we must consider separately the case in which the index of 
\, is 1 and that in which it is 2. In doing that it will be convenient to employ 
the symbol 6;,; to stand for the “Kronecker delta,” namely in the sense 
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1, if ¢=j, 
0, if 
Suppose then, to begin with, that d, is of the index 1. Then 


AQ) 
lim 


rm—r, A— An 


A(An); 


and there is a pair of subscripts ¢, 7, for which the relations (16.10) and (16.11) 
maintain. It will be clear that we may write 


gi(x, A) d2(x, 0 
(21.4) (A — Aa) 4) = ——— | Arad) Ara(d) brea], = 1,2, 


A(a) 
A22(A) 
and that therefore 
$1(%, An) G2(x, An) 
A4,2(Xn) 
Aai(An) A2,2(An) 


A(An) 


(21.5) res, E;(x, \) = 


This formula may be considerably modified. It will be seen that if it is multiplied 
on the right by the determinant 


(n) 


O 
0 ’ 
0 0 1 


and this operation is then compensated for by dividing out the value of this 
determinant, (—1)th™, the resulting form is, by virtue of the relations (14.9) 
and (14.10), 


Un(x) rn) 
O 
Aas(An) 


-1 


The further multiplication of this on the left by the determinant 


0 0 

0 51,0 52,6 ’ 
OM 


and the appropriate subsequent division gives it the form 
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Un(x) (x, An) 
0 


0 0 


1 
(- 


by virtue of the relations (15.14). The expansion of this reduces it, because of the 
equation (16.12), to the evaluation 


and the substitution of these results into the equation (21.3), together with an 
application of the relations (15.16) and (15.15), leads to the conclusion that 


Pain(An) 


Finally with the formulas (21.6) and (21.7) at hand it is a simple matter to recog- 
nize that the value of 


(21.6) res, E;(x, A) = i= 1,2, 


(21.7) res, G(x, s, 4) = 


b 
(21.8) TeSp {f q(s)G(x, s, A)f(s)ds + Ex(x, A) fe + Ex(x, 


is precisely the term 


Tn(f)Un(x) 


of the series (17.6). 


If \,, is of the index 2 and appears in the array (14.7) as Am and An4i, the rela- 
tion (21.4) is more appropriately taken in the form 


$1(x, A) A) 0 
Ai,(A)  A1,2(A) 
 Az,2(A) 
53-i,2 
A—An A—An 
(A — An) Ei(x, A) = 0 0 i 
 A1,2(A) 0 
A—An A—An 
Ae,2(d) 0 


Since in this instance each element A;,;(A) vanishes at Aq, it is seen that 


: 
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Asi) _ 


and hence that 

$1(%, An) An) 0 

—| Ar2(An) 

Aza(An) 
0 0 1 

Ai(\n) 0 


res, E,(x, A) = 


Let each one of the determinants in this ratio be multiplied on the right by 


0 
0 oO 1 


Because of the relations (14.10) and (16.6) it then takes the form 
bm (x) Um+1( x) 0 
— | Ai(Um, An) An) 
Ax(Um Xn) A2(U An) 
0 0 1 
Ax(Um, An) 0 
Ax(Um, Xn) 0 


and if this is again modified by multiplying each of the determinants on the left 
by the factor 


1 0 0 
(m+1) (m+1) 

O pe 
(m) (m) 

pe Mi 


the result, as a consequence of the relations (16.9) and (16.13), is 


Um(x) Um+1(x) 0 
0 
0 
0 0 1 
0 0 
0 0 


lim 
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Upon expansion this yields the formula 


(m) (m+1) 
Um (X) Um+i(X) 


| 


in accordance with which the relation (21.4) assumes the explicit form 


(21.9) res, E,(x, \) = 


a= 1,2, 


Um(x)Vm(S) 
Pm,m(Am) Pint, m+1(Am+1) 


Thus when A, is of the index 2 the residue (21.8) is the sum 


(21 10) G(x, d) 


of the pair of terms which appear in the series (17.6) in association with \,. 

From these conclusions it may now be seen that if the points representing 
the characteristic values are plotted in the complex plane, and if the first 
(N+1) of them taken in the order of increasing distance from the origin are 
enclosed by a curve Cy, the segment Sy(x) of the series (17.6) that is made up 
of terms that are associated with these characteristic values is given by the 
formula 


N b 
Sy(x) = TeSy, {f q(s)G(x, s, A)f(s)ds + Ex(x, A) fe + E2(x, 


Such a sum of residues is, however, familiarly expressible as a contour integral 
with respect to A over the curve Cy, namely as 


1 b 
(21.11) Sy(x) = {f q(s)G(x, s, A)f(s)ds + Ei(x, A) fe + E2(x, ay. 


This is the formula in the general theory that stands in the place of the Fourier 
formula (12.8) 


CHAPTER 22 


The Fourier’s series again. The formula (21.11) is useful, like its specialized 
counterpart (12.8), as the natural medium through which an investigation of the 
‘associated representations of arbitrary functions may be made. As the contour 
of integration Cy is successively taken to include more and more characteristic 
values the formula sums more and more terms of the series, and the approach 
of the integral to a limit reflects the convergence of the series. An identification 
of the function (if any) which the series represents thus becomes accessible 
through a study of the formula’s integrand, more particularly through an 
analysis of the Green’s function and the functions (21.2) when |A| is large. 


a 
& 
‘ 


1g 
st 
re 


Ip 
1e 


al 


er 


THE FOURIER’S SERIES AGAIN 77 


The exposition of the complete analysis that would be requisite for the gen- 
eral case would at this point go well beyond the bounds which have been set for 
the scope of this discussion. The character of the functions ¢1(x, A), ¢2(x, A) as 
these depend upon A would need to be determined. In general, however, no ex- 
plicit formulas for these functions are available, since the differential equation 
(13.3) is not ordinarily solvable. The difficulties which this circumstance inter- 
poses will be manifest, although, to be sure, they are not insurmountable. So- 
called asymptotic forms for the solutions of an equation of the type (13.3) are 
deducible by well established methods, and these are entirely adequate to the 
requirements. This theory of asymptotic solutions, however, we do not propose 
to go into here. Without it, a restriction of the discussion to considerably nar- 
rower confines than have hitherto been observed is in the nature of things in- 
evitable. We shall yield to the necessity by limiting the exposition henceforward 
to the basis of an outright, albeit a wholly typical, specialization. In fact, there- 
fore, the further investigation is to take the form of an extended analysis of the 
example 3 of chapter 18. It will be recalled that the theory of the boundary 
problem basic to that example is none other than the theory of Fourier’s series. 

For the boundary problem (18.5), and with the choice of constants f;=0, 
f2=0, the special form of the formula (21.11) to which the attention is to be given 
is 


(22.1) Sy(x) = —{ fa, s, A) f(s)dsdx. 
Cy 


The functions ¢$:(x, \), 2(x, 4) may be chosen as those that are given by the 
relations (18.6), and these lead through the formulas (15.11) and (15.5) to the 


evaluations 


@ a sin 
= cospr, AQ) = -— 


Azi(A) =psinpx, A3.2(A) = cos pz, 
and 


vi(s, A) = — sin ps, ¥2(s, A) = Cos ps. 
p 


The formula (20.5) is thus explicitly 


1 
COS px — sin px 0 
p 
(x, 5, d) = 2 1 
0 ——sinpr ——sinp(x+s)], 
p p 
2p sin pr 0 cos p(x + s) 


and this with the relations (20.3) and (20.4) yields the form 
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sin p(x — s) 
cos p(x + s — x) when —rSxSs, 
(22.2) G(x, = 
2p sin pr 
0, when 


In the complex \-plane the characteristic values are located at the points 
n’, with n=0, 1, 2, 3, - - -. The contour Cy, which must enclose precisely the 
first (V+1) of these, may, therefore, be chosen as the circle centered at the or- 
igin and of the radius (V+4)*. In the p-plane the semi-circle I'y which is 
centered at p=0, which has the radius (V+4), and upon which 0 Sarg p<z, 
is an image of the circle Cy under the mapping defined by the relation (18.1). The 
formula (22.1) with its integration expressed in terms of p is, therefore, 


1 
(22.3) Sy(x) = s, A) f(s)dsdp. 


When the point x at which the sum Sy(x) is to be considered lies in the in- 
terior of the interval (—7, 7), it is advantageous to take the formula (22.2) in 
the equivalent form 


cos p(s — x) 
(22.4) G(x, s, A) = 4 
2p 2p sin pr 


the symbol |s—x| standing as usual for (s—x) or (x—s) according as x Ss or 
x>s. The formula (22.3) then assumes the form 


(22.5)  Sw(a) -= f + f 


with the function V(s, x, N) defined by the formula 


1 e*°* cos p(s — x) 
(22.6) V(s, x, N) = — dp. 


sin pr 


Since this latter function is explicitly integrable as to s, it is found at once that 
for any choice of (a, 8) as a sub-interval of the range (—7, 7) 


W(s, x, N)ds = eae 08 ) dp 
iJr sin 
(22.7) 
1 e*** sin p(a — x) 
dp. 
p sin pr 


The integrals which appear in this formula closely resemble that in the formula 
(22.6) since the sine and cosine functions are effectively similar in structure. 
The analysis of the one, which is to be set forth, will therefore be found obvi- 
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ously and readily adaptable to the others also. In the event that the point x is 
an end point of the interval (—7, 7) a somewhat different grouping of terms in 
the formulas (22.4) and (22.5) is advantageous. The analysis is, however, es- 
sentially similar to that which applies when —1<x<z, and that being so we 
shall content ourselves with the discussion of this latter case alone. 

When it is expressed entirely in terms of exponentials the formula (22.6) may 
be made to appear in the form 


(22.8) W(s, x, N) = 


—1 4. gl2x—(s—z) pi 
—f dp. 


— 


Even relatively crude appraisals of the functions which appear in the integrand 
of this yield some significant results. Under the resolution p=£+in, with £ and 
n real, the equation of the arc I'y is 


& + 9? = (N + 4)’, n 2 0, 
and from this it is readily seen that each one of its points is either one for which 
or else one for which 
n> iV3. 


At every point of the first one of these categories cos 2r£<0, and the value 
which is explicitly 


(22.9) {1 — cos + 1/2, 


thus clearly exceeds 1. At any point of the second category the value (22.9) ex- 
ceeds (1—e-*" 4/4), and is thus ipso facto greater than 1/r. The relation 


(22.10) w|1— > 1, 


thus maintains over the whole arc I'y. 
Consider now the values 


(22.11) | | 
with s and x both upon the interval (—7, 7) and p still upon I'y. For all such 
s and x the relation 
is easily verified, and since the formula p= (N+ })e*, in which @ stands for arg p, 


shows that the real part of pi is —(N+4) sin 9, it is seen that the values (22.11) 
are both less than 
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Since this latter value is increased by the substitution of 3@ for sin 6 when 
0<0<S7/2 and of 3(7—8), for sin @ when 7/2 S07, we may evidently draw 
from the formula (22.8) the relation 


wis, f 


2 
+ 


x/2 
By explicit integrations this reduces to the inequality 


(22.12) | W(s, x, N)| < 
The companion result 

8 


may be drawn by similar reasoning from the formula (22.7). 

These conclusions are significant. In accordance with them the function 
W(s, x, N) is uniformly bounded as to s and N, and is such that its integral as 
to s converges to zero with 1/N uniformly as to the interval of integration. 
These properties, however, are precisely those which were invoked in chapter 12 
as being sufficient to insure the relation (12.10) for an arbitrary integrable 
function F(s). In the present instance, therefore, we may similarly conclude 
upon the basis of them that the final integral in the equation (22.5) converges 
to zero. Since the remaining integral in that equation may be evaluated thus 


1 in (N + 


< 


8 
f W(s, x, N)ds 


a 


(22. 13) 


it is clear that the equation (22.5) implies for every integrable function f(x) the 


relation 
sin(N +4)(s— 2) 


This is the relation (12.11) already familiar, and through it the reference of the 
sum Sy(x) to the Dirichlet integral has evidently again been accomplished. The 
method of its accomplishment here, however, by contrast with that of chapter 
12, is one of very general applicability. In particular it is one which in no way 
depends upon the special trigonometric combination formulas. With adaptations 
which essentially concern only details, this method is adequate to the analysis 
of the general series (17.6), namely to the representations of arbitrary functions 
in terms of the characteristic solutions of any properly constructed boundary 
problem [20]. 
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I. The solution of the system of equations 


= 0, 
(I. 1) + + = 0, k=1, 1), 
= 0. 


If, in terms of unspecified constants @ and £, a solution of the system is 
sought in the form 


(I. 2) up = ak — BF, 
it is found upon substitution into the equations to be requisite that 
(1.3) + ga + 1] — [67+ +1] =0, 1). 


T hese conditions are fulfilled if a and 8 are roots of the equation 

namely if 

(1.4) ap = 1, 


Since 8 must thus be the reciprocal of a, whereas the form (1.2) must vanish 
when k=7, it is seen to be necessary that 


a®*—a"*= 0, 
namely that 
a= B = 


with an integral value of v. For such an index v the relations (1.4) and (1.2) 
show that g and 1% have the values 


vr 
(1.5) qo = — + — 2 cos—) 
n 
and 
n 


respectively. Since the system under consideration is homogeneous, any multiple 
of a solution is also such. Hence we may write 
(I.6) Uy,~ = Ay Sin—)» k=0,1,2,°°+,%, 
n 


with the coefficient A, arbitrary. 
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For the indices y=1, 2, - - - , (w—1), the characteristic values of g, as given 
by the formula (1.5), are distinct. There can be no others, for by inspection the 
determinant of the system of equations is a polynomial of the degree (n—1) 
in g. 

II. A proof of the relation 
sin v§ = sin p,_1(cos &), 


with p,_; a polynomial of the degree (v—1). 
If in Demoivre’s formula 


(cos vé + 7 sin vt) = (cos — + 7 sin £)’, 


the right-hand member is expanded by the binomial theorem, each resulting 
term that involves 7 to an odd power also involves sin £ to such a power. Upon 
equating the pure imaginary components on the two sides of the equation it is 
thus found that 


sin vé = sin £-Q(cos é, sin? 


each term of Q being of the degree (v—1) in cos £ and sin £, and of even degree 
in sin £. By the substitution of (1—cos?£) for sin?é the form (II.1) results. 


III. A deduction of the identity 


1 j 
(IIT. 1) cosix =—— )> ( ) cos j — 2p)x. 
If in the familiar equality 
(IIT. 2) cos x = 3(e** + e-*), 


each member is raised to the jth power and the one on the right is then ex- 
panded by the binomial theorem, the result is the relation 


1 
(III. 3) cos? x = — 
2) 


u=0 


in which the symbol f ) designates the coefficient of a* in the expansion of 


(1+a)i. Since( -( ), the formula (III.3) may be written alternatively as 


1 j 
cos? x = — ( ) + e~ iz} 


2 i 


with [j/2] denoting the largest integer not exceeding j/2. By invoking the rela- 
tion (III.2) again this may be given the form (III.1). 
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IV. A derivation of the evaluations 


{" —4(1+ cossm), if s=0 (mod 2n), 
a cossm), if s #0 (mod 2n). 


p=1 n 


The first one of these evaluations is obvious. For when s=0 (mod 2”) each 
term in the sum on the left has the value 1 while the right-hand member is 
(n—1). If s#0 (mod. 2m) the relation 


usr 
n 2 


leads at once to the equation 


1 
n 2 tl 
or, as it may equally well be written, to 
n—1 psa n—1 
(IV. 2) cos = 4. 
n p=—n+1 


The sum in the right-hand member of this is a geometric progression. It is 
summed, therefore, by the formula 


= - 
i- 


and since 
(ettt/n)tm = (— 1)* = cos sz, 


its value reduces to —cos sr. The formula (IV.2) thus takes on the form of the 
second evaluation (IV.1). 


V. An evaluation of the determinant 


-~-1) 

(V.1) D(x) =| 15 

42-1 22-1 eee (v 


If in this determinant the elements of each row, beginning with the last 
one and extending in turn back to the second one, are modified by subtracting 
from them x? times the corresponding elements of the preceding row, the value 
of the determinant is unchanged. In its new form, however, all but the first ele- 
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ment in the last column are zeros, and hence, upon an expansion by the ele- 
ments of this column, it is found that 


1(12— 2(22— 2%) — 1]([» — 1]? — 2?) 
13(12 — x?) 23(27 — x?) 1]*([» 1]? — x?) 


1-3(12 — 2%-3(2? — x2). - [yp — — 1]? — x?) 


D(x) = (— 


As it now stands, the factor (n?—x*) is common to the elements of the mth 


column for n=1, 2, - - - , (v—1). Upon factoring these from the determinant the 
evaluation 

v—1 
(V.2) D(x) = (- TI - Din 

n=1 


is obtained, D,,, designating the cofactor of the element in the vth row and 
vth column of the original form of D(x). 


VI. A formal deduction of the relation 


(VI. 1) 


n=1 n 


The roots of the equation 2*=1 are obviously z=+1 and z=et"**/*, for 
n=1, 2,---+, (k—1). The factorization 


— 1 = — {(z — — e-ntilk)}, 
is, therefore, proper. It is, however, clearly equivalent to 
st) (gt — 2200s"), 
n=1 
a relation from which it follows that 


gt — gk k-1 n 
2% k 


2% n=1 


Now if in this the quantity e‘**/* is substituted for z, the equation becomes 


sin rx = sin —J[]]| 2 cos— — 2cos—]; 
or, since 
cos —— — cos—— = 2{ sin? —— — sin? ——}, 
k k - 2k 2k 


| 
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k k-1 mn) 
sin 


1 
1 in? ——> , 
k I 


n=1 


an evaluation by virtue of which the relation (VI.2) may itself be reduced to the 


form 
' sin { —— sin — 
gin rx k 2k 
sin —— 
k 2k 


This is valid for all k, and we may, therefore, permit this index to become in- 


finite. The evaluation (VI.1) formally results. 
VII. Establishment of the formula 


, [2n+1 
sin] 
(VII. 1) + COs sin [30] 


Upon substitution of the relations 


2 cos vO = + 


into the left-hand member of the formula (VII.1), this latter is found to be 
expressible as 


n 


er 


y=—n 


e(ntl) bi 


(VII.1). This latter is, therefore, established. 


This is an identity. Its limiting form as x0 must, therefore, maintain, namely 


for (mod 2z). 


This is a geometric progression whose sum, if 840 (mod 27) is 


If in this the numerator and the denominator are each divided by the factor 
— 2ie*!?, the fraction assumes the form of the right-hand member of the formula 
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extended than any other edition, will continue to meet a need in this country and 
elsewhere. 
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